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The  objectives  of  this  dissertation  are  (1)  to  conduct  a  detailed 
analysis  of  the  grown-in  defects  and  radiation  induced  defects  in  GaAs 
grown  by  the  LEC,  VPE,  LPE,  and  MOCVD  techniques  under  different 
growth  and  annealing  conditions,  (2)  to  identify  the  physical  origins 
of  the  deep-level  trap  in  GaAs,  and  (3)  to  determine  the  potential 
well  of  electron  traps  from  analyzing  the  electric  field  enhanced 
emission  rates  deduced  from  the  nonexponential  DLTS  data. 

A  detailed  theoretical  and  experimental  study  of  the  grown-in 
defects  and  radiation- induced  defects  in  GaAs  has  been  carried  out  in 
this  dissertation;  the  main  research  accomplishments  derived  from  this 
research  are  summarized  as  follows: 

(1)  Theoretical  and  experimental  studies  of  native  point  defects 
in  GaAs  grown  by  the  LEC,  VPE,  LPE,  and  MOCVD  techniques  under 
different  growing  and  annealing  conditions  have  been  made  in  this 
research,  and  the  conclusions  are  listed  as  follows:  (a)  high  purity 
GaAs  material  can  be  grown  for  the  low  arsenic  pressure  case  under 
optimum  cooling  condition,  (b)  GaAs  grown  under  a  high  arsenic 
pressure  condition   will   produce  more  native  point  defects  than   lower 

xi 


arsenic  pressure  conditions,  (c)  arsenic  antisite  (AS--O  defect  only 
can  be  observed  in  GaAs  grown  under  an  As-rich  or  high  arsenic 
pressure  conditions;  this  defect  can  not  be  produced  under  low  arsenic 
pressure  condi  ti  on. 

(2)  A  new  defect  model  supported  by  the  experimental  data  has 
been  developed  in  this  work  to  account  for  the  physical  origins  of  the 
EL2  electron  trap  in  GaAs.  It  is  shown  that  the  EL2  electron  trap  may 
be  attributed  to  two  different  types  of  native  defects:  One  is 
identified  as  the  EL2a  (i.e.,  Ec-0.83eV)  electron  trap,  and  the  other 
is  designated  as  the  EL2b  (Ec-0.76eV)  electron  trap.  The  physical 
origin  for  the  EL2a  level  is  attributed  to  the  doubly-charged  arsenic 
antisite  (i.e.,  AsGg++)  defect,  whereas  the  physical  origin  for  the 
EL2b  electron  trap  is  due  to  the  arsenic-anti site-plus-arsenic  vacancy 
complex  (ie.,  AsGaV,  ).  Based  on  this  model,  the  relationship  between 
the  density  of  EL2a  and  EL2b  trap  levels  and  the  [As]/[Ga]  ratio  in 
the  MOCVD  and  VPE  grown  GaAs  was  established. 

(3)  A  theoretical  model  for  the  nonexponential  DLTS  response  due 
to  the  field  dependent  emission  rate  of  the  trapped  charge  has  been 
developed.  A  comparison  of  the  theoretical  calculation  of  the 
nonexponential  DLTS  response  with  the  DLTS  data  for  each  trap  level 
would  allow  us  to  determine  the  potential  well  of  the  trap  involved. 
This  method  has  been  applied  to  identify  the  physical  origins  of  the 
EL2a  electron  trap  in  GaAs. 
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I .      INTRODUCTION 

Studies  of  deep-level  defects  in  GaAs  have  been  extensively 
published  in  the  literature,  t1-13^  Speculations  on  the  physical 
origins  of  defect  levels  and  defect  complexes  are  still  very 
tentative.  It  is  not  even  possible  to  say  with  any  assurance  that 
simple  gallium  vacancy  (V"Ga)  or  arsenic  vacancy  (V"As)  can  be  related 
to  a  particular  energy  level. 

In  one  respect  understanding  has  been  improved.  Gallium  arsenide 
(GaAs)  specimens  grown  by  various  techniques  such  as  liquid  phase 
encapsulation  (LEC) ,  Bridgmann,  vapor  phase  epitaxy  (VPE) , 
metalorganic  chemical  vapor  deposition  (MOCVD) ,  liquid  phase  epitaxy 
(LPE),  and  molecular  beam  epitaxy  (MBE)  are  now  recognized  as  likely 
to  have  different  properties  in  terms  of  energy  levels  within  the 
bandgap.  In  each  of  these  techniques  the  growth  temperature,  growth 
pressure,  growth  phase,  and  cooling  rate  are  usually  different.  For 
examples,  GaAs  grown  by  the  LPE  technique  from  a  gallium  melt  is 
expected  to  be  low  in  gallium  vacancy  defects  and  possibly  high  in 
arsenic  vacancy.  The  LPE  GaAs  usually  contains  hole  traps  with  energy 
level  of  Ey+0.71eV  (B  center),  while  VPE  GaAs  grown  in  the  As-rich 
conditions  always  contains  an  electron  trap  with  energy  of  Ec-0.83eV 
(EL2  level).  The  reason  for  observing  these  trap  levels  in  the  LPE  or 
VPE  GaAs   is   still  not  clear. 

The  main  objectives  of  this  research  are  (1)  to  make  a  detailed 
analysis  of  the  grown- in  defects  and  radiation- induced  defects  in  GaAs 
specimen   grown   by  the  LEC,    VPE,    LPE,    and   MOCVD   techniques   under 
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different  growth  and  annealing  conditions,  (2)  to  identify  the 
physical  origins  of  the  deep-level  traps  in  GaAs,  and  (3)  to 
determine  the  potential  well  of  electron  traps  by  analyzing  the  field 
enhanced  emission  rates  deduced  from  the  nonexponential  DLTS  response 
data. 

A  detailed  theoretical  and  experimental  study  of  the  grown-in 
defects  and  radiation  induced  defects  in  GaAs  has  been  carried  out  in 
this  dissertation;  these  include  the  following:  (1)  Theoretical 
modelling  of  native  defects  in  GaAs  has  been  developed.  GaAs  specimen 
grown  by  the  LEC,  VPE,  LPE,  and  MOCVD  techniques  under  different 
growth  and  annealing  conditions  have  been  measured  in  order  to  support 
this  model.  The  GaAs-samples  used  in  this  study  include  (i)  VPE  GaAs 
epilayers  grown  on  <100>,  <211A>,  and  <211B>  oriented  semi -insulating 
(S.  I.)  Cr-dopedGaAs  substrates  with  the  gas  phase  controlled  by 
varying  the  [Ga]/[As]  ratios  from  2/1,  to  6/1,  (ii)  LEC  n-GaAs 
samples  prepared  for  study  of  the  effect  of  hydrogen  heat  treatment  on 
deep-level  defects,  (iii)  LPE  GaAs  samples  grown  under  two  different 
temperatures  (i.e.,  700  and  800°C)  and  two  cooling  rates  (0.4  and 
l°C/min.),  and  (iv)  MOCVD  GaAs  epilayers  grown  on  S.  I.  GaAs  substrate 
and  S.  I.  Ge  substrates.  Theoretical  modelling  of  native  point  defects 
is  described  in  chapter  V.  Experimental  evidence  for  some  of  the 
native  point  defects  is  given  in  chapter  VII.  (2)  The  deep-donor 
trap,  commonly  known  as  EL2  center,  with  activation  energy  ranging 
from  E  -0.76  to  Ec-0.83eV,  has  been  observed  in  GaAs  grown  by  LEC, 
VPE,  and  MOCVD  techniques.  The  EL2  level  acts  as  a  recombination 
center  for  lifetime  reduction.  The  physical  origin  of  this  electron 
trap  has  been  a  subject  of  greater   interest  in  recent  years.     Although 


a  large  number  of  papers  have  been  devoted  to  finding  the  physical 
origins  of  the  EL2  electron  trap,  unfortunately,  none  of  these 
published  results  offered  a  consistent  and  unambiguous  explanation  for 
the  observed  EL2  trap  in  GaAs.  A  detailed  theoretical  and 
experimental  study  of  the  EL2  electron  trap  has  been  carried  out  in 
this  work.  The  detailed  results  are  discussed  in  chapter  VI.  (3) 
DLTS  signals  are  often  analyzed  assuming  that  the  capacitance 
transient  is  exponential.^''1  ^  However,  most  of  the  DLTS  signals  are 
nonexponenti al  transient.'-  '^  In  this  dissertation,  we  present  the 
theoretical  analysis  of  nonexponenti al  capacitance  transients  due  to 
the  electric  field  dependent  emission  rate  of  trapped  charge.  The 
emission  rate  of  trapped  charge  carriers  is  enhanced  by  the  Poole- 
Frenkel  and  phonon-assisted  tunneling  effects  in  the  presence  of  an 
applied  electric  field.  Since  the  electric  field  varies  with  position 
within  the  depletion  region  of  a  reverse  biased  p-n  junction,  the 
emission  rate  is  not  constant  within  the  same  region.  The  DLTS 
response  due  to  nonuniform  emission  rates  can  in  general  be  expressed 
as  S(  t)  =  E  exp(-enjt1)  -  E  exp(-en^t2).  Based  on  the 
nonexponenti al  capacitance  transient,  theoretical  calculations  of  DLTS 
response  for  deep-level  traps  in  GaAs  were  made  using  five  different 
potential  wells;  namely,  the  coulombic  well  which  has  a  positive 
charge  state  for  the  empty  state,  the  Dirac  well,  square  well, 
polarization  well,  and  dipole  well  which  all  have  neutral  charge 
states  from  different  physical  origins.  A  comparison  of  the 
theoretical  calculation  of  the  nonexponenti al  DLTS  response  with  the 
DLTS  data  for  each  trap  would  allow  us  to  determine  the  potential  well 
of  the  trap  involved.   This  model  was  applied  to  the  calculations  of 


the  potential  well  for  the  EL2a  level  and  other  electron  traps. 
Details  of  the  results  are  discussed  in  chapter  IV.  (4) 
Characterization  of  low  energy  proton  and  one-MeV  electron  irradiation 
induced  defects  in  the  LPE  grown  (AlGa)As-GaAs  solar  cells  has  been 
studied.  For  the  case  of  low  energy  proton  irradiation,  '■1'  ~  9J  the 
GaAs  solar  cells  were  bombarded  under  different  proton  energies  (50, 
100,  200,  and  290  keV),  proton  fluences  (1010,  1011,  5xl0U,  1012,  and 
1013  p/cm2)  and  Sn-dopant  densities  (2xl016  and  8xl016  cm-3).  For  the 
case  of  one-MeV  electron  irradiation,  l20~^^j  both  undoped  and  Sn- 
doped  GaAs  solar  cells  were  irradiated  under  different  temperatures 
(25  and  200°C) ,  fluences  (1014,  1015,  5xl015,  and  1016  e/cm2),  fluxes 
(2xl09  and  4xl010  e/cm2s) ,  and  annealing  times  (10,  20,  30,  and  60 
min.).  Defects  parameters  for  deep-level  traps  induced  by  low-energy 
proton  and  one-MeV  electron  irradiation  are  presented  in  chapter  VIII. 

Experimental  tools  employed  in  this  study  include  the  current- 
voltage  (I-V),  capacitance- voltage  (C-V) ,  thermally  stimulated 
capacitance  (TSCAP)  and  the  deep-level  transient  spectroscopy  (DLTS) 
measurements.  From  these  measurements,  one  can  determine  the  energy 
levels,  density  of  defects,  defect  profile,  and  the  capture  cross 
section  for  each  trap  level. 

Chapter  II  reviews  the  possible  native  defects  and  impurity 
complexes  in  GaAs.  Chapter  III  presents  the  experimental  details. 
Determination  of  the  potential  well  of  deep  level  traps  using  field 
enhanced  emission  rate  analysis  of  nonexponential  DLTS  in  GaAs  is 
described  in  chapter  IV.  Chapter  V  depicts  the  modelling  of  the 
grown-in  native  point  defects  in  GaAs.  The  physical  origin  of  the  EL2 
electron    trap    in   GaAs    is    explained    in    chapter   VI.      Chapter    VII 


describes  the  study  of  grown-in  defects  vs.  growth  parameters  in  VPE, 
LEC,  VPE,  and  MOCVD  grown  GaAs.  Defects  studies  in  low-energy  proton 
and  one-MeV  electron  irradiated  (AlGa) As-GaAs  solar  cells  are 
discussed  in  chapter  VIII.  Chapter  IX  gives  the  summary  and 
conclusions. 


II.      REVIEW  OF  POSSIBLE  NATIVE  DEFECTS   AND  COMPLEXES 
IN  GALLIUM  ARSENIDE 


The  number  of  possible  native  defects  in  GaAs  is  large  as  may  be 
seen  in  table  2.1.  None  of  these  native  defects  has  been  identified 
with  any  confidence.  This  is  due  to  the  fact  that  experiments  for 
studying  such  defects  tend  to  be  too  uncontrollable.  In  addition  to 
the  defects  shown  in  table  2.1,  other  types  of  defects  such  as 
impurity  defects,  impurity  complexes  defects  may  also  be  expected  in 
GaAs.  The  impurity  complexes  are  shown  in  Table  2.2  and  2.3.  In 
this  chapter,  we  will  focus  our  attention  on  a  few  native  defects 
listed  in  table  2.1,  which  are  believed  to  be  related  to  the  electron 
or  hole  traps  observed   in  this   study. 

Defects  can  be  represented  by  their  chemical  symbols.  For 
examples,  vacancy  is  represented  by  V  and  interstitial  by  i. 
Subscript  indicates  the  lattice  site.  Thus,  VGa  denotes  the  gallium 
vacancy;  Ga^  represents  the  gallium  interstitial  site  defect.  In 
addition  to  these  simple  point  defects,  defect  complexes  (for  example: 
VG  ~AsG  ++VG~)  might  be  expected  to  form  as  the  crystal  cooling  down 
from  high  temperature  to  room  temperature.  Another  type  of  defect, 
namely,  the  antisite  defect  must  also  be  considered.  The  antisite 
defects  (such  as  As^"1""1",  GaAs  ,  and  AsGa++GaAs  )  are  believed  to  be 
important  native  defects  in  GaAs  specimen. 

A  survey  of  the  literature  on  the  subject  of  defects  in  GaAs 
grown  by  the  Liquid  Encapsulation  Czockraski  (LEC) ,  Horizontal 
Bridgmann-Stockbargen  (HB) ,  Vapor  Phase  Epitaxy  (VPE) ,  Liquid  Phase 
Epitaxy   (LPE),   Molecular  Beam  Epitaxy   (MBE) ,   and  Metalorganic  Chemical 


Table  2.1     Possible  native  defects  in  GaAs. 

1.  One  Component  Defects: 

a.  Vacancy:  VGa,  VAs 

b.  Interstitial:  Gaj ,  As^ 

2.  Two  Component  Defects: 

a.  Di vacancy:  VGaVGa,  VGaVAs,  VAsVAs 

b.  Antisite:  Ga^s'  AsGa 

c.  Di-interstitial:  Ga^ASj 

d.  Vacancy-interstitial:  VGaGa.j ,  VGaASj  ,  VAgGaj ,  VAsAs^ 

3.  Three  Component  Defects: 

a.  Anti si te- vacancy  complexes:  AsGaVGa,  AsGgVAs,  GaAsVGa,  GaAgVAs 

b.  Tri vacancy:  V^V^V^,  VAsVGaVAs 

4.  Four  Component  Defects: 

a.  Anti  si  te-di  vacancy  complexes:   Vga^Ga^a'   VAs^aAsVAs 

b.  Di -anti  site:     GaAsAsGa 


Table  2.2  Sane  hypothetical   impurity  complexes. 

1 .  Impur  i  ty  i  nter sti ti  al :  A  j 

2.  Impurity  substitutional:  A^a,  AAg 

3.  impurity- vacancy  complexes:   A^V^,   A^V^,   AAsVGa,   AAsVAs 

4.  impurity-antisite  complexes:   A^As^,   AGaGaAs,   AAsAsGa,  AAsGaAs 

A  =  Group  II    impurities    (Be,  Mg,   Zn,   Cd) 
Group  IV  impurities    (C,   Si,  Ge,   Sn) 
Group  VI    impurities    (0,    S,   Se,   Te) 
Transition   impurities    (Ti ,   V,  Cr,    Mn,   Fe,  Co,  Ni ,  Cu) 


Table  2.3     Impurity  defects   in  GaAs. 


Periodic 


Table 


Group 

III 

Group 

IV 

Group  V 

Ga  atoms 

are 

/    ^ 

s 

As  atoms  are 

replaced 

by 

/ 

V 

replaced  by 

impurities 

impuri  ties 

*  Group  I;   Li,  Cu,   Ag,  Au 

This  group  of  impurities  will  act  as  acceptor  traps  in  GaAs. 
Each    impurity   may   contain    several    deep-level    traps. 

*  Group  II;   Be,   Mg,   Zn,  Cd 

This  group  of  impurities  will  act  as  shallow  acceptor  traps  in 
GaAs. 

*  Transition  Metals;   Sc,  Ti ,  V,  Cr,  Mn,  Fe,  Co,  Ni 

The  energy  level  for  each  of  these  transition  metals  has  a 
downward  trend  from  1.0eV  to  0.15eV  above  the  valence  band  as 
the  metal  elements  change  from  Ti    to  Ni. 

*  Group  IV;   C,   Si,  Ge,   Sn,   Pb 

When  each  of  these  impurities  is  on  a  Ga  site,  it  becomes  a 
shallow  donor  level.  On  the  other  hand,  if  the  atom  i  s  on  a 
As  site,    it  becomes  a  shallow  acceptor   level. 

*  Group  VI;    S,  Se,  Te 

This  group  of  impurities  will  act  as  shallow  donor  traps  in 
GaAs. 
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Vapor  Deposition  (MOCVD)  techniques  showed  that  only  a  few  electron 
and  hole  traps  are  common  point  defects  observed  in  both  the  bulk  and 
epitaxial  grown  GaAs  material,  while  most  of  the  point  defects  can't 
observed  in  GaAs  specimen  with  different  growth  processes.  Lists  of 
the  electron  and  hole  traps  observed  in  this  study  compared  with 
others  l-"-33J  are  gjven  respectively  in  tables  2.4  and  2.5. 

It  is  noted  that  EL2  (same  as  EF10,  ETl,  EB2,  ESI)  is  the 
dominant  electron  trap  in  the  VPE,  MOCVD,  and  bulk  grown  GaAs,  and  is 
absent  in  the  LPE  and  MBE  grown  GaAs  epitaxial  layers.  Trap  densities 
of  EL3,  EL5,  and  EL12  levels  are  varied  by  different  growth 
conditions,  and  density  of  these  traps  can  be  easily  reduced  by 
annealing.  This  suggests  that  they  are  due  to  native  point  defects. 
On  the  other  hand,  ELll  is  not  affected  by  heat  treatment;  therefore, 
the  physical  origin  of  ELll  level  could  be  due  to  an  impurity  complex. 
Crystals  grown  by  different  growth  techniques  are  expected  to  produce 
different  electron  traps,  and  only  a  few  of  them  are  common  and 
independent  of  methods  of  crystal  growth.  For  example,  both  EL6 
(bulk)  and  EL7  (MBE)  could  be  the  same  defect  in  GaAs,  which  is  not 
observed  in  VPE  grown  GaAs.  EL10  (MBE)  and  ELll  (VPE)  could  be  the 
same  level,  yet  not  seen  in  bulk  GaAs  material.  Most  of  these  deep 
level  defects  are  believed  to  be  either  due  to  the  vacancy  related 
defects,  antisite  defects  or  vacancy- i mpur i  ty  complex  defects.  For 
the  hole  trap  level,  HL2,  HL5,  HL7,  HL9,  and  HLll  may  be  related  to 
native  defects;  the  other  trap  levels  are  related  to  the  impurity 
defects. 
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Table  2.4.  Electron  traps  in  n-GaAs. 


Level 

Activation 

Capture  Cross 

Growth  method 

Ref. 

energy  (eV) 

Section  (cm2) 

ET1 

Ec-0.85 
0.30 

6.5xl0~13 

Bulk 

23 

ET2 

2.5xl0~15 

ii 

ii 

ESI 

0.83 

1.0X10"13 

Bulk 

24 

Ell 

0.43 

7.3xl0"16 

VPE 

25 

EI  2 

0.19 

l.ixnr*; 

VPE 

ii 

EI  3 

0.18 

2.2xl0~14 

VPE 

ii 

EB1 

0.86 

3.5x10-^ 
2.2xl0-ii 
3.0x10"^ 

Cr-doped  LPE 

26 

EB2 

0.83 

as-grown  VPE 

ii 

EB3 

0.90 

Electr .irradi . 

27 

EB4 

0.71 

8.3xl0~13 

ii     ii 

ii 

EB5 

0.48 

2.6xl0"13 

MBE 

28 

EB6 

0.41 

2.6xl0-13 

Electr.  irradi . 

27 

EB7 

0.30 

1.7xl0~j4 

MBE 

28 

EB8 

0.19 

1.5xl0"14 

MBE 

28 

EB9 

0.18 

- 

Electr.  irradi . 

27 

EB10 

0.12 

- 

n       H 

ti 

ELI 

0.78 

1.0Xl0~14 

1.2xl0"13 

Cr-doped  bulk 

29 

EL2  (A) 

0.83 

VPE 

30 

EL3  (B) 

0.58 

ii 

1.0xl0"^2 
1.2xl0"13 

M 

ii 

EL4 

0.51 

MBE 

EL5  (C) 

0.42 

VPE 

30 

EL6 

0.35 

1.5x10"^ 

Bulk 

29 

EL7 

0.30 

7.2xl0-1^ 
7.7x10"^ 

MBE 

EL8  (D) 

0.28 

VPE 

30 

EL9 

0.23 

6.8x10"^ 

ii 

ii 

EL10 

0.17 

1.8x10"^ 

MBE 

EL11  (F) 

0.17 

3.0xl0-1° 
4.9xl0-12 

VPE 

30 

EL12  (A) 

0.78 

ii 

ii 

EL14 

0.22 

5.2xl0-1° 

Bulk 

29 

EL15 

0.15 

5.7x10"^ 

Electr.  irradi . 

EL16 

0.37 

4.0X10"18 

VPE 

30 

EF1 

0.11 

Proton  i  rradi . 

18 

EF2 

0.14 

ii      n 

n 

EF3 

0.20 

ii      H 

ii 

EF4 

0.31 

Electr .irradi . 

21 

EF5 

0.35 

n-LEC 

EF6 

0.52 

Proton  irradi . 

18 

EF7 

0.60 

n-LEC 

EF8 

0.71 

Electr.  irradi . 

21 

EF9 

0.76 

LEC 

EF10 

0.83 

LEC,  VPE,  MOCVD 

31 

EF11 

0.90 

Electr.  irradi . 

21 
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Table  2.5     Hole  traps  in  p-GaAs. 


Level 

Activation 
energy  (eV) 

Capture  Cross 
Section  (air) 

Growth  method 

Ref. 

HT1 

Ey+0.44 

1.2xl0"14 

VPE 

23 

HS1 

0.58 

2.0xl0-1^ 

LPE 

24 

HS2 

0.64 

4.1xl0~l° 

ii 

ii 

HS3 

0.44 

4.8xl0"18 

n 

ii 

HBl 

0.78 

5. 2x10"  J- J 

Cr-doped  LPE 

26 

HB2  (B) 

0.71 

1.2xl0~j4 

as-grown  LPE 

ii 

HB3 

0.52 

3.4xl0_1° 
3.4xl0~14 

Fe-doped  LPE 

ii 

HB4 

0.44 

Cu-doped  LPE 

it 

HB5  (A) 

0.40 

2.2x10-^ 
2.0xl0~14 

as-grown  LPE 

H 

HB6 

0.29 

Electr .  irradi . 

27 

HL1 

0.94 

3.7xl0"]-4 

Cr-doped  VPE 

HL2 

0.73 

1.9xl0~14 

as-grown  LPE 

32 

HL3 

0.59 

3.0x10"^ 

Fe-doped  VPE 

33 

HL4 

0.42 

3.0xl0-1^ 
9.0xl0"|4 
5.6xl0~14 

Cu-doped  VPE 

ii 

HL5 

0.41 

as-grown  LPE 

32 

HL6 

0.32 

VPE 

33 

HL7 

0.35 

6.4xl0"15 

MBE 

11 

HL8 

0.52 

3.5xl0-1^ 

MBE 

11 

HL9 

0.69 

1.1x10"^ 

VPE 

11 

HL10 

0.83 

1.7xl0~13 

VPE 

11 

HF0 

0.082 

LEC  (p) 

17 

HF1 

0.13 

Electr.  irradi  . 

21 

HF2 

0.17 

Proton  irradi . 

18 

HF3 

0.20 

ii     ii 

n 

HF4 

0.29 

Electr.  irradi. 

21 

HF5 

0.35 

n      ii 

11 

HF6 

0.40 

ii      n 

11 

HF7 

0.44 

Proton  irradi . 

18 

HF8 

0.52 

ii 

HF9 

0.57 

ii 

HF10 

0.71 

Electr.  irradi. 

21 

III.   EXPERIMENTAL  METHODS 

Experimental  tools  employed  in  this  study  include  current- voltage 
(I-V)  measurement,  a.  c.  admittance  measurement,  capaci tance- voltage 
(C-V)  measurement,  thermally  stimulated  capacitance  (TSCAP) 
measurement,  and  deep  level  transient  spectroscopy  (DLTS)  measurement. 
From  these  measurements,  one  can  determine  the  diode  characteristics 
and  the  defect  parameters  such  as  defect  energy  levels,  defect 
densities,  and  capture  cross  sections.  Each  of  these  measurement 
techniques  is  described  as  follows: 

3.1.     Current-Voltage   (I-V)   Measurement 

Measurement  of  the  current-voltage  (I-V)  characteristics  under 
forward  bias  condition  can  yield  useful  information  concerning  the 
conduction  mechanisms,  recombination  processes  in  the  space  charge 
region  of  a  p-n  junction  diode  or  a  Schottky  barrier  diode.  For  a 
good  p-n  diode  with  no  surface  leakage,  the  total  current  is  composed 
of  the  diffusion  current  in  the  quasi -neutral  region  (QNR)  and  the 
recombination  current  in  the  space  charge  region  (SCR).  When  bulk 
diffusion  current  component  dominates,  the  current  expression  is  given 
by    W 

If  =  Id    [exp(Va/VT)    -  1]  (3.1) 

where  Va  is  the  applied  voltage;  VT  =  kT/q;  and  ld  is  the  magnitude  of 
a   saturation  diffusion  current. 
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Id  =  qn^  A    t(Dn/LnNa)    +    (DpApNd)]  (3.2) 

where  n^  is  the  intrinsic  carrier  density;  A  is  the  diode  area;  Dn 
(D_)  is  the  diffusion  constant  for  electrons  (holes);  L  (L  )  is  the 
diffusion  length  for  electrons  (holes);  and  Na  (Nd)  is  the  acceptor 
(donor)  density.  If  bulk  generation-recombination  current  dominates, 
then  the  current  expression  is  given  by   '■■'4J 

Ir  =  Irg  exp(Va/2VT)  (3.3) 

where  I  is  the  magnitude  of  a  saturation  generation-recombination 
current. 

Irg  =  qnjW  /  2xe  (3.4) 

In  Eq.(3.4),  W  denotes  the  depletion  layer  width;  x  =  (  t  *  )  ■*■/*  is 
the  effective  carrier  lifetime  in  the  space  charge  region;^ -*'  x  (tJ 
is  the  lifetime  of  electrons   (holes),    defined  by 

x  n  =  1  /   (Nta  n  <vfch»  (3.5) 

where  Nfc  is  the  trap  density;  a      (o    )    is  the  capture  cross  section 
for  electrons    (holes);    <vfch>   is  the  average  thermal  velocity. 
The  total  current  can  be  expressed  as 

It  -  If  +  IE  a      IQ  exp(V_/m*VT)  (3.6) 

where  I  is  the  saturation  current;  m  is  the  ideality  factor.  The 
diode  ideality  factor,  m  ,  can  be  used  to  identify  the  dominant 
current  component  in  a  p-n  diode.  Inspection  of  Eqs.(3.1)  and  (3.3) 
shows  that  the  bulk  diffusion  current  depends  more  strongly  on 
temperature  than  the  recombination  current  in  the  junction  SCR.     Since 
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the  recombination  current  in  the  junction  SCR  is  inversely 
proportional  to  the  effective  carrier  lifetimes  (and  hence  directly 
related  to  the  defect  density  in  the  transition  region),  measurements 
of  I-V  characteristics  under  forward  bias  conditions  would  allow  us  to 
determine  the  effective  lifetimes  in  GaAs  cells. 

3.2.     A.  C.  Admittance  Measurement 

A  useful  experimental  tool  for  evaluating  shunt  (R^  and  series 
resistance  (R  )  of  a  p-n  junction  diode  is  by  using  a.  c.  admittance 
measurement  techniques.  *■  '  A  p-n  junction  diode  can  be  represented 
by  a  three  element  device  with  shunt  resistance,  junction  capacitance, 
and  series  resistance,  as  shown  in  Fig.  3.1.  Measurements  of  the 
admittance  as  a  function  of  frequencies  (i.e.,  110  kHZ  to  700  MHZ) 
will  enable  us  to  determine  R  ,  B-,  and  C  of  the  diode.  The  impedance 
of  a  p-n  diode   is  given  by 

Z(w)    =  Rs  +   {Rpd/jwC)    /    [Rp  +    (i/jwC)]} 

=  R^w)    +  jXD(w)  (3.7) 

where 

R^w)    =  Rs  +    [Rp/d+w^Rp2)]  (3.8) 

XD(w)    =  RpWC  /   (l+w^Rp2)  (3.9) 

The  admittance  is  the  inverse  of  the  impedance  in  a  p-n  junction 
diode,  which  is  given  by 

Y(w)    =  1/Z(w)    =   {^(w)    /    [RD2(w)    +  XD2(w)]} 

-  j{XD(w)    /    [RD2(w)    +  XD2(w)]}  (3.10) 
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Figure  3.1  A  three  element  equivalent  circuit  model  for  a  p-n 
junction  diode. 
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In  the  low  frequency  and  high  frequency  limits, 

w  +  0,  Re  Y(w)  =  1  /  [Rp  +  Rs]  (3.11) 

w  -».»,  Re  Y(w)  =  1  /  Rs  (3.12) 

Thus,  from  a  plot  of  the  Im  Y(w)  vs.  Re  Y(w),  the  series  resistance 
and  shunt  resistance  of  the  diode  can  be  determined.  The  a.  c. 
admittance  technique  can  be  very  useful  in  evaluating  the  shunt 
leakage  problems  in  a  p-n  diode  and  allows  accurate  determination  of 
the  Rg,    R^,    and  C  over  a  wide  range  of  frequencies. 

3.3.     Capacitance-Voltage   (C-V)   Measurement 

The  capacitance-voltage  (C-V)  measurement  can  be  used  to 
determine  the  background  doping  concentration  in  the  n-  or  p-GaAs 
epitaxial  layers  using  a  Schottky  barrier  structure  or  a  one-sided 
abrupt  p  -n  (or  n  -p)  junction.  The  depletion  capacitance  across  the 
Schottky  barrier  diode  is  given  by 

C(Vr)    =e  gA  /  W  =  A   {qesNd  /    [2(4^    +  Vr   -  kT/q)  ]  }1/2         (3.13) 

where e s  is  the  dielectric  constant  of  GaAs;  ty:  is  the  built-in 
potential;  and  Vr  is  the  applied  reverse  voltage.  Equation  (3.13) 
shows  that  the  depletion  capacitance  of  a  Schottky  diode  is 
proportional  to  the  square  root  of  dopant  concentration  and  inversely 
proportional  to  the  square  root  of  the  applied  voltage.  If  the 
inverse  of  the  capacitance  square  (C  )  is  plotted  as  a  function  of 
the  reverse  bias  voltage  Vr,  then,  the  background  concentration,  can 
be  calculated  from  the  slope  of  C  vs  Vr  using  the  following 
expression 


C"2(Vr)    =    [2  /    (qe  SA2ND)  ]     ($.    +  Vr)  (3.14) 

The  intercept  of  C~2  vs.  Vr  plot  with  the  voltage  axis  yields  values 
of  (Jjj   which  is  related  to  the  barrier  height  of  a  Schottky  diode  by 

<pBn  =  $j   +  vn  +  kT/9  _A  <&  (3-15) 

where 

Vn  =  Ec  -    (kT/q)    ln(Nd/Nc)  (3.16) 

and  A  (J)  is  the  image  lowering  potential  of  a  Schottky  diode;     N     is 
the  density  of  state   in  the  conduction  band. 

3.4.     Thermally  Stimulated  Capacitance    (TSCAP)   Measurement 

Another  interesting  experiment,  which  is  known  as  the  thermally 
stimulated  capacitance  (TSCAP)  measurement  technique  i37'38J  will  be 
described  in  this  section.  The  TSCAP  experiment  is  carried  out  by 
first  reverse  biasing  a  p-n  diode  or  a  Schottky  diode,  and  then  the 
diode  is  cooled  down  to  the  liquid  nitrogen  temperature  (77°k).  After 
temperature  reaches  77°K,  the  diode  is  momentarily  zero  biased  to  fill 
the  majority  carrier  traps  and  returned  to  reverse  bias  condition,  and 
the  temperature  is  then  raised  from  77  to  400°  K.  The  thermal  scan  of 
capacitance  vs  temperature  plot  is  taken  by  an  X-Y  recorder.  A 
capacitance  step  is  observed  from  the  C  vs  T  plot  if  majority  or 
minority  carrier  emission  is  taking  place  in  a  trap  level  with  a 
certain  temperature  range.  The  amplitude  of  this  capacitance  step  is 
directly  proportional  to  the  trap  density.  The  trap  density  for  n- 
GaAs  can  be  calculated  from  the  following  expression 
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Nt  =  Nd(2AC  /  CQ)  (3.17) 

where  CQ  is  the  depletion  layer  capacitance  and  AC  is  the  capacitance 
change  due  to  the  majority  or  minority  carrier  emission.  Thus, 
knowing  N^  (or  Na)  and  CQ  at  the  temperature  where  the  capacitance 
step  was  observed,  the  trap  density  can  be  calculated  from  Eq.(3.17). 
Note  that  Eq.(3.17)  is  valid  only  for  the  case  when  Nt  <  0.1  Nj.  For 
the  case  of  large  trap  density  with  Nfc  >  0.1  Nj,  a  more  exact 
expression  should  be  used   instead. 

3.5.  DLTS  Measurement 

The  Deep-Level-Transient-Spectroscopy  (DLTS)  experiment  is  a  high 
frequency  (20  MHZ)  transient  capacitance  technique,  which  was 
introduced  first  by  Lang  in  1974.  '  ''14J  The  DLTS  scan  displays  the 
spectrum  of  deep  level  traps  in  the  forbidden  gap  of  a  semiconductor 
as  positive  or  negative  peaks  on  a  flat  baseline  as  a  function  of 
temperature.  Although  this  kind  of  measurement  is  time  consuming,  it 
offers  several  advantages  such  as  being  sensitive,  easy  to  analyze  and 
capable  of  measuring  the  traps  over  a  wide  range  of  depths  in  the 
forbidden  gap. 

Figure  3.2  shows  the  schematic  block  diagram  of  a  DLTS  system 
used  in  this  study.  The  DLTS  system  consists  of  a  sensitive 
capacitance  measurement  apparatus  with  good  transient  response,  a 
pulse  generator,  a  boxcar  with  dual-gated  signal  integrator,  an  X-Y 
recorder,    and  temperature  control   system   for  heating  and  cooling.     By 
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properly  changing  the  experimental  parameters  it  is  possible  to 
measure  the  following  parameters: 

-The  minority  and  majority  carrier  traps. 

-The  activation  energy  of  each  defect  level. 

-The  defect  concentration  which  is  directly  proportional  to 
the  peak  height. 

-The  defect  concentration  profile. 

-The  electron  and  hole  capture  cross  sections  for  each  trap. 

3.5.1.     Principles  of  the  DLTS  Technique 

The  capacitance  transient  is  associated  with  the  return  to 
thermal  equilibrium  of  the  carrier  occupancy  in  a  trap  level  following 
an  initial  nonequi  librium  condition.  The  polarity  of  the  DLTS  peak 
depends  on  the  capacitance  change  after  trapping  the  minority  or 
majority  carriers.  Because  an  increase  in  trapped  minority  carriers 
in  the  SCR  would  result  in  an  increase  in  the  junction  capacitance, 
the  minority  carriers  trapping  will  produce  a  positive  polarity  peak, 
and  vice  verse.  For  example,  in  a  p+n  junction  diode,  the  SCR  extends 
mainly  into  the  n-region,  and  the  local  charges  are  due  to  positively 
charged  ionized  donors.  If  a  forward  bias  is  applied,  the  minority 
carriers  (holes)  will  be  injected  into  this  region.  Once  the  holes 
are  trapped  in  a  defect  level  the  net  positive  charges  in  such  region 
will  increase.  This  results  in  a  narrower  SCR  width  which  implies  a 
positive  capacitance  change.  Thus,  the  DLTS  signal  will  have  a 
positive  peak.  Similarly,  if  the  majority  carriers  are  injected  into 
this  region  and  captured  by  the  majority  carrier  traps,  which  reduces 
the  local  charges,   the  SCR  width  will  be  wider,    implying  a  decrease  of 
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the  junction  capacitance.  Therefore,  the  majority  carriers  trapping 
will  result  in  a  negative  DLTS  peak.  The  same  argument  can  be  applied 
to  the  n+p  junction  diodes.  All  of  the  samples  used  in  the  DLTS 
measurements  are  p+n  diodes  so  that  a  positive  peak  represents  the 
hole  trap  and  a  negative  peak  represents  the  electron  trap. 

3.5.2.  Minority  Carrier  Injection 

Figure  3.3  shows  the  injection  of  minority  carriers.  Figure 
3.3(a)  is  the  injection  pulse  where  V=V  >  0  during  the  pulse  (ta  <  t 
<  tb)  and  V  =  VR  <  0  outside  the  time  interval  (tg,  tb) .  Figure 
3.3(b)  is  a  simplified  energy  band  diagram  (the  band  bending  due  to 
the  junction  electric  field  is  omitted)  in  quiescent  reverse  bias 
condition  (t  <  ta).  Figure  3.3(c)  is  the  saturating  injection  pulse 
(t  <  t  <tb).  During  the  pulse  the  minority  carrier  traps  Efc  are 
filled  by  holes.  Figure  3.3(d)  is  the  transient  process  during  the 
period  (t>tb);  the  capacitance  transient  is  due  to  the  trapped  holes 
begin  to  emit  from  the  trap  centers  to  the  valence  band  by  the  built- 
in  electric  field. 

3.5.3.  Majority  Carrier  Injection 

Figure  3.4  shows  the  injection  of  majority  carriers.  Figure 
3.4(a)  is  the  injection  pulse  where  V  =  V  =  0  during  the  pulse  (ta  < 
t  <  tb)  and  V  =  VR  <  0  outside  the  time  interval  (ta,  tb).  Figure 
3.4(b)  is  a  simplified  band  diagram  in  quiescent  reverse  bias 
condition  (t  <  t  ).  Figure  3.4(c)  is  the  majority  carrier  injection. 
During  the  pulse  period  (tg  <  t  <  tb) ,  the  majority  carriers  were 
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injected  into  the  SCR  and  captured  by  the  majority  trap.  Figure  3.4(d) 
is  the  transient  process  for  the  period  (t  >  th)  in  which  the  captured 
electrons  begin  to  emit  from  the  trap  centers  to  the  conduction  band 
and  then  they  are  swept  out  of  the  SCR  by  the  built-in  electric  field, 
resulting   in  the  observed  capacitance  transient. 

3.5.4.     Defect  Concentration 

The  defect  concentration  is  directly  proportional  to  the  peak 
height  as  described  before,  and  the  peak  height  is  proportional  to  the 
capacitance  change  a  C(0)  [  AC(0)  is  shown  in  Fig.  3.3(d)  and  3.4(d)], 
Therefore,  the  defect  concentration  Nt  is  proportional  to  A  C(0)  to  be 
derived  as  follows:  Let  C(t)  is  the  capacitance  transient^3  ,0-' 
which  proportional  to  electron  (hole)  emitted  to  the  conduction 
(valence)    band.  Then, 

C(t)    =  A{qes    [Nd  -  Nfcexp(-t/T)]    /    [2((Dj    +  Vr   +  kT/q)]}1/2 
=  CQ  {1  -    [Ntexp(-t/x)    /  Nd]}1/2  (3.18) 

where  t  is  time;  t  is  the  carrier  emission  time  constant;  CQ  =  C(Vr), 
as  shown  in  Eq.  (3.13),  is  the  junction  capacitance  at  the  quiescent 
reverse  bias  condition.  Using  binomial  expansion  and  the  condition 
that  Nfc/Nd  <<  1,  Eq.  (3.18)  reduces  to  a  simple  form  as 

C(t)    =  CQ    [1  -  Ntexp(-tA)/2Nd]  (3.19) 

Equation    (3.19)    can  be  rewritten  as 

Ntexp(-t/T)    =    (2AC(t)    /  CQ)    Nd  (3.20) 
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where  A  C(t)  =  C0  -  C(t).  From  the  DLTS  technique,  AC(0)  can  be 
measured.  The  junction  capacitance  CQ  and  the  background 
concentration  can  be  obtained  from  C-V  measurements.  Thus,  the  defect 
concentration  Nt  can  be  calculated  easily  by  using   Eq.(3.20)    at   t=0. 

3.5.5.     Activation  Energy  of  the  Defect  Level 

As  shown  in  Fig.  3.5, A  C(t)  begins  to  decay  after  the  injection 

pulse  is  over.     The  decay  associates  with  a  specific  time  constant 

which    is   equal    to   the   reciprocal    of    the   emission    rate.      For    an 

electron   trap,    the  emission   rate   "en"   is   a   functions   of   temperature, 

capture  coefficient  and  activation  energy,  and  can  be  expressed  by 
[41,42] 

^  =    (an<vth>Nc  /  9)    exp[(Ec  -  Eti)/kT]  (3.21) 

where  Efc;j  is  the  activation  energy  of  the  trap;  g  is  the  degeneracy 
factor;  a  is  the  electron  capture  cross  section  which  is  temperature 
dependent   *■     ' 

°  n  =  °oo  exP("  AEt/kT)  (3.22) 

where  om  is  the  capture  cross  section  at  very  high  temperature;  AEh  is 
the  activation  energy  of  the  capture  cross  section  for  the  trap.  The 
emission  rate  can  be  written  as 

en  =  BT2  exp{[Ec  -    (Efcj    +A^)]    /  kT} 

=  BT2  exp    [(Ec  -     Ejj,)    /  kT]  (3.23) 
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Figure  3.5  The  transient  capacitance  decay  exponentially  after 
the  injection  pulse  is  over. 
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where  B  is  the  proportionality  constant  and  is  independent  of 
temperature.  From  this  relation,  en  increases  with  increasing 
temperature.  The  capacitance  transient  is  rearranged  from  Eq.(3.19) 
as 

C(t)    =       CQ(Nt  /  2Nd)    exp(-tA) 

=    AC(0)    exp(-t  /T  )  (3.24) 

where  i  =  e  ~  .  Equation  (3.24)  implies  a  faster  decay  of  A  C(t),  as 
is  shown  in  Fig.  3.5,  at  higher  temperature  and  a  smaller  decay  time 
constant  which  implies  a  faster  decay  of  AC(t)  as  shown  in  Fig.  3.5. 
The  same  argument  can  also  be  applied  to  the  hole  trap.  The  only 
difference   is  that  AC(t)    has  an  opposite  polarity. 

The  following  is  the  procedures  to  deduce  the  activation  energy 
of  a  defect  level.  We  first  set  t-^  and  t2  in  a  dual-gated  integrator 
boxcar,   then  we  can  write 

C(tx)    =     C(0)    exp(-t1/x)  (3.25) 

C(t2)    =     C(0)    exp(-t2/x)  (3.26) 

The  DLTS  scan  along  the  temperature  axis  is  obtained  by  taking  the 
difference    of   Eqs.   (3.25)    and    (3.26)    which   yields 

S(  T)    =   A  C(0)[exp(-t1  /t  )    -     exp(-t2  /t    )]  (3.27) 

The  maximum  emission  rate,  x  max  ,  is  obtained  by  differentiating 
S(  x )  wi  th  respect  to  t  ,  and  let  i  t  equal  to  zero  [dS(r  )  /  dx  =  0] . 
The  result  yields 

Tmax  "    <fcl  "  fc2>   /  mfVty  (3'28) 
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Under  this  condition,  S(t  )  reaches  its  maximum  value  at  a  specific 
temperature.  The  emission  rate  is  given  by  en  =  1  /  Tmax  for  each  t-^ 
and  t2  setting.  By  changing  t-,  and  t2  several  times,  a  set  of 
temperatures  that  correspond  to  this  set  of  Tmax  (or  emission  rate  en) 
can  be  obtained  as  shown  in  Fig.  3.6.  The  activation  energy  of  the 
trap  can  be  calculated  from  the  slope  of  the  Arrhenius  plot  [i.e., 
ln(en/T2)  vs.  1/kT]. 
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IV.  DETERMINATION  OF  POTENTIAL  WELL  OF  DEEP-LEVEL  TRAPS  USING  FIELD 
ENHANCEMENT  EMISSION  RATE  ANALYSIS  OF  NONEXPONENTIAL  DLTS    IN  GAAS 


Deep-Level  Transient  Spectroscopy  (DLTS)  signals  are  often 
analyzed  assuming  that  the  capacitance  transient  is  exponential. 
However,  most  of  the  DLTS  signals  are  nonexponential  transient. 
[15,16]  Nonexponential  capacitance  transients  may  occur  from  (1)  the 
electric  field  dependent  emission  rate  of  the  trapped  charge,  L  4,  5J 
(2)  multiexponential  decay  due  to  several  trap  levels  with  similar 
emission  rates,  and  (3)  trap  density  of  the  same  magnitude  as  that  of 
the  shallow  dopant  density.  ^46'47^  This  chapter  will  deal  with  cases 
(1)    and    (2). 

To  determine  the  nonexponential  capacitance  transient  due  to 
electric  field  dependent  emission  rate  of  the  trapped  charge,  Makram- 
Ebeid,  t48'49!  Wang  and  Li  [50,51]  ^^q  measured  the  field  enhancement 
emission  rate.  In  this  chapter,  we  present  the  theoretical  analysis 
of  nonexponential  capacitance  transients  due  to  the  electric  field 
dependence  emission  rate  of  trapped  charge  and  multiexponential  decay. 
Emission  rate  of  trapped  charge  carriers  is  enhanced  by  the  Poole- 
Frenkel  effect  and  phonon-assisted  tunneling  effect  in  the  presence  of 
an  applied  electric  field.  Since  the  electric  field  varies  with 
position  within  the  depletion  region  of  a  reverse  biased  p-n  junction, 
the  emission  rate  is  not  constant  within  the  same  region.  The  DLTS 
response  due  to  the  nonuniform  emission  rates  can  in  general  be 
expressed  by  S(t)  =  E  exp(-en:jt1)  -  E  exp(-en;j  t2).  Based  on  the 
nonexponential  capacitance  transient,  theoretical  calculations  of  DLTS 
response    for    the    deep-level    traps    in    GaAs    are    made    using    five 
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different  potential  wells;  namely,  the  Coulombic  potential  well  which 
has  a  positive  charge  state  for  the  empty  state,  the  Dirac  well, 
square  well,  polarization  well  as  well  as  dipole  well  which  all  have  a 
neutral  charge  state  from  different  physical  origins.  A  comparison  of 
the  theoretical  calculations  of  nonexponential  DLTS  response  with  the 
DLTS  data  for  each  trap  level  would  allow  us  to  determine  the 
potential  well  for  the  trap  involved.  This  model  was  applied  to  the 
calculations  of  potential  well  for  the  EL2  level  and  other  electron 


traps    in    GaAs.       We    found    that    the    EL2    electron    trap   was   due    to 

5Ga 


Coulombic  potential  well  with  a  double  charge  state  (i.e.,   As~        ). 


Details  of   the  results  are  discussed  in  this  chapter. 

4.1     Capture  and  Emission  process  at  a  Deep-Level  Trap 

In  the  study  of  electric  field  dependence  of  emission  rate,  it 
is  important  to  know  the  charge  state  of  a  trap  so  that  the  type  of 
potential  well  for  such  a  trap  can  be  determined.  The  charge  state  of 
a  deep-level  trap  may  be  either  positively  charged,  neutral,  or 
negatively  charged  ^10^  (i.e.,  Nt+,  Nt°,  or  Nfc~).  For  example,  if  an 
electron  is  captured  by  a  positive  donor  trap,  then  the  kinetic 
equation  for  the  capture  process  requires  that 

Nt+  +  e  =   Nt°  (4.1) 

In  this  case,  the  capture  process  is  Coulomb  attractive.   Similarly, 
the  emission  process  for  the  same  trap  can  be  written  as 

Nt°  =  Nt+  +  e  (4.2) 
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which  again  is  an  attractive  process.  This  attractive  capture- 
emission  process  is  the  feature  of  a  Coulombic  potential  well.  In 
GaAs,  for  example,  the  deep-level  defects  such  as  arsenic  vacancy 
^As+)  aTV^  arsen:>c  antisite  L13J  (AsGa++)  defects  have  a  Coulombic 
potential  well. 

For  an  empty  neutral  trap  with  Nt°,  the  capture  process  is 
given  by 

Nt°  +  e  =  Nt"  (4.3) 

and  the  emission  process  for  the  same  trap  is  given  by 


Nt~  =  Nt°  +  e  (4.4) 


The  potential  wells  for  a  deep-level  neutral  trap  in  a  semiconductor 
may  include  Dirac  delta  potential  well,  *■  *  square  potential  well, 
polarization  potential  well  and  dipole  potential  well.  *■  '  The  Dirac 
delta  potential  well  does  not  show  Poole-Frenkel  effect.  Neutral 
arsenic  vacancy  (VAs°)  defect  in  a  GaAs  crystal  behaves  as  a  Dirac 
delta  potential  well  in  the  capture-emission  process.  For  the  square 
potential  well,  a  potential  well  has  a  depth  of  VQ  for  r  <  rQ  and 
zero  for  r  >  rQ.  Deep-level  defects  such  as  neutral  arsenic 
interstitial  (ASj°)  have  these  characteristics.  The  potential  of  a 
polarization  potential  well  has  the  form  V(r)  =  -  A^  /  r  ,  where 
A-l  =  e  Oj  /  8e0es  *;  a.  is  the  polar  i  zabi  li  ty  of  a  neutral  impurity 
atom,  and  e  is  the  relative  dielectric  constant  of  the  host  crystal. 
Lax  I5  J  used  a  polarization  potential  well  to  model  capture  of 
electrons  and  holes  by  neutral  impurities.  Tasch  and  Sah  l  '  used 
this  potential  well  to  model  the  observed  field  dependence  of  Au  in 
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silicon.  A  neutral  impurity  trap  such  as  AlGa°  has  the  properties  of 
a  polarization  potential  well.  A  complex  of  two  oppositely  charged 
centers,  each  with  charge  magnitude  of  Zq,  will  form  a  dipole 
potential  well.  The  potential  is  represented  by  a  dipole  potential, 
V(r)  =  qr-P  /  4eQesr  ,  with  a  dipole  moment  given  by  P  =  Z  q  R  , 
where  RQ  is  the  distance  separating  the  two  charge  centers.  The 
antisite  pair  defect  (AsGa++GaAs  )  behaves  as  a  dipole  potential 
well   in  the  capture-emission  process. 

For  the  process  Nt~  +  e  =  Nt  ,  the  trapping  process  is  repulsive 
and  the  trap  in  Nt~  state  has  a  Coulombic  energy  barrier  that  the 
electron  must  surmount  thermally  or  tunnel  through.  Thus  one  may 
expect  its  capture  cross  section  to  be  very  small    (^  10  cm^  or 

less)    and  with  a  strong  temperature  dependent  property. 

In  DLTS  technique,  l'»"i  electrons  captured  by  trap  level 
occurred  in  the  quasi -neutral  region  (QNR)  of  a  p-n  diode  which  has 
zero  electric  field,  as  shown  in  Fig.  4.1(a).  On  the  other  hand, 
electrons  emission  from  a  trap  level  usually  take  place  in  the  space- 
charge  region  (SCR)  of  a  p-n  diode  in  which  high  electric  field 
prevails  in  this  region,  as  shown  in  Fig.  4.1(b).  The  effect  of 
electric  field  on  the  emission  process  is  to  enhance  the  emission  rate 
by  either  Poole-Frenkel  effect  or  phonon-assi sted  tunneling  effect  to 
be  discussed  in  next  section.  The  emission  rate  enhanced  by  the 
electric  field  will  cause  the  DLTS  response  to  be  nonexponential 
transient. 


a:  Capture    Process 


b:   Emission 


Figure  4.1  (a)  Electrons  captured  by  trap  level  are  occurred  in  the 
quasi-neutral  region  of  a  p-n  junction. 

(b)  Electrons  emission  from  a  trap  level  usually  take 
place  in  the  space  charge  region  of  a  p-n  junction. 


36 


4.2     Overview  of  the  Theory 

The  emission  rate  is  based  on  the  detailed  balance  expression, 
141, 42J  an(j  uncjer  zero  electric  field  condition,  Eq.  (3.23)  is 
rearranged  as 

eno  =  en°°  exP(-  Eti/kT)  (4,5) 

where  Efcj  is  the  ionization  energy  of  the  trap  level,  and  e^  contains 
the  matrix  elements  of  the  transition.  Figure  4.2  shows  a  trap  level 
under  high  electric  field.  There  are  three  basic  mechanisms  I-54-! 
which  affect  the  emission  rate  enhancement  under  high  electric  field 
conditions. 

(1)  Poole-Frenkel    effect,     in    which    electrons    climb    over    a 

barrier  lowered  by  the  presence  of  applied  electric  field. 

(2)  Phonon-assi sted  tunneling  effect,  in  which  electrons  absorb 

thermal  energy  from  the  lattice,    and  then  tunnel   through  the 
barrier  at  a  higher  energy. 

(3)  Pure  quantum  mechanical  tunneling  effect.  Since  pure 
tunneling  becomes  important  only  at  very  high  electric  fields  L  J 
(i.e.  F  >  10  V  /  cm),  we  will  not  consider  this  effect  in  our  DLTS 
analysis. 

4.2.1  Poole-Frenkel  effect 

The  Poole-Frenkel  and  Schottky  lowering  effects  are  due  to  the 
lowering  of  a  Coulombic  potential  barrier  by  an  applied  electric 
field.  The  Schottky  effect^4J  is  associated  with  the  lowering  of 
barrier   height   of     a  metal-semiconductor   contact.      The   Poole-Frenkel 
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Figure  4.2  There  are  three  basic  mechanisms  which  affect  the  emission 
enhanced  under  high  field  conditions. 


38 


effect  is  associated  with  the  lowering  of  potential  barrier  in  a  deep- 
level  trap  of  a  bulk  semiconductor.  Both  donor  and  acceptor  traps 
would  show  the  Poole-Frenkel  effect.  The  enhancement  of  emission  rate 
from  a  Coulombic  potential  well  due  to  the  Poole-Frenkel  effect  was 
first  done   in  a  one-dimensional  model  by  Frenkel,    l*'i    and  later 

extended  to  three-dimensional  case  independently  by  Harke   '      ■*   and 

r  59 1 

Jonscher.  l   J  A  samalar  calculation  for  the  polarization  well  and 


dipole  well  was  done  by  Martin  et  al.. 


[54] 


In  this  section  we  will  consider  the  general  case  for  calculating 
the  emission  rate  enhancement  due  to  Poole-Frenkel  effect.  The 
potential  for  a  deep-level  trap  in  an  electric  field,  F,  can  be 
expressed  by 

VT(r)    =   V(r)    -  q   F   r   cos (6)  (4.6) 

where  V(r)  is  the  potential  of  a  deep-level  trap.  For  0  <  6  <  tt  /2,  the 
minimum  potential  was  found  by  setting  dVT(r)/dr  =  0  at  r  =  rm  . 
The  lowering  of  potential  barrier  due  to  the  presence  of  an  applied 
electric   field   is   found  by  evaluating  Vm(r)   at  r  =  rmax.   Thus, 

AEti(6)   =VT(rmax)  (4.7) 

A  one  dimensional  result  is  given  by  setting  6=0.  The  Poole-Frenkel 
effect  leads  to  a  decrease  of  the  ionization  energy,  and  the  one 
dimensional  thermal  emission  rate,  eni,  in  the  presence  of  an  electric 
field  can  be  written  as 

enl  =  en°°  exp[-(Eti   -AEtj)   /  kT] 

=  enQexp(    A  EtiAT)  (4.8) 
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where  eno  is  the  thermal  emission  rate  at  zero  electric  field  given  by 
Eq.  (4.5).  The  three-dimensional  calculation  requires  an  integration 
over  e  due  to  the  spatial  variation  of  AEtj(9).  For  0  <  9  <tt/2,  the 
emission  rate  is  proportional  to  8,  and  for  tt  /2  <6<  tt,  it  is 
assumed  that  the  emission  rate  does  not  change  with  electric  field. 
Thus, 

2  TT  /2  2  it  TT 

en3/eno=(l/4Tr)[  /  d'pJ       sin(e)exp( A  Ej/kT)d     +1  d«j,    I      sin(6)do] 
0       0  '  0       t/2 

(4.9) 

where  eno  is  the  three  dimensional  thermal  emission  rate. 

4.2.2       Phonon-assisted  tunneling  effect 

There  are  several  different  types  of  potential  wells  which  may 
exist  in  a  deep-level  trap,  depending  on  the  charge  state  of  a 
particular  trap.  The  potential  well  may  be  affected  by  an  external 
electric  field.  Figure  4.3  shows  a  potential  well  for  a  trap  level 
with  energy  of  Efc:j  as  a  function  of  the  radial  coordinate.  A  trapped 
electron  can  absorb  a  phonon  and  tunnel  through  the  barrier.  The 
tunneling  process  can  be  treated  by  using  the  WKB  [®\b]  (wentzel- 
Kramers-Bri llouin)  approximation  for  the  potential  barrier.  The 
probability  ratio,  ^61^  (or  the  barrier  penetration  factor,  P)  for  an 
electron  impinging  on  the  barrier   is  given  by 

<2 

y 

(i 


X 

P  =  exp[-2  I    ~k(x)    dx] 
x 

x2 


=  exp{-[(8m)1/2  /  Jflj       [V(x)-E]1/2dx}  (4.10) 

xl 


4C 


Figure  4.3  A  trapped  electron  can  absorb  a  nhonon  and  tunnel  through 
the  barrier. 
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where  k(x)  =  [(2m)1/2/H]  x  [V(x)-E]1/2  is  positive  when  V(x)  >  E.  The 
potential  well  for  the  trap  level  is  V(x).  The  electron  energy  is  E. 
The  two  turning  points,  x-^^  and  x2,  are  the  fixed  points  for  which 
electron  energy  is  equal  to  the  barrier  energy  of  the  potential  well 
for  the  trap. 

A  trapped  electron  can  absorb  a  phonon  with  energy  of  exp(- 
Eth/kT)  and  tunnel  through  the  barrier  at  a  higher  energy.  The 
probability  of  this  composite  event  is 

Pc  =  exp(-Eth/kT)   P  (4.11) 

where  Et^  is  the  thermal  excitation  energy  or  phonon  energy;  and  P  is 
given  by  Eq.  (4.10).  The  emission  rate  due  to  phonon-assi  sted 
tunneling  effect  is  given  by  integrating  Eq.  (4.11)  over  the  phonon 
energy,    Efch. 

Eti/kT 

ent=J  pc  dEth  /  kT  (4.12) 

0 

where  1/kT  is  a  normalization  factor.  The  normalized  emission  rate 

enhancement  is  obtained  by  dividing  Eq.(4.12)     by  the  zero  electric 

field  emission  rate  e^,   which  yields 

EtjAT 
ent  /  eno  =  exp(EtJAT)'  Pc  dEfch  /  kT 

0 

Eti /kT 
=  exp[(Etj-Eth)AT]  ;  exp{[-(8m)1/2/K] 

0 
x2 
J       [V(x)-E]dx     dE^  /  kT}  (4.13) 


X^ 


*1 

Equation     (4.13)     is    a    general    expression    for    the    emission    rate 
enhancement  due  to  phonon-assi sted   tunneling  effect. 
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The  total  emission  rate  enhancement  is  obtained  from  the  sum  of 
Poole-Frenkel  effect  and  phonon-assi sted  tunneling  effect.  From 
Eqs.  (4.9)  and  (4.13),  one  obtains 

^H  =  eno  <  ^3  +  ent  >  <4-14) 

Vincent  et  al.  *■  *  and  Martin  et  al.  *■  ••  showed  that  for  a  Coulombic 
well  and  other  potential  wells,  both  the  Poole-Frenkel  effect  and 
phonon-assi sted  tunneling  effect,  are  important  over  the  field  range 
of  interest  (104  -  106  V/cm  at  300°  K).  Our  results  are  in  good 
agreement  with  this  observation,   as  will  be  shown  next. 

4.3     Theoretical  calculations  of  the  emission  rate  enhancement 
for  different  potential  wells 

4.3.1     Coulombic  potential  well 

A  donor  type  electron  trap,  which  has  a  Coulombic  potential  well, 
is  positively  charged,  Nt+  when  it  is  empty.  The  capture  (emission) 
process  for  a  such  trap  is  given  by  Eqs.  (4.1)  and  (4.2).  Assuming  an 
electron  is  at  a  distance,  r,  from  the  trap,  the  attractive  force 
between  a  positive  trap  and  the  electron  is  given  by 

F(r)   =   -q2   /   4  e^r2  (4.15) 

The  potential    for    the   trap   is   obtained  by   integrating   Eq.    (4.15)    from 

infinity  to  r.   Thus, 

r 
V(r)    =  I    F(r)    dr      =   -q2   /  4  it  e0esr  (4.16) 
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When  an  external  field,  F,  is  applied,  the  total  potential  energy  as  a 
function  of  distance  is  given  by 

VT(r)  =  -q2/4  7re0Esr  -  qFrcos(e)  (4.17) 

Electron  traps  such  as  As  vacancy  (V"As+)  and  As  antisite  (AsGa++) 
defects  have  Coulomhic  potential  well  for  electrons.  Equation  (4.17) 
can  also  be  applied  to  an  acceptor  type  hole  trap,  which  also  has  a 
Coulombic  potential  well.  In  this  case,  the  Coulombic  potential  well 
is  negatively  charged  when  it  is  empty.  The  capture  process  for  such  a 
hole  trap  is 

Nt~  +  h+  =  Nt°  (4.18) 

and  the  emission  process  is 

N^°   =  Ni."   +   h+  (4.19) 


Hole  traps  such  as  Ga  vacancy  (VGa~)  and  Ga  antisite  (Ga«  )  defects 
have  a  Coulombic  potential  well  for  holes.  Compared  to  other  types  of 
potential  wells,  a  Coulombic  potential  well  shows  a  very  large  Poole- 
Frenkel  effect. 

4.3.1.1  Poole-Frenkel  effect 

The  Poole-Frenkel   lowering  potential     Etj   and  the  location  of 
lowering  rmax  are  given  by  the  condition  dVm(r)/dr  =  0,   or 

rmax  =     (q  /   [4  tt  e  e  F  cos(6)]}1/2 

(4.20) 


AEfci   =  q   [    q  F  cos(6)   /   (  vzQzs)]1/2 
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The  Poole-Frenkel  effect   leads  to  a  decrease  of  ionization  energy 
Efcj,   and  the  one-dimensional  thermal  emission  rate  can  be  written  as 

enl  =  enoexp(AEtiAT)  (4.21) 

For  a  three-dimensional  case,    the  emission  rate  of  a  trapped  electron 

may  be  obtained  by  assuming  that  emission  rate  is  field  dependent  for 

0  <    Q  <  tt/2  and  is  independent  of  electric  field   for     t/2  <  9    <  tt  . 

Thus, 

2  it     tt/2 
en3    =    (enc/47r  )     [J        J     sin(e)    exp(AEti/kT)de  d<j> 
0       0 
2tt     tt/2 
+  J      j  sin(e)de   d  4>    ]  (4.22) 

0       0 

Integrating    Eq.   (4.22)  ,  and   the   result   is 

en3/eno  =   (1/R2)  [eR(r-l)  +  1]    +  1/2  (4.23) 

where 

R  =    q    (qF/TrE0£s)1/2    /  kT  =    A  Efci    /kT  (4.24) 

4.3.1.2       Phonon-assi sted  tunneling  effect 

The    phonon-assi sted    tunneling    effect    may    be    evaluated    by 

substituting    Eq.    (4.17)    to    Eq.    (4.13).      Vincent    et    al.     ^53^     have 

obtained  a  closed   form  for   the  phonon-assi sted  tunneling  effect  for  a 

Coulombic  potential  well,  which  reads 

(Efci-  AEti)/kT 
enfc/eno    =j  exp{z   -   7.i/l  [4  (2m) l/z  (kT)  i/z 

0 

/   3qtfF]     [1   -    (AEti/zkT)5/3]}   dz  (4.25) 


45 


The  total  emission  rate  is  composed  of  Poole-Frenkel  effect  and 
phonon-assi sted  tunneling  effect 

enHC   =   eno   (^3    +   ^  (4'26> 

Figure  4.4  shows  the  normalized  enhanced  emission  rate  vs.  electric 
field  for  the  Coulombic  potential  well  for  the  Ec  -  0.83  eV  electron 
trap  oberved  in  GaAs.  The  Poole-Frenkel  effect  dominates  when 
electric  field  is  lower  than  2  X  10  V/cm,  while  the  phonon-assi sted 
tunneling  effect  becomes  important  for  electric  fields  higher  than  2  x 
104  V/cm. 

4.3.2     Dirac  delta  potential  well 

The  potential  for  a  Dirac  delta  well   is  given  by 


V(r)   =  -V  for  r  =  0 

(4.27) 
=0  for  r  *  0 


In  the  presence  of  an  applied  electric  field,  the  potential  is  given 
by 

V(r)   =  -V     -  qFr  cos (6)      for    r  =  0 

(4.28) 
=  0  -  qFr  cos (6)  for   r  p  0 

A  neutral  vacancy  trap  such  as  neutral  arsenic  vacancy  (V"As°)  may  have 
a  Dirac  potential  well.  For  this  potential  well,  the  Poole-Frenkel 
lowering  potential  is  equal  to  zero,  and  only  the  phonon-assi sted 
tunneling  effect  is  considered.     For  the  case  of  a  one-dimensional 
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Figure  4.4  Normalized  emission  rate  enhancement  vs.  electric  field 

for  the  Coulombic  potential  well  for  EL2a  electron  trap 

for  phonon-assisted  tunneling  effect  and  Poole-Frenkel 
effect. 
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calculation,  cos(  )  is  equal  to  1,  and  the  emission  rate  enhancement 

is  obtained  by  substituting  Eq.    (4.28)    to   Eq.    (4.13),    which  yields 

E    -/kT 
ent/enQ  =       {exp[(EtJ-Eth)/kT]};       *     {exp  [-(8m)  1/2/Kl 

0 

[j      (-qFx  +     Etj  -   Efch)1/2   dx]       dEfch/kT}         (4.29) 
xl 
Equation    (4.29)    gives   an  analytical    form   for   the   transparency  of 

triangular  barrier  of  a  Dirac  well.     For  the  triangular  well,  x^  =  0 

and    X2    =    (Efcj    -   Efch)    /   qF.       Integrating    Eq.     (4.29)    over    Efch    and 

letting  z  =   (Efc;j   -  Eth)  /  kT,    we  obtain  the  enhanced  emission  rate  for 

the  phonon-assisted  tunneling  effect 

E   -/kT 
ent  /  eno  =  '   t3        exP{z  "  z3/2t4(2m)1/2(kT)3/2/3qMF]}  dz 
0 

(4.30) 

The  total  emission  rate  due  to  Dirac  well  can  be  expressed  by 


enHR  "  eno 


(  1  +  e^)  (4.31) 


4.3.3     Square  Potential  Well 

The  same  treatment  can  be  applied  to  the  square  well.  The 
potential  of  a  square  well  in  the  presence  of  an  applied  electric 
field  is 

V(r)  =  V    -  qFr  cos  (e)  for   r  <  r 

(4.32) 
=  0  -  qFr  cos(9)  for  r   >   rQ 

where  rQ   is  the  radius  of  the  trap.      A  neutral   interstitial  trap  such 
as  neutral  arsenic  interstitial   (ASj°)   may  have  such  a  potential   well. 
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The  barrier  lowering  due  to  Poole-Frenkel  effect  is  equal  to  qFr  . 
The  enhanced  emission  rate  for  a  one-dimensional  and  a  three- 
dimensional  Poole-Frenkel  effect  is  given  respectively  by 

enl/eno  =   exP(R> 

(4.33) 

en3   /  eno  =    (X/2R)    x    texp(R)    -  1]    +    1/2 

where  R  =  qFr^kT. 

For  the  enhanced  emission  rate  due  to  phonon-assi sted  tunneling 
effect,  Eq.  (4.29)  is  still  valid.  In  a  square  well,  x-^  =  rQ  and  x2  = 
(Etj  -  Et^)/kT.  The  range  of  phonon  energy,  Eth  is  from  0  to  (E.j  - 
AEtj)  /  kT.  Thus,  we  obtain  an  expression  for  the  normalized  enhanced 
emission  rate  due  to  phonon-assi sted  tunneling  effect: 


(Etj-  AEti)AT 
snt/  eno   =!  exp{z-z-i/^[4(2m)i/^(kT):!/^ 


■/3qKF] 


[1  -    (qFr^kT)3/2]}     dz  (4.34) 


The   total   emission  rate   is    the   sum   of   Eqs.   (4.33)    and    (4.34) 

enHS    =   ^o    (^3   +   ^W  <4-35) 


4.3.4     Polarization  potential  well 

An  electron  with  charge,  e,  at  a  distance  r  from  a  trap  center 
with  polar i  zability  ,  a  j,  will  induce  a  dipole  moment,  p  =  a-e/e  r  . 
This  dipole  will  produce  an  attractive  force  on  the  charge  of 
pe  /  (2  tt  e  Q  e  sr3)  =  cx^e2  /  (2  it  eQes2r5).  Thus,  the  attactive 
potential   is  given  by 

V(r)    =    -A1  /   r4  (4.36) 
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where 

Al   =     «je2  /  (8  tie  Qes2) 

The  polarizability  of  the  atom  can  be  expressed  by 

Oj   /aH  =    (mQ/m)    (  EH/Efcj)2  (4.37) 

where  Etj  is  the  ionization  energy  of  the  trap;  EH  =  13.6  eV  is  the 
ionization  energy  of  a  hydrogen  atom;  a  „  =  0.666  X  10  cm3  is  the 
polarizability  of  a  hydrogen  atom;  a ^  is  the  polarizability  of  the 
trap;  mQ  is  the  free  electron  mass;  and  m  is  the  effective  electron 
mass.  Lax    >-^->*    used   the  polarization  potential    well    to   model 

capture  of  electrons  and  holes  by  neutral  impurities.  Tasch  and  Sah 
*■  *  used  this  potential  well  to  model  the  observed  field  dependence 
emission  rate  of  Au  in  silicon.  A  neutral  impurity  trap  such  as  AlQa° 
has  the  properties  of  a  polarization  potential  well.  In  the  presence 
of  an  electric  field,  the  total  potential  of  the  polarization  well  is 
given  by 

VT(r)    =  -Ax   /  r4  -   qFr  cos(e)  (4.38) 


4.3.4.1  Poole-Frenkel  effect 

The  minimum  potential  is  obtained  by  setting  dV"T(r)/dr  =  0  at  r  = 
rmax'  an<^  t'"ie  resulfc  is  given  by 

rmax  =    (4A1  /  qFcos(0)]1/5 

1 .-  ,/t.  (4.39) 

AEtj    =   -1.649   A^^    [qFcos(e)]    7 

Emission  rate  is  field  dependent  for  0  <  9   <  ^/2  and  is  independent  of 
electric  field  for  tt/2  <  0  <  it  .  Thus,  the  normalized  field  enhanced 
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emission  rate  due  to  Poole-Frenkel  effect  can  be  written  as 

enl/eno  =  exp(R)  (4.40.a) 

for  the  one  dimensional  case,  and 

en3/eno  =  (1/2)  x  [1  +  (5/4)  x  j  Rt1/4exp(R)  dt]   (4.40.b) 

0 

for  the  three  dimensional  case.      Where         R  =   1.649  A1//5    (qF)4'5  /  kT 

4.3.4.2  Phonon-assi sted  tunneling  effect 

Substituting  Eq.  (4.38)  into  Eq.  (4.13)  and  letting  z  =  (Etj  - 
EtvJ/kT  yields  the  normalized  enhanced  emission  rate  due  to  phonon- 
assi  sted  tunneling  effect: 

(Eti-AEti)AT 

ent/eno    =    j  exp{z    -    [(Sm)1/2/*]       (-A^r4 

0 

-  qFx    +    zkT)1/2    dx}    dz  (4.41) 

Numerical  integration  can  be  carried  out  by  assuming  x-^  =  (A-^  / 
zkT)1'      and  x^  =  zkT/qF.     The  total  enhanced  emission  rate  is 

^HP   =   eno  <en3   +   ent>  (4'42) 

Note  that  integrals  given  in  Eqs.  (4.40)  and  (4.41)  can  not  be  solved 
analytically.  The  numerical  calculations  of  the  enhanced  emission 
rates  as  a  function  of  the  electric  field  were  shown  in  Fig.  4.5  for 
Ec  -  0.83  eV  trap  in  GaAs  at  300°  K.  The  Poole-Frenkel  effect  is 
dominant  when  electric  field  is  lower  than  4  x  10  V/cm,  while  the 
phonon-assi sted  tunneling  effect  becomes  important  for  electric  field 
greater   than     4  x  10     V/cm. 
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Figure  4.5  Normalized  emission  rate  enhancement  vs.  electric  field 

for  EL2a  electron  trap  in  GaAs  for  phonon-assisted  tunneling 
and  Poole-Frenkel  effects. 
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4.3.5  Dipole  potential  well 

The  dipole  potential  can  be  represented  by  a  trap  center  with  two 
oppositely  charged  ions,  each  with  charge  Zq. 

V(r)  =  r-p  /  4  it  e0esr2  (4.43) 

where  p  is  the  dipole  moment,  p  =  ZqRQ;  RQ  is  the  distance  between  two 
opposi  tely  charges.  Antisite  pair  defect  (AsGa++GaAs  )  behaves  as 
the  dipole  potential  well  in  the  capture-emission  process.  As  is 
shown  in  Fig.  4.6,  if  the  applied  electric  field,  F,  forms  an  angle 
8  p  with  z  axis,  and  is  polarized  along  the  z  direction,  then  the 
total  potential  becomes 

VT(r)     =    -qPcos(6)/(4iT  e0esr2)    -   qFr   si  n  (8p)  si  n  (0  )  cos  (<(>) 
-  qrF  cos(eF)cos(6)  (4.44) 

4.3.5.1     Poole-Frenkel  effects 

The  mininum  potential   is  obtained  by  setting  dVT(r)/dr  =  0,  which 
yields 

rmax  =  *[p  cos(9)   /  (2lT£oes)]     [F  sinfe^sinP)  cos(<|>) 

+  Fcos(8EJcos(8)    ]}1//3  (4.45) 

and  the  change  in  barrier  height  is  given  by 

Eti(9,(fr)    =  -(3   x  2_2/3)    x    [P  cos(e)    /    (4tt  e0es)]1/3 

x    [F  sin(eF)sin(e)cos(<j))   +  F  cos(ep)  cos  (6)  ]  2/3 

(4.46) 
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Figure  4.6  In  dipole  potential  well,  external  electric  field  form 
an  angle  6..  with  the  dipole  moment. 
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Unlike  other  potentials,  the  dipole  potential  experiences  a  barrier 
lowering  which  is  field  dependent.  The  integrals  fall  into  two  cases 
according   to  the   incident  angle    Q  „. 

1.  Case  1:        0  <  ep  <     tt/2 

The  barrier   is  lowered  for     -  tt/2  <  <(>  <  tt/2. 

and  for      tt/2  <  <J>    <3  tt/2,        0  <  6<   tan-1  [-cot  (9p)  /cos  (<}>)] . 
The  normalized  field  enhanced  emission  rate  for  this  case  is  given 

by 

tt/2         tt/2 
en3/eno  =   (1/2  tt  )   j  d^    ;     sin(e)    exp[-  g  a  Vtj  (q  ,$)  ]de 

-tt/2        0 

3 tt/2       tan-   [-cot  (e  F)  /cos  ($  )  ] 
+   (1/2tt  )     r       d  4>    j  '  sin(e)exp[-gAEti  (6,(j))]de 

tt/2         0 

377/2  2      1/7 

+     (1/2  tt);  {1    +     [cot(QF)/cos(<p)]z}~l/z    d|. 

tt/2 

(4. 47. a) 

2.  Case  2:      tt/2  <9p<  tt 

The  barrier   is   increased  for  tt/2  <  <}>    <3  tt/2. 

and  for  -  tt/2  <  <f>  <  tt/2,  0  <  6  <  tan-1  [-cot  (3  F) /cos  #  )  ] . 
The  barrier  is  lowered  for  -tt/2  <<r)  <tt/2  ,  tan-1  [-cot (6p) /cos ($)] 
<9  <  tt/2. 

The  normalized  field  enhanced  emission  rate  for  this  case  is 
given  by 

tt  /2 

{] 

-TT    /2 


en3/eno  =  1  -    (1/2^  )  j         {1  +    [cot(eF)/cos(^)]2}-1/2  d<| 


tt/2        tt/2 
+    (1/2tt);  d*  sin(6)exp[-  BAEti(6,<j))]de 

-tt/2       tan~1[-cot(eF)/cos(*)]  (4.47.b) 
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For   &p  =  0,   the  e  dependence  disappears,  and  the  problem  can  be  solved 
analytically.     The  barrier  lowering  simplifies  to 

Efci    =    -1.9    (p/4  tt   E(f  s)1/3    F2/3    cos(0)  (4.48) 

and  the  normalized   Poole-Frenkel   effect  emission   rate  enhancement   is 
given  by 

en3/eno  =    (1/R)    x    [exp(R)  -  1]  (4.49) 

where 

R  =  1.9   (p^^^eg)1/3  F2/3  /kT 

4.3.5.2     Phonon-assi sted  tunneling  effect 

For  the  potential  wells  studied  here,  the  one-dimensional 
analysis  was  done  for  9  =  0,  where  the  Poole-Frenkel  barrier  lowering 
is  dominant.  For  the  dipole  well,  the  maximum  barrier  lowering  occurs 
at  e  =  0m,  where  0  <  em  <8F  and  <j>  =  0.  Now,  differentiating  Eq. 
(4.46)  for  AEt;j  with  respect  to  9  at  <j>  =0,  and  setting  the  result 
equal  to  zero,   we  obtain  a  solution  for  9m 

0  m  =   tan"1    [    {[8  +  cos2(6F)]1/2 

-  3  cos(0F)}   /  2sin(0F)]  (4.50) 

Substituting      Eq.     (4.50)       into    Eq.     (4.44)       and      Eq.     (4.44)       into 
Eq.    (4.13),    and    letting    z    =    (Et;j    -   Efch)    /   kT,    the    field    enhanced 
emission  rate  due  to  phonon-assi sted  tunneling  effect  can  be  expressed 
as 
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(Eti-AEtj)/kT  r2 

ent/eno  =j  exp{z  -    [(8m)1/z/M]j       [-qPcos(6m) 

0  rx 

/(4^eQesr2)-qrF  sj^e^sinie^    -   qrF  cos^  cos  (9^ 

+   zkT]1/2  dr}  dz  (4.51) 

Here  r1  =   [  (qPcos6m  /  (47r^eszkT)  ] 1//2     and     r2  =  zkT  /  {qF   [  sin(^) 
sin(6m)   +  cos  (8F)  cos(9m)]  }.     The  normalized  field  enhanced  emission 
rate  due  to  phonon-assisted   tunneling  effect  and   Poole-Frenkel   effect 
can  be  calculated  by     numerical     method.     The  results  are  shown  in 
Fig.   4.7.     Thus,    the  total   field  enhanced  emission  rate  is 

enHD  =  eno    <en3  +  ^  (4-52> 

Table  4.1  summarizes  the  potential  wells  for  Colombic  well,  Dirac 
well,  square  well,  polarization  well  and  dipole  well.  Table  4.2 
summarizes  the  one-dimensional  and  three-dimensional  emission  rate  due 
to  Poole-Frenkel  effect  for  five  different  potential  wells.  Figure  4.8 
shows  the  emission  rate  vs  electric  field  for  the  Dirac  well,  square 
well,  polarization  well  and  coulombic  well  as  compared  with  the  zero 
electric  field  values.  The  enhanced  emission  rate  for  the  neutral 
trap  with  Dirac  well,  square  well  or  polarization  well,  has  almost  the 
same  value.  However,  the  enhanced  emission  rate  for  the  Coulombic 
well  depends  very  strongly  on  the  electric  field. 


4.4     Theoretical  calculations  of  the  nonexponential  DLTS 
transient  for  different  potential  wells 


From  the  analysis  of    field   enhanced   emission   rate,    the  electric 
field  dependent  DLTS  response  can  be  calculated.      This   is  discussed  as 
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Figure  4.7  Normalized  enhanced  emission  rate  vs.  electric  field  for 

EL2a  electron  trap  in  GaAs  for  phonon-assisted  tunneling  and 
Poole-Frenkel  effects. 
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Figure  4.8  Enhanced  emission  rate  vs.  electric  field  for  EL2a  electron 
trap  in  GaAs  for  Coulombic  well,  Dirac  well,  polarization 
well ,  and  dipole  well . 
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follows:  consider  a  p+n  abrupt  junction  diode,  the  electric  field  is 
in  general  spatial  dependent  within  the  depletion  region,  and  can  be 
expressed  by 

F  =  Fmax(l  -  x/W)  (4.53) 

where  x  is  the  distance  from  the  junction;  Fmax  is  the  maximum 
electric  field  occurring  at  the  metalallurgical  junction  of  the  p+n 
diode  and  is  given  by 

Fmax  =  <3Ndw  /=oes  <4-54> 

where  W  is  the  depletion  width  under  reverse  bias  condition. 

w=    [2£o£s   <Vbi    +  V   /  ^Nd^/2 

=  LD[2(Vbi+Vr)AT  -  2]1/2  (4.55) 

where  Lp.   is  the  extrinsic  Debye  length  and  is  given  by 

Lp  =  [eoeskT  /  (q2ND)]1/2  (4.56) 

where  V^j  is  the  built-in  potential,  which  reads 

Vbj  =  (kT/q)  lntNj/np  (4.57) 

where  n^  is  the  intrinsic  carrier  density, 

nd=  (NCNV)1/2  exp(-Eg/2kT) 

=   4.9x  1015(mdemdh/mo2)3/4T3/2  exp(-Eg/2kT)     (4.58) 

and  Eq  is  the  energy  bandgap,   which  is  a  function  of  temperature.    *■     ■■ 

E    (T)    =  E    (0)    -Yl  T2  /    (T  +Y2)  (4.59) 
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where  E  (0)  =  1.519eV,  y  1  =  5.405  x  10~4  eV/K,  and  y  2  =  204  for  GaAs. 
Electrons  which  are  located  in  the  region  between  WQ  (i.e.,  zero  bias 
depletion  layer  width)  and  W  (i.e.,  reverse  bias  depletion  layer 
width)  will  be  emitted  into  the  conduction  band  when  an  applied  bias 
is  increased  from  0  to  -V"r,  as  is  shown  in  Fig.  4.9.  The  depletion 
width  in  the  junction  space  charge  region  is  first  divided  into 
equally  spaced  small  segments.  In  each  small  segment,  the  electron 
trap  density  and  electric  field  are  assumed  constant.  For  each 
electric  field  (Fj)  strength,  there  is  a  corresponding  enhanced 
emission  rate,  en:j(Fj).  If  one  assumes  that  the  electron  emission  is 
exponentially  transient  within  each  segment,  namely,  exp[e  j  (Fi)]  t, 
then  the  total  emission  transient  in  the  depletion  width  is  equal  to 
the  sum  of  the  individual  components,   which  read 

I  exp[-eni(Fj)t]  (4.60) 

Therefore,  the  DLTS  signal  can  be  expressed  by 

S(x)  =Z  expt-e^F^]  -  E  exp[-eni  (Fj)  t2]        (4.61) 

where  en^  is  the  electron  emission  rate  within  each  small  segment  in 
the  depletion  region,  as  is  shown  in  Eq.  (4.14).  Values  of  e  :  can  be 
calculated  from  Eqs.  (4.26,  4.31,  4.35,  4.42,  and  4.52)  for  different 
potential  wells.  From  Eq.  (4.61),  it  is  noted  that  the  DLTS  spectral 
response  is  nonexponential.  Figure  4.10  illustrates  the  DLTS  response 
for  the  EL2a  (Ec-0.83eV)  electron  trap  in  GaAs  calculated  from  Eq. 
(4.61)  for  different  potential  wells  and  for  zero  electric  field.  The 
theoretical  calculations  showed  that  the  location  of  DLTS  signal  peak 
for  attractive  or  neutral  trap  is  different  along  the  temperature 
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Figure  4.10  DLTS  response  for  the  EL2a  level  for  Coulombic  well 
Dirac  well,  with  zero  electric  field. 
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axis.  For  Coulombic  potential  well,  the  DLTS  signal  peak  will  occur 
at  the  lowest  temperature,  while  for  the  neutral  trap  with  Dirac  well, 
square  well,  and  polarization  well  will  produce  identical  DLTS  signals 
at  the  same  temperature.  Figure  4.11  shows  the  DLTS  spectral  response 
for  the  EL2a  electron  trap  for  the  case  of  Coulombic  well  with 
single-  and  double-charge  state  along  with  the  experimental  results. 
Figure  4.12  illustates  the  DLTS  signal  for  the  EL2a  electron  trap  in 
GaAs  for  Coulombic  well  with  a  double  charge  states  for  different 
window  rates.  Compared  the  theoretical  with  experiment  DLTS  spectral 
signal  will  enable  us  to  determine  the  potential  well  for  the  EL2a 
electron  trap  level.  Our  results  showed  that  the  most  probable 
potential  well  for  the  EL2a  trap  is  the  Coulombic  well  with  a  double 
charge  states. 

4.5  Summary  and  Conclusions 

We  have  analyzed  the  nonexponential  capacitance  transients  due  to 
electric  field  dependent  emission  rate  of  trapped  charge  of  deep- level 
traps  in  GaAs.  A  comparison  of  the  theoretical  calculations  of  the 
nonexponential  DLTS  response  with  the  DLTS  data  for  each  trap  level 
would  allow  us  to  determine  the  potential  well  of  the  trap  involved. 
This  is  one  of  the  methods  which  may  be  used  to  identify  the  physical 
origins  of  the  trap  in  a  semiconductor.  From  the  theoretical  analysis 
of  DLTS  response,  the  conclusions  are  given  as  follows: 
(1).  The  location  of  DLTS  signal  peak  for  attractive  or  neutral 
trap  is  different  vs  temperature.  For  Coulombic  well,  the 
DLTS  signal  peak  will  occur  at  the  lowest  temperature,   while 
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Figure  4.11  Nonexponential  DLTS  response  for  the  EL2a  trap  for 
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for  neutral  traps  such  as  Dirac  well,  square  well, 
polarization  well,  and  dipole  well,  the  DLTS  signals  have  almost 
the  same  shape  and  the  DLTS  peak  will  occur  at  the  same 
temperature . 

(2).  For  Coulombic  well,  the  DLTS  signal  peak  for  the  double  charge 
state  will  occur  at  a  lower  temperature  than  that  of  the  single 
charge  state. 

(3).  The  shape  of  the  DLTS  signal  is  affected  by  the  capture  cross 
section.  For  large  capture  cross  section,  the  narrow  shape  of 
the  DLTS  response  will  occur  in  the  low  temperature  region. 
For  small  capture  cross  section,  the  shape  of  DLTS  response 
will  become  broader  and  move  to  the  higher     temperature  region. 

(4).  The  DLTS  response  is  proportional  to  the  background 
concentration,  and  reverse  bias.  Higher  background 
concentration  and  higher  reverse  bias  will  result  in  a  bigger 
emission  rate  enhanced  due  to  high  electric  field     dependence. 

(5).  Nonexponential  capacitance  transient  due  to  electric  field 
dependence  is  expected  in  the  high  reverse  bias  condition.  Only 
in  the  very  low  electric  field  (<103  V/cm)  or  in  the  uniform 
electric  field     condition  will  DLTS  response  be  exponential. 


V     MODELLING  OF  GROWN-IN  NATIVE  POINT  DEFECTS    IN  GaAs 

Theoretical  modelling  of  grown-in  native  point  defects  in  GaAs  is 
presented  in  this  chapter.  Based  on  chemi cal- thermodynamic 
principles,  expressions  for  the  equilibrium  defect  concentration  as 
functions  of  temperature  and  arsenic  pressure  during  crystal  growth 
are  derived.  Thermal  kinetic  equations  are  then  employed  to  predict 
the  possible  native  defects  in  GaAs  after  crystal  cooling.  For  GaAs 
grown  under  As-rich  or  high  arsenic  pressure  condition,  it  is  shown 
that  several  native  defects  such  as  gallium  vacancy  (VG  ) ,  arsenic 
interstitial  (As.j),  arsenic  antisite  (AsGa) ,  arsenic  anti si te-plus- 
arsenic  vacancy  (v"AsAsGa)  and  their  complexes  may  be  observed  in  GaAs 
material. 

5.1     Introduction 

GaAs  specimens  grown  by  various  techniques  such  as  Liquid 
Encapsulation  Czockraski  (LEC) ,  Vapor  Phase  Epitaxy  (VPE) ,  Liquid 
Phase  Epitaxy  (LPE),  Metalorganic  Chemical  Vapor  Deposition  (MOCVD) , 
and  Molecular  Beam  Epitaxy  (MBE)  are  known  likely  to  produce  different 
defect  properties  in  terms  of  energy  levels  within  the  bandgap  of  a 
semiconductor  due  to  differences  in  the  native  defects  and  trace 
impurities  or  impurity  complexes.  In  each  of  these  techniques,  the 
growth  temperature,  growth  pressure,  growth  phase,  and  cooling  rate 
are  usually  different.  For  example,  GaAs  grown  by  LPE  technique  from 
a  gallium  melt  is  expected  to  be  low  in  gallium  vacancy  defects  and 
high  in  arsenic  vacancy.     The  LPE  GaAs  usually  contains  hole  traps 
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with  energy  level  of  Ey  +  0.71  eV  (B-center) ,  while  the  VPE  GaAs  grown 
in  high  arsenic  pressure  conditions  always  contains  an  electron  trap 
with  energy  of  Ec  -  0.83  eV  (EL2a  level).  The  reason  for  observing 
these  trap  levels  in  LPE  or  VPE  GaAs   is  still  not  clear. 

In  this  chapter,  we  shall  introduce  the  chemical-thermodynamic 
principles  and  thermal  kinetic  equations  to  interpret  grown-in  defects 
in  GaAs.  The  use  of  chemical-thermodynamic  principles  to  analyze 
point  defects  in  GaAs  was  proposed  first  by  Logan  and  Hurle. LJ  They 
considered  the  shallow  vacancy  as  well  as  shallow  interstitial  levels. 
Bublik  *■  *  modified  their  model,  and  applied  it  to  the  deep  vacancy 
level.  More  recently,  Hurle  *■  *  used  Frenkel  defects  instead  of 
Schottky  defects  for  the  arsenic  vacancy  and  arsenic  interstitial 
defects.  During  crystal  growth,  both  VGa  and  VAs  Schottky  vacancy 
pairs  are  important  defects  for  the  undoped  GaAs.  In  the  present 
work,  we  use  VGa  and  V\s  Schottky  vacancy  pairs  with  As-Frenkel 
defect  to  calculate  the  vacancy  and  interstitial  concentrations  during 
crystal  growth.  During  crystal  cooling  from  high  temperature  to  room 
temperature,  thermal  kinetic  is  expected  to  be  the  dominant  mechanism. 
Since  the  formation  energy  for  antisite  defect  is  very  low  in  the 
after-growth  condition  (V"Ga  +  ASj  =  AsGa,  a  H  =  0.35  eV) ,  the  antisite 
defect  is  likely  to  be  the  dominant  defect  after  crystal  growth,  as 


was  proposed  by  Van  Vechten. 
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From  the  theoretical  analysis  of  GaAs  material  grown  under  As- 
rich  or  high  arsenic  pressure  condition  [such  as  LEG,  VPE,  and  MOCVD 
techniques],  defects  such  as  arsenic  interstitial  (ASj),  gallium 
vacancy  (VjO,  arsenic  antisite  (AsGa),  arsenic  antisi te-plus-arsenic 
vacancy   (VacASqJ,    and  their  complexes  may  be  the  dominant  defects  in 
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undoped  GaAs,  while,  gallium  vacancy  (Vq  ),  arsenic  vacancy  (V^g)  and 
their  comlexes  may  be  the  dominant  grown  -in  defects  in  MBE  and  LPE 
GaAs  under   low  arsenic  pressure  conditions. 

Section  5.2  explains  the  theoretical  calculations  of  vacancy  and 
interstitial  defects  in  GaAs  during  crystal  growth.  Thermal  kinetic 
equations  are  described  in  section  5.3.  Section  5.4  discusses  the 
possible  grown-in  point  defects  in  GaAs  for  the  As-rich  or  high 
arsenic  pressure  cases.  The  possible  grown-in  point  defects  in  GaAs 
for  the  Ga-rich  or  low  arsenic  pressure  cases  is  depicted  in 
section  5.5.      Summary  and  conclusions  are  given   in  section  5.4. 

5.2     Theoretical  Calculations  of  Vacancy  and  Interstitial 
in  Undoped  GaAs 

In  this  section,  the  chemical -thermodynamic  principles  L&^-6yJ 
are  used  to  derive  expressions  for  thermal  equilibrium  defect 
concentration  as  functions  of  temperature  and  arsenic  pressure  in 
GaAs.  Hurle  '■°^*  calculated  the  point  defects  by  considering  only  As- 
vacancy  and  As-interstitial  defects  in  undoped  GaAs.  However,  it  is 
known  that  Ga-vacancy  related  defects  may  play  an  important  role  in 
antisite   formation. 

In  this  chapter,  we  consider  several  types  of  point  defects  such 
as  arsenic  monovacancy,  V*  ,  positively  charged  arsenic  monovacancy 
V"as  ,  arsenic  interstitial  As^,  positively  charged  arsenic 
interstitial  As^  ,  gallium  monovacancy  V"Ga,  and  negatively  charged 
gallium  monovacancy  VG~. 

To  deal  with  the  problem  of  native  point  defects  in  GaAs,  one  can 
write  down   the   reaction  equations   for   formation   of  each   type   of 
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defects  as  well  as  formation  of  electrons  and  holes  in  the  crystal. 
There  is  an  additional  reaction  equation  which  represents  the  transfer 
of  atoms  between  gas  and  solid  phases.  For  each  of  these  reactions, 
there  is  a  mass  action  which  applies  in  equilibrium.  The  mass  action 
can  be  written  in  terms  of  concentrations,  and  Boltzmann  statistics  is 
used  for  electrons.  To  these  mass  action  equations  one  adds  the 
condition  of  charge-neutrality,  and  the  resulting  set  of  equations  can 
then  be  solved  as  a  function  of  arsenic  partial  pressure.  Considering 
the  defects  cited  above,  one  has  the  following  reaction  equations. 

ASAS  +  Vj  =  ASi  +  VAs  (5.1) 

VAs    =   VAS+    +  e  "  (5.2) 

Asi    =   ASj"1"   +   e"  (5.3) 

0  =  e~  +  h+  (5.4) 

(l/2)As2(g)    +  Vd   =  Asd                                                             (5.5) 

0  =  VGa  +  VAs  (5.6) 

VGa  =  vGa~  +  h+  (5.7) 

Equations  (5.1)  to  (5.4)  represent  the  reactions  for  forming  the 
ionized  arsenic  Frenkel  defects.  Equation  (5.5)  denotes  the  transfer 
of  arsenic  atoms  between  solid  and  gas  phases.  Equation  (5.6)  shows 
the  formation  of  Schottky  pairs.  Equation  (5.7)  represents  the 
ionization  of  a  Ga-vacancy.  The  mass  action  relationships 
corresponding  to  the  above  reactions  are  given  as  follows    [62,63] 
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Kfa  =    [ASj]     [VAS] 

=  2.92  x  106  exp(-4.845AT)  (5.8) 

Kav  =  n  tVAs+^  /  tVAsl 

=  442  exp(-0.27AT)  (5.9) 

Kai   =  n   [vi+]   /   [As0 

=  4.9  x  109   x  T3/2  exp(-0.4/kT)  (5.10) 

Kcv  =  "P 

=  1  x  10"12  T3  exp(-1.62/kT)  (5.11) 

KAs2i    ■      ^    PAS2_1/2 

=  16.4  exp(-1.125AT)  (5.12) 

Ks  =    ^1    tVAsl 

=  3.286xl04  exp(-3.6AT)  (5.13) 

Kgv  "  P   tVGa"J   /   tVGa^ 

=  3.7  x  10~8  T3/2  exp(-0.66AT)  (5.14) 

Square  brackets  in  the  above  equations  are  used  to  indicate  the 
concentration.  The  AsAs  represents  arsenic  atom  at  the  arsenic  site, 
and  is  taken  as  unity.  The  partial  pressure  of  AS2  in  the  gas  phase 
is  denoted  by  PAs2-  The  equilibrium  constants  appearing  on  the  right 
handside  of  Eqs.  (5.8)  through  (5.14)  have  the  general  form. 

K  =  exp(ASA)  exp(-AHAT)  =  KQ  exp(-AH/kT)        (5.15) 

where  AS  and  AH  are  the  entropy  and  enthalpy  changes  for  each 
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reaction,  respectively.  The  charge  neutrality  condition  is  obtained 
from  Poisson  equation 

n  +     [VGa~]    =  P  +     [VAs+]    +   [Asi+]  (5'16) 

Equations  (5.8)  through  (5.16)  are  usually  solved  by  Brouwer 
approximation.    *■'    *      In   the  present  case,    we  solved  Eqs.    (5.8)    to 

9 

(5.16)   directly  and  obtained  an  expression   for  n^  as 

n2  =   {Kcv  +    (KavKfa/KAs2j)    PAs2-1/2     +  K^K^P^1/2} 

/   (1  +  KqvKsKAs2iPAs21/2  '  KfaKcv>  <5-17> 

Thus,  the  electron  concentration  can  be  calculated  from  Eq.  (5.17)  and 
other  defect  density  can  also  be  deduced  from  n  via  Eqs.  (5.8)  to 
(5.16).  Figure  5.1  shows  the  defect  concentration  vs.  PAs2  f°r  T  = 
1000  K.  Note  that  Ga  vacancy  is  the  dominant  point  defect  in  the 
entire  As2  pressure  range  shown.     The  LPE  GaAs  is  usually  grown  in  the 

Q 

lower  As2  pressure  range  with  PAs2  =  10  atm,  and  defects  such  as 
VG  ,  VGa~,  and  V"As+  are  the  dominant  defects.  For  VPE  GaAs  with  a 
corresponding  partial  pressure  of  10  atm,  defects  such  as  VQg,  Vq3-, 
VAs  ,  ASj  ,  and  As^  are  the  dominant  defects.  In  the  pressure  range 
between  LPE  and  VPE  growth,  the  concentration  of  VGa~  is  equal  to  the 
concentration  of  V\,  +.  Figure  5.1  also  shows  that  the  grown-in 
defects  for  the  VPE  and  MOCVD  techniques  are  much  more  than  LPE  and 
MBE  growth  techniques.  Figure  5.2  shows  defect  density  as  a  function 
of  growth  temperature  for  PAs2  =  5  x  10  J  atm.  The  concentrati  on  of 
V^-  monovacancy  and  VAs+  monovacancy  increases  with  increasing  growth 
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Figure  5.1     Defect  density  vs.  As  partial   pressure  for  GaAs  grown 
at  1000°K. 
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Figure  5.2  Defect  density  vs.  temperature  for  GaAs  grown  at  an 
arsenic  pressure  of  5x1 0"3  atm.. 
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temperature,    but  the  density  of  ASj+  mono-intersti tial  decreases  with 
increasing  growth  temperature. 

5.3     Thermal  Kinetic  after  Crystal  Growth 

In  general,  defects  will  migrate  as  the  crystal  cools  down  from 
the  growth  temperature.  Thus,  one  would  expect  defect  concentrations 
to  reach  a  new  equilibrium  condition  at  lower  temperature.  [64,71] 
The  enthalpy  of  a  single  vacancy  migration  in  GaAs  is  1.6  eV.  [64,72] 
For  examples,  if  the  jump  attempt  frequency  is  equal  to  Debye 
frequency,  then  the  jump  rate  will  be  4  x  107  s-1  at  1400  K  and 
5  x  104  s_1  at  1000  K.  For  example,  for  a  typical  crystal  to  cool 
down  from  1050  K  to  950  K  in  10  min.,  simple  vacancy  would  pass 
through  more  than  108  lattice  sites  in  that  period.  The  migration  of 
vacancy  will  produce  antisite  defects,  antisite  pairs  or  impurity 
complexes.  In  LPE  grown  GaAs  samples,  we  have  observed  that  samples 
with  faster  cooling  rate  (1°  C/min.)  would  produce  more  defects  than 
slower  cooling  rate  (0.4°  C/min.).  Most  of  the  simple  vacancies  are 
ionized,  and  may  encounter  with  other  defects  to  form  bound  complexes. 
At  room  temperature,  almost  all  the  vacancies  present  should  be  tied 
up  with  other  defects  to  form  complexes.  To  consider  the  ultimate 
fate  of  a  simple  point  defect  introduced  during  crystal  growth, 
several  types  of  defects  and  defect  complexes  should  be  considered. 
It  seems  likely  that  the  most  common  defect  complexes  will  be  those 
which  have  no  net  charge  because  these  would  have  the  greater  binding 
energy.      For  example,    a  negatively  charged  gallium   vacancy  and  a 


78 


positively  charged  arsenic  vacancy  may  produce  a  neutral  gallium- 
arsenic  di vacancy.     This  can  be  expressed  by 

VGa~  +  VAs+  "  VGa'VAs+  (5.18) 

In  addition,  the  anti  si  te-di  vacancy  complexes  may  be  formed  by  the 
reaction  equations  shown  below 

2VAS+  +  ^As"  =  VAs+GaAs~VAs+  (5.19) 

and 

2VGa~  +  AsGa++  =  VGa~AsGa++VGa~  (5-20) 

The  anti  site  pair  defect  complexes  can  be  written  as 

^As"  +  AsGa++  =  GaAs"AsGa++  (5-21) 

The   interaction  of  VAs+    (VGa~)    with  a   single   acceptor    (or   donor) 
impurity,   A~   (D+),   can  be  written  as 

A~  +  VAs+  =  AAs  (5-22) 

7Ga~  =  DGa  (5.23) 

which  should  have  about  the  same  binding  energy  as  the  gallium-arsenic 
di  vacancy  complex   given  by   Eq.     (5.18). 

From  Fig.  5.1,  it  is  noted  that  the  concentration  of  VG     is  very 
high,  and  the  following  reactions  are  prevailed 

VAs+  +  VGa  +  e"  "  VGaVAs  (5-24) 

vGa  +  D+  =  DGa+  (5.25) 
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In  general,  the  concentration  of  positively  charged  and 
negatively  charged  simple  point  defects  introduced  during  crystal 
growth  are  not  equal.  Therefore,  the  concentration  of  one  type  of 
simple  defects  may  be  exhausted  as  shown  by  the  reaction  Eqs.  (5.18) 
to  (5.25).  This  can  be  illustrated  by  the  following  kinetic  equations. 
(1).  Arsenic  antisite  can  be  formed  by  ASj    ,   As^ ,   V>.a~,   and  V^. 

Asi+  +  VGa"  =  AsGa++  +  2e~  (5'26> 

Asi+  +  VGa  =  ^Ga**  +  G"  <5'27) 

ASj    +  VGa~  =  AsGa++  +  3e"  (5.28) 

Asi  +  vGa  =  ^Ga**  +  2e"  <5-29> 

(2) .  Anti  si  te-vacancy  complexes  can  be  expressed  by 

VGa"  +  AsAs  =  AsGaVAs+++  +  4e~  <5-30> 


'Ga  T  "aAs  *  rtaGavAs 


(5.31) 

W  +  ^Ga  +  4e~  =  ^As^a""  <5'32) 

which  will  occur  everytime  a  vacancy  migrates  to  a  nearest  neighbor 

site. 

(3).  Impurity- vacancy  complexes  can  be  written  as 

VGa"  +  AAs  =  AAsVGa~  (5'33> 

and 

VGa"  +  DGa+  +  e"  =  DGavGa~  (5.34) 


Equations  (5.33)  and  (5.34)  are  the  interaction  of  VG~  with  a 
acceptor  (donor)  substitutional  sites.  For  examples,  TeA  VGa~  defect 
complex  may  be  formed  in  Te-doped  GaAs,  while,  SnGavGa~  and  ^Ga^Ga" 
defect  complexes  may  be  formed  in  Sn-doped  and  Ge-doped  GaAs, 
respectively. 

The  VAs+  may  migrate  to  a  donor  substitutional  site,  and  the 
reaction  equations  is  given  by 

VAs+  +  DGa  "  DGaVAs+  <5'35> 

Obviously,  such  reactions  act  as  to  compensate  the  dopant  and  to 
prevent  the  Fermi -level  from  approaching  either  band  edge  during 
crystal  cooling. 

5.4     The  Possible  Grown-in  Point  Defects  in  GaAs  for  the 
As-rich  or  High  Arsenic  Pressure  Case 

Applying  defect  modelling  to  the  high  arsenic  pressure 
conditions,  the  V"Ga,  V"Ga~,  VAg  ,  ASj  +  ,  and  As^  are  found  to  be  the 
dominant  defects  when  PAs2  *s  greater  than  10  atm  [such  as  LEC,  VPE, 
and  MOCVD  techniques];  this  is  shown  in  Fig.  5.1.  The  migration  of 
VG  ,  V"Ga~,  and  V"As  defects  may  result  in  forming  complexes  given  by 
Eqs.  (5.18)  to  (5.20)  and  Eqs.  (5.22)  to  (5.35)  after  crystal  growth. 
If   we  neglect   the  neutral   defect  complexes,    then  native  defects  such 

as  AsGaVAs+++'  AsGa++'  vGa"'  GaAsVGa  '  VAs+'  Asi+'  VGa'  Asi '  and 
impurity    complexes    DGa+,    AAsVGa~,    DGaV"Ga~,    DGaV"As+    are    the   possible 

defects  which  may  be  performed   in  high  arsenic  pressure  conditions. 

In  As-rich  condition,   the  concentration  of  VA_+  may  be  decreased,  As.j  + 
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may  be  increased,  and  AsGaVAs+++,  AsGa++,  VGa~,  ASj+,  VGa,  As^,  and 
impurity  complexes  are  the  possible  defects. 

5.5     The  Possible  Grown- in  Point  Defects  in  GaAs  for  the 
Ga-rich  or  Low  Arsenic  Pressure  Case 

Apply  defect  modelling  to  the  low  arsenic  pressure  conditons. 
VG   ,  VGa~,  and  VAs     are  found  to  be     the  dominant  point  defects  when 

PAs2  -"s  less  tnan  10  atK1'  tn:"s  ^s  the  case  fot  the  LPE  and  MBE  grown 
GaAs.  The  vacancies  migrate  to  form  complexes  given  by  Eqs.  (5.18)  to 
(5.20),  Eqs.  (5.22)  to  (5.24),  and  Eqs.  (5.30)  to  (5.35)  after  crystal 
growth.  If  we  neglect  the  neutral  defect  complexes,  then  native 
defects   such  as  VGa",    VAs+,    GaAsVGa       ,    AsGaVAs+++,    VGa   and    impurity 

complexes  such  as:  DGa+,  AAsvGa~'  DGavGa~'  DGaVAs+  are  the  P°ss^ble 
defects  in  low  arsenic  pressure  conditions.  The  concentration  of  VGa~ 
is  almost  equal  to  that  of  VA  ,  as  is  shown  in  Fig.  5.1,  and  the 
binding  energy  for  the  gallium-arsenic  divacancy  is  very  high,  as 
given  in  Eq.  (5.18).  Thus,  high  purity  GaAs  can  be  grown  in  this  case 
under  optimum  cooling  condition.  In  Ga-rich  and  low  arsenic  pressure 
conditions,  the  concentration  of  VG~  and  VGa  may  be  decreased, 
gallium  interstitial  may  be  react  with  VAs  to  form  GaA  .  Thus  the 
possible  defects  are  VAs+,   GaAs      ,   and  GaAgVGa  in  undoped  GaAs 

grown  under  Ga-rich  condition. 

5.6     Summary  and  Conclusions 

In  this  chapter,    we  presented  a  new  defect  model   for   predicting 
native  point  defects  in  GaAs.     During  crystal  growth,  the  chemical- 
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thermodynamic  principles  are  used  to  derive  the  density  of  vacancy  and 

interstitial  defects  under   thermal  equilibrium  condition.      After 

crystal  growth,    thermal  kinetic  equations  are  employed  to  predict  the 

antisite,    complexes,    and   impurity  complexes.     Conclusions   are    listed 

as  follows 

(1).  High  purity  GaAs  material  can  be  grown  for  the  low  arsenic 
pressure  case  under  optimum  cooling  condition. 

(2).  GaAs  grown  under  higher  arsenic  pressure  condition  will  produce 
more  native  point  defects  than  under  lower  arsenic  pressure 
condition. 

(3).  Native  defects  such  as  AsGa,  AsGaV"As,  V"G  ,  and  As^  and  impurity 
complexes  such  as  DGa+,  AAsVGa~,  DGaVGa~  are  the  possible  defects 
for  GaAs  layers  grown  under  the  As-rich  or  high  arsenic  pressure 
condition. 

(4).  Native  defects  such  as  VGa,  VAs,  GaAsVGa,  AsGgVAs  and  impurity 
complexes  such  as  DGa,  AAsVGa,  DGaVGa,  DGaV"As  are  the  possible 
defects  for  GaAs  layers  grown  under  low  arsenic  pressure  case. 

(5).  Native  defects  such  as  V"As,  GaAg,  Ga^V^  and  impurity  complexes 
such  as  DGaV"Ag+  are  the  possible  defects  for  GaAs  layers  grown 
under  Ga-rich  and  low  arsenic  pressure  conditions. 

(6).  Arsenic  antisite  (AsGa)  defect  only  observed  in  GaAs  grown  under 
As-rich  or  high  arsenic  pressure  conditions;  this  defect  can  not 
be  produced  under  low  arsenic  pressure  and  Ga-rich  conditions. 


VI.      ON  THE  PHYSICAL  ORIGINS  OF  EL2  ELECTRON  TRAP  IN  GaAs 
(THEORETICAL  AND  EXPERIMENTAL  EVIDENCE) 

The  activation  energy  for  the  EL2  electron  trap  in  GaAs  reported 
in  the  literature  ranges  from  Ec-0.76eV  to  Ec-0.83eV.  The  physical 
origin  of  this  trap  is  a  subject  of  great  interests  in  recent  years. 
A  large  number  of  papers  has  been  devoted  to  this  subject.  Based  on 
our  theoretical  model  and  experimental  results,  we  found  that  EL2 
level  is  formed  by  two  electron  traps.  One  is  identified  as  the  EL2a 
(Ec-0.83eV)  electron  trap,  the  other  is  denoted  as  the  EL2b  (E  - 
0.76eV).  The  physical  origin  of  EL2a  level  is  due  to  the  arsenic 
antisite  (AsGa)  defect,  whereas,  the  physical  origin  of  the  EL2b 
level  is  due  to  arsenic  anti si te-plus-arsenic  vacancy  (AsG  V,  ) 
complex . 

Based  on  the  native  defect  modelling  presented  in  last  chapter, 
modelling  of  EL2  electron  trap  is  derived  in  GaAs  material  grown  under 
As-rich  or  high  arsenic  pressure  case.  Furthermore,  from  calculations 
of  the  trap  density  for  the  MOCVD  and  VPE  grown  GaAs  with  As  to  Ga 
mole  fraction  ratio  greater  than  one,  it  is  found  that  the  density  of 
EL2a  trap  is  proportional  to  the  mole  fraction  ratio  of  (r-l)1//2; 
whereas,  the  density  of  EL2b  trap  is  proportional  to  (r-l)1//4,  where 
r   =    [As]/[Ga]. 

The  EL2a  trap  is  best  fitted  to  a  Coulombic  potential  well  with 
double  charge  state  (i.e.,  AsGa++) ,  this  was  verified  by  the  Poole- 
Frenkel  effect  and  phonon-assisted  tunneling  effect  as  observed  in  the 
field  dependent  electron  emission  rates  using  nonexponential  DLTS 
method. 


83 


84 


DLTS  measurements  were  performed  on  VPE  GaAs  layers  grown  on 
different  orientations  [e.g.,  (100),  (211A),  (211B)]  and  LEC  GaAs 
layers  annealed  in  hydrogen  gas  at  different  temperatures  [e.g.,  200, 
300,  500°C].  From  the  results  of  our  annealing  study,  it  was  found 
that  EL2b  trap  would  disappear,  and  EL2a  trap  would  emerge  in  the  DLTS 
scan  at  annealing  temperature  of  500°C.  This  result  may  be 
interpreted  by  the  model  of  EL2  electron  trap  to  be  presented  in  this 
chapter.  The  experimental  evidence  which  may  be  used  to  support  the 
modelling  of  EL2  will  be  presented  in  chapter  VII. 

Section  6.1  review  the  EL2  electron  trap  in  GaAs.  Theoretical 
model  of  EL2  electron  trap  is  explained  in  section  6.2.  Section  6.3 
discusses  the  method  of  determining  the  potential  wells  for  EL2  trap 
levels  from  analyzing  the  field  enhanced  emission  rates  in  the 
nonexponenti al  DLTS  experiment.  Summary  and  conclusions  are  given  in 
section  6.4. 

6.1.     Review  of  the  EL2  Electron  Trap  in  GaAs 

Activation  energy  of  EL2  electron  trap  covers  the  energy  range 
from  E  -0.75eV  to  Ec-0.83eV  as  reported  by  many  previous 
investigators.    l'3-91J      This   trap  has  been   observed   in  GaAs  grown  by 

the  Bridgmann,    [77'781    LEC,    [79~85]   VPE,    I86'87!    and  MOCVD   t88"91] 

f921 
methods   as   well   as   high   temperature   heat-treated   GaAs    samples.     L     J 

However,    this   level   was  not  observed  in  the  LPE    ^93^    and  MBE    '      > 

grown  GaAs  epitaxial  materials.     Recent  studies  of  the  LEC  bulk  grown 

GaAs    reported    by  Taniguchi    et   al.     *■      *    have    found    that   Ec-0.77eV 

electron  trap  exists   in  the  front  section  of  LEC  GaAs   ingot;   whereas, 
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Ec-0.82eV  level  was  the  dominant  trap  level  appears  in  tail  section  of 
LEC  GaAs  ingot.  The  same  result  was  reported  in  MOCVD  grown  GaAs 
epi  layers  by  Watanabe  et  al..  ^75^  They  found  that  EL2a  (Ec-0.83eV) 
was  the  dominant  trap  level  for  GaAs  grown  at  720  to  740°C,  while, 
EL2b  (E  -0.76eV)  was  the  dominant  trap  level  for  GaAs  grown  at  630  to 
660°C.  In  their  annealing  studies,  Day  et  al.  ^94^  reported  that  EL2b 
was  observed  in  as-grown  MBE  n-GaAs  material;  this  level  can  be 
annealed  out  at  800°C  for  1/2  hour  or  at  700°C  for  1  hour.  Whereas, 
EL2a  was  not  observed  in  the  as-grown  MBE  n-GaAs.  The  EL2a  level  can 
be  created  by  thermal  annealing  process  and  its  trap  concentration  can 
be  enhanced  by  high  temperature  annealing.  These  results  indicate 
that  EL2  electron  trap  may  be  due  to  two  different  trap  levels.  One 
of  them  is  the  EL2a  (Ec-0.83eV)  electron  trap  and  the  other  is  EL2b 
(E  -0.76eV)  electron  trap.  Our  model  and  experimental  evidence  have 
showm  that  these  two  electron  traps  have  different  physical  origins. 

Studies  of  the  LEC  grown  GaAs  t79'95!  have  shown  that  n-type 
(S.I.)  material  can  be  grown  only  from  melts  above  a  critical  As 
composition,  and  the  EL2a  level  was  native  defect  observed  in  this 
material.  Ga-rich  melts  were  found  to  yield  p-type,  low  resistivity 
GaAs  crystal.  Ta  et  al.  ^79^  reported  that  EL2a  level  was  observed  in 
As-rich  GaAs  material.  GaAs  samples  prepared  by  LPE  method  are  grown 
from  a  Ga-rich  solution,  whereas,  VPE  GaAs  samples  are  commonly 
prepared  in  an  As-rich  gas  ambient,  in  which  growth  rate  and  surface 
morphology  are  optimized.  ^95^  The  EL2a  level  was  not  observed  in  as- 
grown  MBE  n-GaAs  material.  ^94^  However,  it  was  found  that  EL2a  level 
can  be  created  and  its  trap  density  can  be  enhanced  by  high 
temperature  annealing  (above  500°  C) .   ^94^     There  is  clear  evidence 
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that  the  EL2a  trap  can  only  be  observed  in  GaAs  grown  under  As-rich  or 
high  arsenic  pressure  case.  The  high  arsenic  pressure  will  enhance 
the  formation  of  VGa,  and  As-rich  case  will  increase  the  number  of 
As^.  ^89^  Thus,  the  EL2a  level  is  associated  with  gallium  vacancy, 
arsenic  interstitial  or  their  complexes,  ^96^  and  is  not  related  to 
oxygen  impurity  complex.    ^    '* 

In  undoped  GaAs  material,  Watanabe  et  al.  t75'89^  have  shown  that 
the  concentration  of  EL2b  level  is  proportional  to  ([As^l/fTMG])1'  . 
However,  the  concentration  of  EL2a  level  is  proportional  to 
([AsH3]/[TMG])1//2.  They  also  found  that  densities  of  EL2a  and  EL2b 
level  would  increase  with  increasing  growth  temperature  for  MOCVD 
grown  GaAs  material.  Bhattacharya  et  al.  *•  '*  observed  a  linear 
dependence  of  the  EL2a  trap  density  on  [As]  /  [Ga]  ratio  in  MOCVD 
GaAs  samples.  Lagowski  et  al.  ^78^  have  found  that  the  concentration 
of  EL2a  level  was  increased  with  increasing  As  pressure  during 
Bridgmann  bulk  growth,  while,  Miller  et  al.  ^86^  found  that  the 
density  of  EL2a  level  was  increased  with  increasing  [AsH^] /[GaCl] 
ratio  in  VPE  grown  GaAs.  Li  et  al.  ^21^  also  found  that  the  density 
of  EL2a  level  was  decreased  with  decreasing  [AsC^l/tGa]  ratio  in  the 
Ga-rich  VPE  GaAs. 

In  the  S-doped  GaAs  material,  Watanabe  et  al.  t75'89^  found  that 
density  of  EL2b  level  decreases  with  increasing  dopant  concentration 
of  surfur  impurity.  The  density  of  EL2a  level  was  found  to  decrease 
with  increasing  concentration  of  shallow  donor  dopants  (Si,  Se,  Te)  as 
was  observed  by  Lagowski  et  al..  *■'  *  Donor  concentration  above  a 
threshold  value  (1  x  1017  cm-3)  led  to  a  rapid  elimination  of  EL2 
trap,     t78-'      This   is  consistent   with  our  observation    '     *    in  which  no 


87 


EL2  trap  level  was  detected  in  MOCVD  GaAs  samples  with  Sri  dopant 
density  higher  than  3  x  10x'  cm-  .  If  group  VI  (Se,  Te,  S)  elements 
occupy  the  As  lattice  site  and  combined  with  gallium  vacancy  to  form 
anti site-vacancy  (AAsVGa~)  complexes,  then  the  concentration  of  V"Ga 
will  be  decreased.  Thus,  the  probability  of  forming  an  arsenic 
anti  site  (AsGa)  defect  will  be  greatly  reduced.  '•8'J  Johnson  et  al. 
*•  '  used  photoluminescene  upconversion  method  to  observe  EL2a  level, 
and  concluded  that  this  level  is  due  to  arsenic  anti site  defect.  The 
fact  that  EL2a  level  is  a  donor  type  defect  was  also  supported  by  the 
observed    field   dependence    of   emission    rates.  L 4° r / 3 J  This    is 

consistent  with  the  modelling  of  native  point  defect  described  in  the 
previous  chapter.  We  conclude  that  arsenic  anti  site  defect  is  the 
only  grown-in  defect  which  is  observed  in  As-rich  or  high  arsenic 
pressure  case,  but  is  not  observed  in  Ga-rich  and  low  arsenic  pressure 
case.     From  the  experimental  evidence,    it  can  be  shown  that  EL2a  trap 

is  due   to  anti site  defect,    AsG        ,    formed  during    the   post-grown 

T94  991 
cooling,     L     '     J    as  will  be  discussed  further   in  next  section. 

6.2.       Theoretical   Modelling  of  the  EL2   Electron  Trap  in  GaAs 

We  shall  next  present  a  new  model  for  explaining  EL2  level  vs 
different  [As]/[Ga]  ratio  for  the  MOCVD  and  VPE  grown  GaAs  samples. 
Assuming  that  the  mole  fraction  ratio  of  [As]  to  [Ga]  is  equal  to  r, 
it  can  be  shown  that  for  r  >  1,  the  concentration  of  EL2a  level  is 
proportional  to  (r-1)1'  ,  and  the  concentration  of  EL2b  level  is 
proportional  to  (r-1)1'  .  In  this  section,  the  kinetics  of  EL2 
formation  in  GaAs  are  described.  In  the  growth  process,  vacancies, 
intersti tials  and  anti sites  are  formed  in  high  temperature  thermal 


equilibrium.  In  the  cooling  process,  vacancies  migrate  to  form  EL2a 
and  EL2b  electron  trap.  In  the  annealing  process,  EL2b  trap  level  is 
transferred  to  EL2a  trap  level  after  thermal  annealing.  These  three 
different  processes  of  defect  formation  in  GaAs  are  discussed  next. 

6.2.1.     Growth  process 

Defect  formation  under  thermal  equilibrium  for  the  case  of  r  >  1. 
(a)      MOCVD  grown  GaAs  epitaxial  material. 

Arsine  (AsH-,)  and  trimethygallium  (TMG)  were  used  as  sources 
for  As  and  Ga  in  MOCVD  grown  GaAs.  Let  the  mole  fraction  ratio  of 
arsine  to  TMG     be  r,  then 

[AsH3]    /   [Ga(CH3)3]    =  r  (6.1) 

If  we  assume  [Ga(CH3)3]  equal  to  1,  then  [AsH3]  will  be  equal  to 
r,  where  the  square  bracket  in  Eq. (6.1)  represents  the  mole  fraction 
of  arsine  and  TMG  gas.  The  reaction  of  arsine  is  to  decompose  into 
As4 (g)  or  As2(g),  depending  on  the  growth  temperature.  In  general, 
arsine  will  decompose  into  As4(g)  in  the  epitaxial  growth  temperature 
(below  1000  °C),  while  arsine  will  decompose  into  As2(g)  at  the  melt 
growth  temperature  (above  1400  °C).  ^64^  The  growth  temperature  for 
MOCVD  process  is  usually  below  1000  °C,  and  thus  the  decomposed 
reaction  of  arsine  can  be  written  as 

AsH3(g)    =    (l/4)As4(g)    +  3/2  H2  (g)  (6.2) 

The  chemical  reaction  of  As4(g)  and  TMG  is  shown  below 
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Ga(CH3)3  +   (r/4)As4(g)   +   (3/2)H2(g) 

=  GaAs(s)    +    [(r-l)/4]As4(g)   +  3CH4 (g)  (6.3) 

If  GaAs(g)  were  completely  deposited  on  surface  of  GaAs  substrate, 
then  the  extra  As4(g)  will  he  decomposed  into  AsAs(s)  or  ASj ,  as  shown 
in  Fig.    6.1,     by  the  following  reactions    l100J 

(l/4)As4(g)    =  AsAs(s)    +  VGa  (6.4) 

(l/4)As4(g)    =  ASj  (6.5) 

(l/2)As4(g)    =  AsAs(s)    +  As^  (6.6) 

Equations  (6.4)  and  (6.5)  may  be  obtained  from  Eqs.  (5.1)  to  (5.7). 
These  two  equations  are  the  dominant  reactions  during  crystal  growth. 
The  reaction  of  Eq.  (6.6)  can  be  neglected  compared  to  Eqs.  (6.4)  and 
(6.5).     The  law  of  mass  action  for  Eqs.   (6.4  -   6.6)    is  given  by 

kl  =    [AsAs(s)]     [VGa]   /   [As4(9H1/4  (6.7) 

k2  =    [ASj]    /    [As4(g)]1/4  (6.8) 

k3  =    [AsAs(s)]     [AsGa^    /    tAs4(g)]1/4  (6.9) 

From  Eqs.    (6.7)   -    (6.9),    yields   the  concentration  of  VGa,   As^,    and 

As^ 


3Ga* 


[vGa]  =  kl  tAs4 (g) ] 1/4  /  [ASas(s)] 

=  ^[(r-D/4]1/4  =  ^'(r-1)1/4  (6.10) 

[ASj]  =  k2'(r-l)1/4  (6.11) 

[AsGa]  =  k3'(r-l)1/2  (6.12) 
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(b)     VPE  grown  GaAs  epitaxial  material. 

Arsine  (AsH^)  and  gallium  chloride  (GaCl)  were  used  as  sources 
for  As  and  Ga  in  the  VPE  grown  GaAs.  Let  the  mole  fraction  ratio  of 
arsine  and  gallium  chloride  be  r,  then 

[AsH3]   /   [GaCl]    =  r  (6.13) 

If  we  assume  [GaCl]  equal  to  1,  then  [AsHo]  will  be  equal  to  r.  The 
reaction  of  As4(g)    and  GaCl  can  be  expressed  by 

GaCl  +    (r/4)As4(g)    +    (l/2)H2(g) 
=  GaAs(s)    +   [(r-l)1/4]   As4 (g)    +  HC1 (g)  (6,14) 

From  Eq.  (6.14),  it  is  noted  that  there  is  (r-1)1/4  extra  As4  (g)  in 
the  reaction  tube,  and  thus  As  atoms  will  be  deposited  on  As  sites  or 
As  interstitial  according  to  Eqs.  (6.4)  to  (6.6).  Therefore,  the 
concentrations  of  Vga,   ASj ,   and  As^  can  be  expressed  by 

[VGa]    =  k4[As4(g)]1/4  /   [AsAs(s)]    =  k4'(r-l)1/4        (6.15) 
[ASi]    =  k5'(r-l)1/4  (6.16) 

[AsGa]    =  k6'(r-l)1/2  (6.17) 

From  Eqs.  (6.10)  -  (6.12)  and  (6.15)  -  (6.17),  it  is  noted  that  for 
both  MOCVD  and  VPE  grown  GaAs,  the  density  of  V^  is  proportional 
to  (r-1)1'  ;  the  density  of  ASj  is  proportional  to  (r-1)  '  ;  and  the 
density  of  As^   is  proportional    to    (r-1)   '    . 

6.2.2.     Cooling  process 

Formation  of  arsenic  antisite  (AsG  )  and  arsenic  anti si te-plus- 
arsenic  vacancy  (Asq3VAs)   defects  occurs  when  V^  migrates  into  As^   or 
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As,  sites  under  thermal  nonequi librium  condition.  We  believe  that 
EL2a  level  is  due  to  As^,  and  EL2b  level  is  due  to  AsQaVAs. 
Figure  6.2  shows  the  cooling  process  of  GaAs  layers. 

The  As^  defect  can  be  formed  by  V^  and  As;,    which  is  written  as 

Asi   +  vGa  =  AsGa  <6'18> 

The  AsnaVAs  can  ^e   f°riIie<3  by  V*Ga  migration   to    it's   neighbor 
(As^J    site,    and  is  shown  as 

AsAs<s>    +  VGa  "  AsGaVAs  <6-19> 

The  mass  action  law  for  the  reaction  equations  given  in  Eqs.  (6.18) 
and   (6.19)   is  shown  as 

Ka  =    [AsGa]    /    [ASj]    [VGa]  (6.20) 

Kb=    fAsGaVAS]    /   ^[VGa^     tAsAS(s^}  <6-21> 

The  As^  and  AsGaV"As  defect  concentrations  are  given  respectively  by 

[AsGa]    =  Ka    [ASj]     [VGa]    =  Ka"(r-l)1/2  (6.22) 

tAsGaVAs]    =  Kb    ^1    =  ^"(r-l)1/4  (6.23) 

Figure  6.3  shows  the  densities  of  EL2a  and  EL2b  trap  as  a  function  of 
[As]  /  [Ga]  mole  fraction  ratio  as  predicted  by  Eqs.  (6.22),  (6.23), 
and  the  published  experimental  data.  l'5,87,101J  T]ie  resuits  show 
that  for  r  >  1,  the  density  of  EL2a  trap  is  proportional  to  (r-1)1'  , 
and  the  density  of  EL2b  trap  is  proportional  to  (r-1)  '  ,  as  is 
predicted  by  Eqs.    (6.22)  and   (6,23). 
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Figure  6.3  Defect  density  of  EL2a  trap  vs.  different  As/Ga  mole 

fraction  ratios  along  with  the  published  data  for  GaAs. 
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6.2.3.     Annealing  process 

EL2a  level  may  be  formed  from  EL2b  via  thermal  annealing.  Figure 
6.4  shows  the  formation  mechanisms  from  EL2b  to  EL2a.  Note  that  EL2b 
level  may  gain  sufficient  thermal  energy,  and  decomposes  into  Vq3 
according  to  the  inverse  reaction  of  Eq.  (6.19)  when  the  annealing 
temperature  is  higher  than  500°C  &sQa  *s  formed  when  the  Vq3  migrate 
to  the  ASj  site,  according  to  Eq.  (6.18).  This  is  supported  by  the 
experimental  data  of  Taniguchi    et  al.    '•'*•'    and  Day  et  al..    *■     ' 


6.3.     Determination  the  potential  well  for  EL2  electron  trap  from 
field  enhanced  emission  rate  analysis 


As  discussed  above,  the  most  probable  physical  origin  for  the 
EL2a  trap  is  due  to  arsenic  antisite  defect,  and  the  possible  origin 
of  EL2b  trap  is  attributed  to  arsenic  anti  si  te-plus-arsenic  vacancy 
complex  defect.  The  physical  origins  of  these  electron  traps  may  be 
verified  if  the  type  of  potential  well  is  known.  This  can  be  done  by 
studying  the  field  dependent  emission  rates  deduced  from  the 
nonexponenti al  DLTS  response  for  these  electron  traps,  as  is  presented 
in  chapter    IV. 

The  type  of  potential  well  for  the  EL2a  electron  trap  was 
analyzed  from  the  theoretical  calculations  of  nonexponenti al  DLTS 
response  in  GaAs.  A  comparison  of  the  calculated  nonexponenti al  DLTS 
response  with  the  measured  DLTS  data  enable  us  to  determine  the 
potential  well  for  this  electron  trap.  We  have  calculated  the 
nonexponenti al  DLTS  response  for  four  different  types  of  potential 
wells,    namely,    Coulombic   well,    Dirac   well,    square   well    and 
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polarization  well  for  the  EL2a  trap,  and  the  results  are  shown  in 
Fig.  4.10.  Figure  4.11  shows  the  calculated  DLTS  response  comparing 
with  the  experimental  DLTS  data  for  the  Coulombic  potential  well  with 
single  and  double  charge  states.  The  capture  cross  section  (a  ) 
equals  8xl0~14  cm2  was  used  in  the  calculations.  A  comparison  of  the 
calculated  nonexponenti al  DLTS  response  with  the  measured  DLTS  data 
shows  that  the  observed  EL2a  trap  is  best  fitted  to  a  Coulombic 
potential  well  with  double  charge  state,  as  shown  in  Fig.  6.5.  Figure 
4.12  shows  the  theoretical  nonexponenti al  DLTS  response  of  Coulombic 
potential  well  with  double  charge  state  for  different  emission 
constants   [1/T  =  34.4,  86.6,  172,  344,  and  866  s-1]. 

The  field  dependent  emission  rates  for  the  Coulombic  potential 
well  can  be  derived  from  the  three-dimensional  Poole-Frenkel  effect 
and  phonon-assisted  tunneling  effect,   and  can  be  expressed  by 


2tt     tt/2  2tt     tt/2 

enHC/en0={[j        j        sin(e)exp(AEt;iAT)ded(})  +  J        I        sin(e) 
0       0  "  0       0 

(Etj-  AEfcj)AT 
/4  tt  }+J  exp{z-zJ/^  [4  (2m) l/l  (kT)  i/2/3qtfF] 


[1-    (  AEtj/zkT)5/3]}  dz  (6.24) 

where  en0  is  the  emission  rate  at  zero  electric  field,  A  Etj  is  the 
Poole-Frenkel  barrier  due  to  the  external  electric  field,  a  Et-  = 
q[qFcos(G)  /  (  tt  eQ  e  s)  ] 1'2,  and  F  is  the  applied  electric  field.  The 
first  term  in  Eq.  (6.24)  is  due  to  the  three-dimensional  Poole-Frenkel 
effect,  and  the  second  term  is  due  to  the  phonon-assisted  tunneling 
effect.  Since  the  electric  field  varies  with  position  within  the 
depletion  region  of  a   reverse  biased  p-n  junction,    the  emission   rate 
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Figure  6.5  DLTS  response  of  EL2a  trap  in  GaAs.  Solid  line: 
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is  not  constant  within  the  same  region.  The  DLTS  response  due  to  the 
nonuni formed  emission  rates  can  in  general  be  expressed  by  S(t  )  = 
I  exp(-enit1)  -  E  exp(-enit2). 

A  comparison  of  our  theoretical  calculations  of  nonexponential 
DLTS  response  due  to  the  field  dependent  emission  rates  with  the 
measured  DLTS  data  reveals  that  the  potential  well  for  the  EL2a 
electron  trap  is  due  to  Coulombic  well  with  double  charge  state,  and 
hence  the  most  likely  physical  origin  for  this  electron  trap  is 
attributed  to  Asga  .  This  conclusion  is  consistent  with  the  EL2 
modelling  based  on  the  kinetic  equations  and  defect  structure  as 
reported  in  last  section. 

6.4.  Summary  and  Conclusions 

In  this  chapter,  we  presented  a  new  modelling  for  the  EL2 
electron  trap  in  undoped  GaAs.  The  EL2  level  is  due  to  two  electron 
trap.  One  of  them  is  EL2b  (Ec-0.76eV)  electron  trap,  the  other  is 
EL2a  (Ec-0.83eV)  electron  trap.  The  physical  origin  of  EL2a  is 
attributed  to  arsenic  antisite  (AsGa)  defect;  whereas,  the  physical 
origin  of  EL2b  level  is  attributed  to  arsenic  antisi te-plus-arsenic 
vacancy  (AscaVAS)«  Base(3  on  this  model,  densities  of  EL2a  and  EL2b 
trap  levels  vs  [As]  /  [Ga]  ratio  in  the  MOCVD  and  VPE  grown  GaAs  was 
derived.  The  result  shows  that  density  of  EL2a  trap  level  is  directly 
proportional  to  the  mole  fraction  ratio  of  (r-1)1'2;  but  the  density 
of  EL2b  trap  level  is  proportional  to  the  mole  fraction  of  (r-1)1'  . 
This  prediction  is  supported  by  experimental  data  for  the  MOCVD  and 
VPE  grown  GaAs. 


VII.    STUDY  OF  GROWN-IN  DEEP  LEVEL  DEFECTS  VS  GROWTH  PARAMETERS 
IN  VPE,   LEC,   LPE  and  MOCVD  Grown  GaAs 


Studies  of  grown-in  deep  level  defects  in  GaAs  grown  by  VPE,  LEC, 
LPE,  and  MOCVD  techniques  under  various  growth  conditions  will  be 
described   in  this  chapter. 


7.1.     Study  of  Grown-in  Deep  Level  Defects  vs  Growth  Parameters 
in  the  VPE  GaAs  Layers 


Detailed  study  of  the  grown-in  deep  level  defects  vs  substrate 
orientation  and  gas  phase  stoichiometry  in  the  VPE  GaAs  layers  grown 
by  a  novel  Ga/AsClo/F^  reactor  [which  is  under  high  arsenic  pressure 
and  Ga-rich  conditions]  has  been  made,  using  DLTS  and  C-V  methods. 
Density  of  electron  traps  vs.  [Ga]/[As]  ratios  of  2/1,  3/1,  4/1,  5/1, 
and  6/1  was  determined  for  the  VPE  GaAs  layers  grown  on  <100>,  <211A>, 
and  <211B>  oriented  semi-insulating  (S.  I.)  Cr-doped  GaAs  substrates. 
Two  electron  traps  with  energies  of  Ec-0.71  eV  (i.e.,  EB4)  and  Ec-0.83 
eV  (i.e.,  EL2a)  were  observed  in  the  samples  studied.  Results  showed 
that  density  of  both  electron  traps  depends  strongly  on  the  [Ga]/[As] 
ratio  for  the  <21lA>  and  <211B>  oriented  samples,  and  less  strongly 
for  the  <100>  oriented  samples.  For  example,  in  the  <21lA>  (i.e.,  Ga- 
rich  face)  oriented  samples;  the  density  of  EB4  and  EL2a  electron 
traps  was  found  to  decrease  with  increasing  [Ga]/[As]  ratio;  while 
reverse  trend  was  found  for  the  <21lB>  (i.e.,  As-rich  face)  oriented 
samples.  As  for  the  <100>  oriented  samples,  the  density  of  EL2a  level 
may  either  increase  or  decrease  with  increasing  [Ga]/[As]  ratio.  The 
physical  origins  of  this   trap   is   attributed   the   arsenic   antisite 
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(AsG  ++) ,  as  was  discussed  in  chapter  VI.  As  for  the  EB4  level,  the 
results  are  less  conclusive,  and  further  study  is  needed. 


7.1.1.  Introduction 

Although  the  AsClo  epitaxial  growth  system  is  the  most  widely 
used  method  for  the  growth  of  microwave  device  materials,  only  limited 
information  is  available  on  how  variation  of  the  [Ga]/[As]  ratio  in 
the  gas  phase  composition  can  affect  the  materials  produced.  In  most 
AsCl^  reactors,  this  ratio  is  not  easily  controlled  being  a  function 
of  source  efficiency,  wall  deposit  extensi veness,  flow  rates,  etc.. 
If  stoichiometry  related  defects  or  defect  complexes  are  important  in 
the  VPE  GaAs  growth,  it  might  be  expected  that  the  [Gal/fAsC^]  ratio 
in  the  growth  environment  is  important.  Miller  et  al.  l  J  have 
studied  the  effect  of  gas-phase  stoichiometry  on  deep  levels  in  the 
VPE  GaAs  grown  by  the  hydride  system  with  [AsH3]/[GaCl]  ratio  varying 
from  1/3  to  3/1.  A  deep-level  electron  trap  (EL2a)  with  energy  of  Ec- 
0.83  eV  was  observed  in  their  study;  its  density  was  found  to  depend 
on  the  mole  fraction  ratio  of  Ga  to  As  atom..  Similar  studies  were 
also  reported  by  Merenda  t102J  in  the  VPE  GaAs  grown  by  the  chloride 
system  and  the  MOCVD  grown  GaAs.  i5'90'9-1-]  None  of  these  studies  have 
looked  into  the  effects  of  varying  substrate  orientation  and  [Ga]/[As] 
ratio  on  the  deep- level  defects  in  the  VPE  GaAs  epi layers  as  will  be 
reported  in  this  chapter.  To  facilitate  our  study,  a  novel  AsCl^ 
vapor  phase  epitaxial  reactor  l31'103J  was  used  to  grow  the  GaAs 
epitaxial  layers;  the  system  is  capable  of  operation  with  a  variable 
controlled    [Ga]/[As]    ratio.     The  GaAs     epi layers  were     grown  on     three 
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different  orientation  [i.e.,  <100>,  <211A>,  and  <21lB>]  substrates. 
DLTS  and  C-V  methods  were  employed  to  determine  the  density  of  deep- 
level  defects  vs.    [Ga]/[As]    ratio   in  the  VPE  GaAs  layer. 

7.1.2.     Experimental  Details 

The  VPE  GaAs  epitaxial  layers  used  in  this  study  were  grown  by  a 
novel  AsClo  reactor  system.  [31,103]  The  basjc  design  of  this  reactor 
is  like  the  usual  two  bubbler  AsCl3  reactor  used  for  the  MESFET 
growth,  [l^4!  The  novel  feature  of  our  reactor  is  the  addition  of  a 
third  bubbler  which  feeds  a  cracking  furnace  and  then  a  second  source 
chamber  to  allow  operation  in  a  GaCl  rich  mode.  It  uses  a  six-zone 
clamshell  furnace  which  rolls  off  sideways  allowing  an  end  to  end  gas 
flow  system.  The  third  or  control  bubbler  feeds  a  cracking  furnace 
that  is  maintained  at  900°C  when  in  use. 

The  GaAs  epi layers  were  grown  on  <100>,  <211A>,  and  <211B> 
oriented  S.  I.  Cr-doped  GaAs  substrates  with  growth  temperatures 
ranging  from  710°C  to  750°C  To  facilitate  our  study  of  the  effect  of 
orientation  on  the  grown-in  defects  in  the  VPE  GaAs,  seven  growth  runs 
were  made  for  this  study  (e.g.,  77A,  77B,  and  77BB  were  grown  in  the 
same  run,  etc.).  The  gas  phase  stoichiometry  was  controlled  by 
varying  the  [Ga]/[As]  ratios  from  2/1  to  6/1.  Al-GaAs  Schottky 
barrier  structures  were  fabricated  from  these  VPE  GaAs  epi layers  for 
the  DLTS  and  C-V  measurements.  Density  and  energy  levels  of  the 
grown-in  defects  in  these  VPE  GaAs  samples  were  determined  from  the 
C-V  and  DLTS  data.     The  results  are  presented  next. 
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7.1.3.  Results  and  Discussions 

In  this  section,  we  present  the  results  of  our  DLTS  and  C-V 
measurements  on  the  VPE  GaAs  samples  described  in  preceding  section. 
The  DLTS  scans  of  electron  traps  as  a  function  of  [Ga]/[As]  ratio  for 
different  runs  and  substrate  orientations  are  shown  in  Fig.  7.1 
through  Fig.  7.9.  Figure  7.1  shows  the  DLTS  scans  of  electron  traps 
as  a  function  of  the  [Ga]/[As]  ratio,  for  samples  with  <211A> 
orientation.  In  general,  one  to  two  electron  traps  (i.e.,  EB4  =  E  - 
0.71  eV,  and  EL2a  =  Ec-0.83  eV)  were  observed  in  these  samples.  Note 
that  sample  37HB  with  the  highest  [Ga]/[As]  ratio  (i.e.,  6/1)  has  the 
lowest  density  for  EL2a  electron  trap,  while  sample  35CB  with  the 
lowest  [Ga]/[As]  ratio  (i.e.,  2/1)  has  the  highest  density  of  EL2 
electron  trap.  Thus,  for  the  <211A>  orientated  samples  (i.e.,  Ga-rich 
face),  the  density  for  EL2a  level  decreases  with  increasing  [Ga]/[As] 
ratio;  this  same  trend  was  also  observed  in  EB4  electron  trap.  Figure 
7.2  shows  the  DLTS  scans  of  electron  traps  as  a  function  of  the 
[Ga]/[As]  ratio  for  epi layers  grown  on  the  <100>  oriented  substrates. 
Results  showed  that  the  density  of  EL2a  level  may  either  increase  or 
decrease  with  increasing  [Ga]/[As]  ratio,  depending  on  the 
concentrations  of  VGa  or  As^  formed  during  crystal  growth.  To  explain 
this,  modelling  presented  in  chapter  VI  appears  to  be  adequate  for 
this  high  arsenic  pressure  and  Ga-rich  condition.  We  suggested  that 
EL2a  level  is  attributed  to  the  As^g  arsenic-antisi te  defect  which  is 
formed  by  the  migration  of  VG  defect  to  ASj  site  after  crystal 
cooling,  as  shown  in  Eq.  (6.18).  VPE  GaAs  layers  are  grown  under  high 
arsenic   pressure   condition.      In   the  case  of   <211A>   oriented   samples, 
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Figure   7.2     DLTS  scans  of  electron  traps   for  VPE  GaAs   samples   grown 
on   (100)  oriented  substrates  with  Ga/As   ratio   varying 
from  2/1    to  6/1. 
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Figure  7.5  DLTS  scans  of  electron  traps  for  two  VPE  samples  with 
(100)  and  (211A)  orientations  and  Ga/As  ratio  equal  to 
4/1. 
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Figure  7.9     Density  of  EL2a   level    vs.    Ga/As   ratio   in   VPE   GaAs 
samples   grown  on    (211A),    (21  IB),   and   (100) 
oriented  substrates. 
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Vq  should  be  suppressed  by  the  following  reaction  in  the  Ga-rich 
condition. 

VGa  +  Gad  =  Ga^  (7.1) 

Thus,  the  density  of  EL2a  level  will  decrease  with  increasing 
[Ga]/[As]  ratio,  as  is  the  case  of  <211A>  samples  shown  in  Fig.  7.1. 
On  the  other  hand,  in  the  case  of  <211B>  oriented  samples,  arsenic 
face  will  prevent  the  reaction  of  Eq.(7.1)  in  the  Ga-rich  case.  At 
this  case,  the  concentrations  of  As^  and  VV^  levels  shall  be  enhanced 
with  increasing  [Ga]/[As]  ratio.  Therefore,  the  density  of  EL2a  level 
will  increase  with  increasing  [Ga]/[As]  ratio;  this  is  the  case  for 
<21lB>  oriented  samples.  This  defect  model  seems  adequate  for 
predicting  the  orientation  dependence  of  the  density  of  EL2a  vs 
[Ga]/[As]  ratio,  as  is  shown  in  Fig.  7.9.  Figure  7.3  shows  the  DLTS 
scans  of  electron  traps  for  sample  35CA  <100>  and  sample  35CB  <211A> 
with  [Ga]/[As]  ratio  equal  to  2/1.  These  two  samples  were  grown  in 
the  same  run  side  by  side  at  715°C.  The  DLTS  results  showed  that 
density  of  EL2a  electron  trap  for  the  <211A>  sample  is  significantly 
higher  than  that  of  the  <100>  sample  (i.e.,  1.24xl014  cm-3  vs.  7.6  x 
10X^  cm  J).  The  reason  for  this  may  be  attributed  to  the  fact  that 
during  crystal  growth  in  high  arsenic  pressure  condition,  the  density 
of  Vq3  is  higher  in  the  <211A>  face  than  that  in  the  <100>  face  for 
[Ga]/[As]  ratio  equal  to  2/1.  The  EB4  electron  trap  was  not  observed 
in  both  samples  shown  in  Fig.  7.3.  Figure  7.4  shows  the  DLTS  scans  of 
electron  trap  for  sample  75A  <100>  and  sample  75B  <211A>  with 
[Ga]/[As]  ratio  equal  to  3/1.  Note  that  the  EL2a  level  was  observed 
only  in  sample  75A  while  EB4  level  was  observed  only  in  75B  sample. 
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The  growth  temperature  for  both  samples  was  held  at  725°C.  Figure  7.5 
shows  the  DLTS  scans  of  electron  traps  for  samples  76A  <100>  and  76B 
<211A>.  The  results  showed  that  both  EL2a  and  EB4  electron  traps  were 
observed  in  these  two  samples.  Density  of  EL2a  and  EB4  levels  was 
higher  for  sample  76B  <21lA>  than  sample  76A  <100>.  Figure  7.6  shows 
the  DLTS  scans  of  electron  traps  for  samples  77B  <211A>,  77A  <100>, 
and  77BB  <211B>.  These  three  samples  were  grown  side  by  side  at  750°C 
and  with  [Ga]/[As]  ratio  equal  to  4/1.  Both  EL2a  and  EB4  electron 
traps  were  observed  in  these  three  samples;  the  density  of  EL2a  and 
EB4  levels  was  highest  for  sample  77BB,  followed  by  samples  77B,  and 
77 A.  Figure  7.7  show  the  DLTS  scans  of  electron  traps  for  samples  78B 
<211A>  and  78BB  <211B>  grown  at  725°C  and  with  [Ga]/[As]  ratio  equal 
to  5/1.  The  results  showed  that  EL2a  level  was  the  dominant  electron 
trap  for  sample  77BB  while  only  EB4  level  was  observed  in  sample  78B. 
Figure  7.8  shows  the  DLTS  scans  of  electron  traps  for  samples  37HA 
<100>  and  37HB  <211A>  grown  at  725°C  and  with  [Ga]/[As]  ratio  equal  to 
6/1.  The  results  showed  that  both  EL2a  and  EB4  level  was  observed  in 
sample  37HB.  To  sum  up  the  results  shown  in  Fig.  7.1  through  Fig. 7. 8 
concerning  the  effects  of  orientation,  [Ga]/[As]  ratio,  and  the  growth 
temperature  on  the  density  of  EL2a  and  EB4  levels,  we  plot  in  Fig. 
7.9,  the  density  of  EL2a  level  vs  [Ga]/[As]  ratio  for  three  substrate 
orientations.  In  Fig.  7.9,  it  is  clearly  shown  that  density  of  EL2a 
level  was  found  to  decrease  with  increasing  [Ga]/[As]  ratio  for  the 
<211A>  oriented  samples  and  to  increase  with  increasing  [Ga]/[As] 
ratio  for  the  <211B>  oriented  samples;  for  the  <100>  oriented  samples 
the  density  of  EL2a  level  may  increase  or  decrease  with  increasing 
[Ga]/[As]    ratio,  depending  on  the  V"Ga  or  As^  defects  formed  during 
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crystal  growth.  Table  7.1  summarizes  the  defect  parameters  deduced 
from  our  DLTS  and  C-V  measurements  on  twelve  VPE  GaAs  samples  prepared 
under  different  growth  conditions  and  substrate  orientations. 

7.1.4     Summary  and  conclusions 

A  study  of  the  grown- in  deep  level  defects  in  VPE  GaAs  vs. 
substrate  orientations,  [Ga]/[As]  ratios,  and  growth  temperature  has 
been  carried  out  in  this  chapter.  One  to  two  electron  traps  (i.e., 
EB4  =  Ec-0.71  eV,  and  EL2a  =  Ec~0.83  eV)  have  been  observed  in  these 
VPE  GaAs  samples.  From  the  results  of  our  DLTS  and  C-V  analysis,  it 
is  found  that 
(1) .     background   carrier   concentrations   are   two   to   five  times  higher 

for  the  <211A>  samples  than  the  <100>  samples; 
(2).     defect  density     was    in  general  higher   for   the  <211>  samples than 

the  <100>  oriented-samples; 
(3).  density  ofEL2a   level   was      found   to  decrease   with   increasing 

[Ga]/[As]    rati  of or    the   <211A>   samples,    and    increases    with 

[Ga]/[As]    ratio  for  the  <211B>  samples; 
(4).     for  the  <100>  samples,    the  density  of  EL2a  level  may  increase  or 

decrease  with     [Ga]/[As]     ratio    depending    on    the   As^    or    VGa 

defect  density  formed  during  crystal  growth; 
(5).   the   EL2a level    is   attributed   to   the  arsenic  antisite  defect 

(As^g-^) ,   formed  after  crystal  cooling; 
(6).     for   the  EB4  electron   trap, which   is  not  commonly  observed   in   the 

VPE  GaAs, its  origin  is  not  known  and  further   study  is  needed   in 

order   to  gain  a  better  understanding  on  this  defect  level. 
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7.2.     Study  of  Grown-in  Deep  Level  Defects  vs  Growth  Parameters 
in  the  LEC  GaAs  Layers 


In  this  section,  we  present  the  results  of  our  study  of  the 
native  defects  observed  in  the  LEC  grown  bulk  GaAs.  Both  p  type  GaAs 
and  n  type  GaAs  have  been  studies.  An  acceptor  level  located  at  80 
meV  from  the  valence  band  edge  has  been  observed  in  the  LEC  p-GaAs 
materials.  For  n-GaAs,  a  study  of  the  effect  of  heat  treatment  on  the 
native  defects  on  samples  annealed  at  200,  300,  and  500°C  for  one  hour 
in  hydrogen   (H2)   ambient  has  also  been  carried  out  in  this  work. 

7.2.1     Intrinsic  double  acceptor  level  in  LEC  grown  p-GaAs 

Two  samples  of  p-type  conducting  GaAs  crystals  grown  by  the  LEC 
method  were  chosen  for  this  study.  Ohmic  contacts  were  made  by 
evaporation  of  Ag-Mn  alloy,  and  an  Al-Schottky  barrier  structure  was 
used  for  C-V  and  DLTS  measurements. 

Figure  7.10  shows  the  DLTS  scan  of  hole  trap  observed  in  both 
sample  A  and  B.  The  apparent  activation  energy  was  obtained  from  the 
Arrehenius  plot  of  the  hole  emission  rate  e  /T  vs.  1/kT.  The 
measured  activation  energy  for  this  hole  trap  is  130  meV  from  the 
valence  band.  The  true  activation  energy  of  this  hole  trap  is  80  meV 
if  correction  due  to  the  temperature  dependence  capture  cross  section 
is  taken  into  account  [i.e.,  a  p  =  7.1x10  X3  exp(- AE^/kT)  cm  ,  where 
A  EL  =  50  meV].  Therefore,  it  is  evident  that  this  hole  trap  observed 
by  the  DLTS  experiment  is  the  same  center  observed  in  the  1.441  eV 
emission  under  photoluminescence  experimental  excitation  and  Hall 
measurements.   ^'^       This     hole     trap     is     a     double  charged     acceptor 
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center,    and    is   due   to    the   gallium   antisite   defect,    Ga,    ~  .      The 
density  of  this  hole  trap  was  found  to  vary  from  lid  ^  to  10-*-"  cm     . 


7.2.2     Study  of  deep  level  defects  vs  annealing  temperature  in  Hq 
ambient 


In  the  following,  the  results  of  the  TSCAP  and  DLTS  measurements 
on  the  LEC  n-GaAs  materials  annealed  in  hydrogen  ambient  are 
discussed.  Figure  7.11  shows  the  TSCAP  scans  in  the  LEC  grown  GaAs 
Schottky  diodes  at  a  reverse  bias  equals  to  -2  volts.  From  this 
figure,  it  is  noted  that  the  capacitance  values  for  samples  annealed 
at  200  and  300°C  for  one  hour  in  H?  ambient  show  little  change  from 
the  unannealed  sample.  However,  as  the  annealing  temperature  raises 
to  500°C,  the  capacitance  value  decreases  by  more  than  50%  from  the 
unannealed  sample,  indicating  a  large  reduction  in  free  carrier 
density  in  the  500°C  annealed  sample.  This  may  be  due  to  the 
introduction  of  the  EL2a  electron  trap  by  the  500°C  annealing  process, 
as  is  evidenced  by  our  DLTS  data.  The  capacitance  steps  observed  in 
the  TSCAP  scans  such  as  those  occurred  at  130  and  240°K  for  sample 
annealed  at  500°C  are  due  to  the  electron  emission  from  the 
corresponding  electron  traps,  the  exact  energy  level  for  these 
electron  traps  can  be  determined  from  the  DLTS  technique. 

The  DLTS  scans  of  electron  traps  for  the  LEC  grown  GaAs 
unannealed  and  annealed  at  500°C  in  hU  ambient  for  one  hour  are  shown 
in  Figs.  7.12  and  7.13.  The  defect  density  deduced  from  the  DLTS 
measurements  are  summarized  in  table  7.2.  For  samples  unannealed  or 
annealed  at  200  and  300°C  for  one  hour,  there  are  three  electron  traps 
with  energies  of     Ec-0.35eV,  Ec-0.61eV,  and  EL2b  =  Ec-0.76eV  oberved 
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Table  7.2  Defect  parameters  for  the  LEC  grown  GaAs  for 
different  annealing  temperatures. 


Sample 

ND(cm~3)     ET(eV)       NT(cm"3)        a  n(cm2) 

Unanneal 

7.1xl015     Ec-0.35      3.3x10^       6.0xl0~14 
Ec-0.61      1.3x10*3       1.2xl0~j3 
Ec-0.76      1.6xl015       7.1xl0-14 

200°C 
Annealing 

5.0X1015     E  -0.35      5.0xl014       6.0xl0~14 
E  -0.61      4.0xl0|4       1.2xl0"j3 
Ec-0.76      2.5xl014       7.1xl0"14 

300°C 
Annealing 

4.3xl015     Ec-0.35      1.1x10^       6.0xl0~14 
Ec-0.61      1.0xl014       1.2x10"^ 
Ec-0.76      1.3xl014       7.1xl0"14 

500°C 
Annealing 

3.3xl015     Ec-0.35      6.0X10?-4       6.0xl0~?-4 
Ec-0.61      2.6xl0J3       1.2xl0~|3 
Ec-0.83      7.2xl014       5.0xl0"13 
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in  these  samples.  As  the  annealing  temperature  increases  to  500°C,  a 
new  electron  trap  with  energy  of  Ec-0.83eV  (EL2a  level)  is  emerged 
into  the  DLTS  scans  along  with  two  other  electron  traps  (i.e.,  Ec-0.35 
and  0.61eV),  while  EL2b  level  disappeared;  this  is  shown  in  Fig.  7.13. 
The  Ec-0.35eV  level  is  the  main  electron  trap  observed  in  these  LEC 
grown  GaAs  samples.  This  trap  has  also  been  reported  by  Martin  et  al. 
t77l  in  the  bulk  GaAs  and  by  Fairman  et  al.  t106J  in  the  VPE  grown 
GaAs.  The  reduction  in  free  carrier  density  with  increasing  annealing 
temperature  along  with  large  capture  cross  section  for  these  electron 
traps  indicating  that  they  are  acceptor  type  electron  traps.  Another 
interesting  study  on  these  traps  is  concerned  with  the  spatial 
dependence  of  the  defect  density  for  each  electron  trap.  This  is 
obtained  by  performing  the  bias  dependence  of  the  DLTS  scans,  and  the 
results  are  illustrated  in  Fig.  7.14.  From  the  DLTS  scans  as  a 
function  of  the  reverse  bias  voltage,  the  spatial  dependence  of  trap 
density  is  determined  for  each  trap.  Figure  7.15  show  the  spatial 
dependence  of  the  electron  traps  observed  in  the  200°C  annealed 
samples.  It  is  shown  that  for  all  the  samples  studied  here,  the 
defect  density  is  highest  at  the  surface,  and  decreases  as  it  moves 
away  from  the  interface  of  the  metal  GaAs  Schottky  contact;  the  trap 
density  becomes  constant  deep  into  the  bulk  region  (i.e.,  2um)  of  the 
GaAs  substrate,  indicating  that  all  these  electron  traps  are  bulk 
related  defects.  Mircea  et  al.  ^30^  have  also  observed  a  similar 
defect  density  profile  for  the  Ec-0.61eV  level  in  which  the  density  of 
Ec-0.61eV  was  increased  by  the  mechanical  damage  near  the  surface  of 
GaAs  specimen. 

Figure  7.16  shows  the  spatial  dependence  of  the  Ec-0.35eV 
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Figure  7.15  Density  of  electron  traps  in  the  200  C  annealed  GaAs 
sample  as  a  function  of  depletion  layer  width. 


127 


<S) 

CD 

c 

cr. 

CL 

S_ 

f- 

+J 

(13 

l/l 

C 

O 

■n 

J. 

0) 

+-> 

i — 

U 

nr 

q; 

aj 

c 

cd 

c 

>,  a; 

•r-    <d 

C  4- 
c:  T3 


128 


electron  trap  in  the  unannealed  and  annealed  samples.  The  results 
showed  that  density  of  the  Ec-0.35eV  trap  was  decreased  with 
increasing  annealing  temperature.  The  EL2b  trap  was  annealed  out 
followed  by  a  500°C  annealing  in  H2  ambient  for  one  hour.  Recently, 
Day  et  al.  ^94^  has  observed  the  Ec-0.35eV  and  EL2b  (Ec-0.76eV)  traps 
in  the  MBE  grown  n-GaAs,  and  found  that  both  Ec-0.35eV  and  EL2b  traps 
can  be  annealed  out  at  800°C  for  1/2  hour.  To  explain  their  results, 
the  EL2  model  proposed  in  previous  chapter  appears  to  be  adequate  for 
this  case.  The  physical  origin  of  EL2b  trap  is  attributed  to  AsGaVAs 
complex.  At  500°C  annealing  temperature,  the  EL2b  level  will 
decompose  into  V"Ga,  and  its  reaction  equation  is  given  by 

AsGaVAs  =  AsAs  +  VGa  (7'2) 

Note  that  at  high  annealing  temperatures  (e.g.,  Ta  >  500°C) ,  the  trap 
may  be  formed  by  the  migration  of  VGa  defect  to  the  ASj  site  to  form 
an  arsenic-antisi  te  defect  as  was  predicted  by  Eq.  (6.18).  This 
defect  model  seems  adequade  for  explaining  change  of  the  EL2b  trap  to 
EL2a  level  at  high  annealing  temperatures  as  shown  above. 


7.3.     Study  of  Grown-in  Deep  Level  Defects  vs  Growth  Parameters 
in  the  LPE  n-GaAs  Layers 


Samples  of  LPE  GaAs  were  prepared  for  studying  grown-in  defects 
vs  growth  temperature  and  growth  rate  in  the  high  purity  n-GaAs 
epi layers.  These  samples  were  prepared  by  infinite  solution  melt  LPE 
technique  grown  at  700  and  800°C  with  cooling  rate  of  0.4  and 
l°C/min.,  respectively.  Figure  7.17  shows  the  DLTS  scans  of  electron 
trap    (EB4=E„-0.71eV)    for    samples   grown   at   700   and   800°C  with   cooling 
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rate  of  l°C/min..  Note  that  for  the  same  cooling  rate,  sample  grown 
at  800°C  has  defect  density  six  times  lower  than  that  of  the  700°C 
grown  sample.  On  the  other  hand,  for  sample  grown  at  700  and  800°C 
with  a  cooling  rate  of  0.4°C/min.,  no  electron  traps  were  detected  in 
these  samples.     The  results  of  our  DLTS  and  C-V  measurements  on  the 

LPE  GaAs   were   summarized   in   table  7.3.      The  background   concentration 

1  fi         3 
in  these  samples  is  around  10       cm     . 


7.4.     Study  of  Grown-in  Deep  Level  Defects  vs  Growth  Parameters 
in  the  MOCVD  n-GaAs  Layers 


Characterization  of  native  defects  in  n-GaAs  epi layers  grown  on 
S.  I.  GaAs  and  S.  I.  Ge  substrates  by  the  MOCVD  technique  has  been 
studied.  The  Sn-doped  GaAs  epi  layer  on  S.  I.  GaAs  substrate  was  grown 
at  700°C  for  30  min.  with  [Ga]/[As]  ratio  of  11/1,  and  the  dopant 
density  was  4.1xl017  cm-3.  As  for  the  Sn-doped  GaAs  epi  layer  on  the 
S.  I.  Ge  substrate,  the  sample  was  prepared  at  706°C  for  30  min.  with 
[Ga]/[As]  ratio  of  7/1;  the  background  concentration  was  found  to  vary 
between  1016  to  1.9xl017  cm-3.  Figure  7.18  shows  the  DLTS  scan  for  n- 
GaAs  layers  on  S.  I.  GaAs  substrate.  The  result  shows  that  there  is 
one  electron  trap  with  energy  of  Ec-0.6eV.  This  trap  level  is 
corresponding  to  the  Ec-0.61eV  for  the  LEC  grown  GaAs  samples.  As  for 
the  GaAs  layer  on  S.  I.  Ge  substrate,  one  hole  trap  with  energy  of 
Ev+0.75eV  and  the  density  of  1.7xl013  cm-3  was  observed  in  this 
sample . 
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7.5.     Summary  and  Conclusions 

From  the  results  of  our  study  of  the  grown- jn  defects  in  the  LEC, 
VPE,  LPE,  and  MOCVD  GaAs  samples,  it  is  clear  that  grown-in  defects 
and  the  background  concentration  are  sensitive  to  the  growth 
conditions,  stoi chi ometry,  and  substate  orientations.  Table  7.4  shows 
the  trap  level  observed  in  this  study.  The  main  conclusions  are 
listed  as  follows 
(1).  The  EL2a   (Ec-0.83eV)   electron  trap  is  attributed  to  the  AsGa++ 

arsenic-anti si te  defect  which  was  observed   in  VPE  GaAs,    and 

500°C  annealed    (H2)    LEC  GaAs  layers. 
(2).    The     EL2b    (Ec~0.76eV)    electron    trap   is   attributed    to    the   AsGaVAs 

complex   which   was   observed   in   the  LEC  GaAs  layers.      This   level 

was  annealed  out  by  a  500°C  hydrogen  gas  heat  treatment. 
(3).   The     Ey+0.08eV  hole  trap  was  observed  in  the  LEC  grown  p-GaAs 

samples.     The  physical  origin  of  this  hole  trap  is  believed  to 

be  due  to  Ga,s  gallium-antisi te  defect  as  confirmed  by  the  Hall 

effect  and  photoluminescence  measurements. 
(4).     Reducing    the   temperature   cooling   rate  during  crystal   growth   is 

beneficial    for   reducing   the  density  of  grown-in  defects    in   LPE 

GaAs  epilayes. 
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VIII      DEFECT   STUDIES    IN  LOW-ENERGY   PROTON  AND  ONE-MEV  ELECTRON 
IRRADIATED  AlGaAs-GaAs   SOLAR  CELLS 


The  (AlGa) As-Ga As  solar  cells  with  AMO  conversion  efficiently  of 
18-19%  have  been  reported  by  Hovel  and  Woodall,  I10 'J  and  more 
recently  by  Loo  et  al..  I-10°-I  The  GaAs  solar  cell  has  shown  better 
radiation  resistance  characteristics  than  the  silicon  solar  cell,  and 
hence  it  is  potentially  a  better  candidate  for  space  applications.  It 
is  well  known  that  degradation  in  solar  cell  performance  occurs  if  it 
were  constantly  bombarded  by  high  energy  electrons  and  protons.  Thus, 
it  is  essential  to  obtain  a  better  understanding  of  the  fundamental 
limitations  on  solar  cell  performance  under  different  irradiation 
conditions.  Although  a  considerable  amount  of  work  on  radiation 
damage  in  silicon  solar  cells  has  been  reported  in  the  literature, 
L10°J  only  a  handful  of  publications  have  dealt  with  radiation  damage 
study  in  GaAs  solar  cells,    f108-112! 

8.1     Introduction 

The  purpose  of  this  chapter  is  to  report  the  results  of  our  DLTS 
analysis  of  the  low-energy  proton-irradiation  and  the  one-MeV  electron 
irradiation  induced  deep  level  defects  and  their  effects  on  the 
recombination  parameters  and  performance  in  the  liquid  phase  epitaxy 
(LPE)  grown  (AlGa) As-GaAs  solar  cells.  For  the  case  of  low-energy 
proton  irradiation,  lJ-B»J-i0J  the  GaAs  solar  cells  were  bombarded  under 
different  proton  energies  (50,  100,  200,  and  290  keV),  proton  fluences 
(1010,  1011,  5X1011,  1012,  and  1013  p/cm2)    and  Sn-dopant  densities 
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(2xl016  and  8X10-1-6  cm-3).  For  the  case  of  one-MeV  electron 
irradiation,  [20/21]  both  undoped  and  Sn-doped  GaAs  solar  cells  were 
irradiated  under  different  temperatures  (25  and  200°C) ,  fluences 
(1014,  1015,  5xl015,  and  1016  e/cm2),  fluxes  (2xl09  and  4xl010),  and 
annealing  times  (10,  20,  30,  and  60  min.).  Carrier  removal  rate  was 
derived  from  the  theoretical  bombardment  method  and  compared  with  the 
experimental  results.  The  proton  irradiation  was  carried  out  by  using 
hydrogen  implantation,  and  the  proton  profile  is  followed  by  the 
Gaussian  distribution.  L113/H4 J  The  penetration  depth  of  proton 
particles  is  proportional  to  the  proton  energy.  Theoretical  analysis 
of  the  radiation  damage  due  to  electron  and  proton  bombardment  is 
presented   in    section   8.3. 

The  DLTS  technique  was  employed  to  measure  the  defect  energy 
levels  and  densities  as  well  as  carrier  capture  cross  sections  and 
lifetimes.  The  trap  levels  produced  by  low-energy  proton  and  one-MeV 
electron  irradiation  are  found  to  be  quite  different.  The  defect 
spectrum  for  the  low-energy  proton  irradiation  is  more  complicated 
than  that  of  the  one-MeV  electron  irradiation.  Minority  carrier 
diffusion  lengths  for  the  one-MeV  electron  irradiated  cells  were 
derived  from  the  DLTS  data.  Effective  diffusion  lengths  of  the  low- 
energy  proton-irradiated  samples  was  carried  out  by  the  Wysocki  et 
al.'s  model  [l*-5]  $ue  to  n0nuni  formly  diffusion  length  of  the 
irradiated  GaAs  cells.  Defect  parameters  for  the  deep-level  traps 
induced  by  low-energy  proton  irradiation  are  presented  in  section  8.4. 
Section  8.5  discusses  defect  parameters  due  to  one-MeV  electron 
irradiation.     Summaries  and  conclusions  are  given   in  section  8.6. 
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8.2     Experimental  Details 

(AlGa)As-GaAs  mesa  diodes  used  in  this  study  were  fabricated  by 
infinite  solution  melt  LPE  techniques.  ll°»20r21J  These  test  diodes 
have  identical  device  structure  (see  Fig.  8.1)  and  were  processed 
under  identical  conditions  and  with  same  device  structure  as  the  large 
area  (2x2  cm2)  GaAs  solar  cells  reported  by  Loo  et  al..  1-108 J  The 
p  (AlGa)As  window  layer  has  a  dopant  density  of  10  cm  J  and 
thickness   of   0.3   um   or    less;    the   p+GaAs    layer    was    formed   by   Be 

IP  O 

diffusion  with  a  dopant  density  of  10  cm  J  and  thickness  of  0.5  um, 
and  the  n-GaAs  active  layer  of  10  um  thick  was  grown  subsequently  on 
the  n+GaAs  Substrates  with  2°  off  the  (100)  orientation.  The  AuZn 
ohmic  contacts  on  the  window  layer  are  about  0.3  to  0.4  um  thick.  The 
n-substrate  contact  is  AuGeNi  of  about  0.5  um  thickness  with  Ag 
overlay. 

For  proton  irradiation  study,  Sn-doped  GaAs  with  dopant  densities 
of  2xl016  and  8x10  cm-3  were  used.  The  proton  energies  used  are  50, 
100,  200,  and  290  keV,  and  proton  fluences  vary  from  1010  to  1013 
p/cm  .  The  proton  penetration  depth  is  around  0.5um  for  50  keV,  1.0um 
for   100keV,    1.8um  for   200keV,   and  2.5um  for  290  keV. 

For  one-MeV  electron  irradiation,  electrons  penetrated  through 
the  entire  GaAs  solar  cells,  and  produced  uniform  damages.  Two  types 
of  GaAs  solar  cells  were  fabricated  for  the  electron  irradiated  study. 
The  undoped  GaAs  with  carrier  density  of  1.5xl0X3  cra"J  and  the  Sn- 
doped  GaAs  with  dopant  density  of  5xl016  cm-3.  For  Sn-doped  GaAs 
cells,  electron  irradiation  was  performed  at  room  temperature  for 
fluences  of  12      ,  10      ,  and  1016  cm-2,  and  subsequently  annealed  at 


138 


< 


0 


o 

o 

O 

~ 

£ 

5 

5 

II 

II 

tl 

<0 

lO 

< 

< 

< 

o 

o 

0 

0 

0 

0 

O 

\s 

**s 

+ 

c 

4- 
C 

Q 

X 

O       *z       -*- 


2.         C  =. 


N> 


0 

In 


Q 

< 

■   0 

X 

< 

a 

<  i  <  ! 

o  •    a    1 

OiO    i 

+■ 

a  1    c    1 

139 


230  C  for  10,  20,  30,  and  60  minutes.  In  addition,  irradiation  was 
also  made  on  the  Sn-doped  GaAs  for  flux  rates  of  2x10  and  4x10 
e/cm2-s.,  for  a  total  fluence  of  1015  e/cm2  and  at  150  and  200°C 
cell's  temperature.  For  the  undoped  GaAs,  electron  irradiation  was 
performed  at  200°C  cell's  temperature  for  fluence  of  1014  and  1015 
e/cm2,  and  at  room  temperature  for  four  di  f  ferent  electron  f luences 
(i.e.,  1014,  1015,  5xl015,  and  1016  e/cm2). 

DLTS  and  C-V  measurements  technique  were  employed  to  determine 
the  defect  and  recombination  parameters  in  both  low-energy  proton  and 
one-MeV  electron   irradiated    (AlGa)As-GaAs  solar   cells. 


8.3     Theoretical  Analysis  of  Lattice  Damages 
Created  by  Particles  Irradiation 


When  high  energy  particles  (electrons  or  protons)  enter  the  GaAs 
lattice,  elastic  collisions  between  particles  and  GaAs  lattice  atoms 
are  encounted  which  results  in  lattice  damage.  Moreover,  the  energy 
imparted  to  the  atom  must  be  greater  than  the  displacement  energy  (V"D) 
for  the  damage  to  take  place.  The  exact  value  of  the  displacement 
energy  will  depend  on  colli  son  mechanisms,  lattice  structure,  and  the 
recoil  direction.  An  average  value  for  the  displacement  energy  in 
GaAs  is  25  eV.  [116,117]  Each  time  a  particle  collides  with  a  lattice 
atom,  the  maximum  energy  that  can  be  transferred  to  the  lattice  atom 
is  given  by 


[118] 


Tm  =    [4  M-l  M2  /    (M-l  +  M2) 2]    Ej  (8.1) 

where  E^   is  the  incident  particle  energy.     Tm  is  the  maximum  energy 
transferred   to   the   lattice  atom.      Mi,    and  M2  are  the  mass  of  the 
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incident  particle  and  the  lattice  atom.     For  electron  bombardment, 

relativistic  effects  are  important,  and  we  have  used  the  expression 
[118] 

Tm  =  2Ej    (Ej    +  2mQc2)    /  M2c2  (8.2) 

where  mQ  is  the  electron  mass,  and  c  is  the  velocity  of  light. 
According  to  Eq.  (8.1),  a  proton  particle  can  transfer  up  to  5%  energy 
to  the  lattice  atom,  while,  an  electron  particle  can  only  transfer  a 
small  amount  of  its  energy  to  the  lattice  atom. 

An  energetic  particle  penetrate  into  a  solids,  the  major  energy 
loss  will  be  from  collisions.  The  energy  loss  i11^  is  due  either  to 
the  nuclear  stopping  (i.e.,  elastic  collisions  between  electrons  and  a 
screened  Ga  or  As  nuclear)  or  the  electronic  stopping  (i.e.,  the 
inelastic  energy  loss  due  to  particle  passing  through  free  electrons 
and  bound  electrons).  Nuclear  stopping  is  the  dominant  process  when 
the  incident  particles  have  low  energy.  Electronic  stopping  is 
important  if  the  particles  have  high  energy  or  the  mass  ratio  of 
incident  atom  to  lattice  atom  is  small.  For  the  case  of  low  energy 
proton  irradiation,  the  mass  ratio  is  about  1/72,  both  nuclear  and 
electronic  stopping  will  be  important.  The  proton  particle  will  be 
implant  to  the  lattice  atom.  However,  for  one-MeV  electron 
irradiation,  the  mass  ratio  is  very  small,  and  hence  the  electronic 
stopping  is  the  dominant  mechanism  for  the  energy  loss. 

For  the  case  of  proton  irradiation,  the  distribution  of  protons 
at  rest  in  the  solid  is  approximately  a  Gaussian  distribution,  and  is 
given  by    [133,114] 
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N(x)    =   {<$>p/[(2*  )1/2  Rp]}exp{-0.5[(x-Rp)/A  ty2}      (8.3) 

where  $_  is  the  proton  fluence,  EL  is  the  projected  range,  and  ARp 
is  the  straggle  range.  The  Gaussian  distribution  for  protons  of 
different  energies  are  shown  in  Fig.  8.2.  From  this  figure,  it  is 
clear  that  even  at  a  proton  fluence  of  10  p/cm  ,  the  peaks  of  the 
distribution  can  be  as  high  as  10  p/cm  .  Therefore,  it  is  obvious 
that  even  at  this  proton  fluence  the  lattice  defects  can  be  easily 
created.  In  general,  the  defect  distribution  is  more  skew  '•J-20J  and 
the  defect  peak  is  shallower  '  121,122 J  than  the  proton  Gaussian 
distribution.  The  estimated  defect  profile  is  shown  in  Fig.  8.3. 
Figure  8.4  shows  the  carrier  concentration  profile  for  a  200  keV,  10 

n 

p/cm*  proton  irradiated  GaAs  solar  cell.  We  can  see  that  most  of  the 
carriers  have  been  removed  within  the  proton  range.  The  carrier 
removal  profile  is  similiar  to  the  opposite  shape  of  Gaussian 
distibution.  Near  the  end  of  the  proton  travel,  the  carrier  density 
shows  a  peak  profile.  This  may  be  attributed  to  the  ex  i  stance  of 
protons   in  this  region.    1113/123] 

For  electron  irradiation,  the  relativistic  expression  given  by 
Mott  lJ""'J"«J  for  the  differential  scattering  cross  section  may  be 
used  for  calculation.  The  total  cross  section  is  given  by  using  the 
a  approximation  given  by  Mckinley  et  al.  l126J  and  the  number  of 
displacement  proposed  by  Grimshaw  et  al.,    ^^'^   which  reads 

<5d  =A2{[Tm/(2Vd)]+2Tr  a    g{  1-  [  (1/(2) 1/2]  (Tm/Vd) 1/2-  (  B  2+ira   (3  ) 
In  (2)}     +     0.56  (Tm/Vd)    [ln(Tm/2Vd)    -    (  B2  -  ua  0  ) 
+    (2Vd/Tm)  (B2  +  tt  a  3  )  -2  tt  a  B  (2Vd/Tm) 1/2]  }  (8.4) 
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where  A2  =  (1/ir  )  [Ze2/(4e  sm0c2)]  2,  a  =  Z  /  137,  3  =  v  /c,  Z  is  the 
atomic  number  of  the  lattice  atom,  v  is  the  velocity  of  the  electron. 
The  electron  energy  loss  rate  ,  -dE^/dR  ,  deduced  from  empirical 
formula,    [128,129]    js  given  by 

dEj/dRg  =    (Ej/Rg)    [1.265  -  0.191  logCE^]-1       for  Ej   <  2.5 
and 

dEj/dRe  =  0.0189  for  E  >  2.5  (8.5) 

where  Ej  is  the  incident  electron  energy  in  MeV,  Re  is  the  proton 
range  of  electron  particle  in  mg/cm  .  From  Eq.  (8.5),  the  penetrate 
depth  of  one-MeV  electron  particles  is  about  9mm,  which  is  much 
greater  than  the  cell  thickness.  The  carrier  removal  due  to  electron 
bombardment  is  therefore  uniform  throughout  the  cell.  The  carrier 
removal  rate  is  expressed  by 

AN     /  <j>e     =  N     o-d  d  (8.6) 

where  AN  is  the  total  number  of  displacement  particles,  (J)  is  the 
electron  fluence,  N  is  the  number  of  GaAs  atom  per  unit  volume,  d  is 
the  cell  thickness.  Figure  8.5  shows  the  carrier  removal  rate  vs. 
electron  fluences  for  the  undoped  and  Sn-doped  GaAs  irradiated  by  one- 
MeV  electrons  as  well  as  the  results  published  recently  by  authors. 
121, 130-132 J  The  carrjer  removal  rate  shows  a  wide  range  of  variation 
from  0.5  to  5  cm"  .  The  reason  for  the  large  variation  in  the  value 
of  carrier  removal  rate  observed  by  different  researches  is  not  due  to 
the  impurity  effects  HHJ  rather  than  in  the  electron  beam  flux.  It 
is   important   to  point  out  that  the  large  carrier  removal  rate  is 
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Figure  8.5     Carrier  removal   rate  vs.  electron  fluence  for  one-MeV 
electron  irradiated  GaAs  solar  cells  and  the  published 
data  for  GaAs. 
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produced  by  large  beam   flux  while  small  beam  flux  produces  small 
:ier   removal  rate.    [111/131,132] 


cam 


8.4     Defect  Parameters  for  the  Low-Energy  Proton 
Irradiated  GaAs  Solar  Cells 


In  this  section,  we  present  the  results  of  our  studies  of  the 
deep-level  defects  in  (AlGa)As-GaAs  solar  cells  produced  by  low-energy 
proton-irradiation  for  different  proton  energies  (i.e.,  50,  100,  200, 
and  290  keV)  and  fluences  (1010,  1011,  5xl0U,  1012,  and  1013  p/cm2) , 
and  Sn-dopant  densities  (2x10  and  8x10  cm-3).  Figure  8.6  and 
table  8.1  summarize  the  defect  levels  observed  in  the  low-energy 
proton-irradiated  (AlGa)As-GaAs  cells.  The  main  hole  trap  observed  in 
these  proton  irradiated  cells  is  the  Ev+0.71eV  level,  know  as  B 
center,  tlHJ  which  was  observed  in  both  unirradiated  and  proton 
irradiated  cells.  The  F^+0.52  eV  hole  trap  is  another  important  hole 
trap  which  was  observed  in  cells  irradiated  with  proton  fluences 
equal  or  greater  than  10  p/cm  .  The  main  electron  traps  observed 
in  these  proton  irradiated  cells  are  the  Ec-0.71  eV  and  Ec-0.31  eV. 
The  Ec-0.71eV  electron  trap  is  believed  to  be  the  same  level  as  the 
Ev+0.71eV  hole  trap  since  the  energy  bandgap  for  GaAs  is  1.42  eV  at 
300°K.  In  fact,  the  observed  increase  in  the  recombination  current  in 
the  unirradiated  and  proton  irradiated  cells  is  due  to  the 
recombination  of  electron-hole  pairs  via  this  midgap  trap.  The  Ec- 
0.31  eV  electron  trap  is  believed  to  be  related  to  the  gallium  vacancy 
as  reported  by  Lang.'-11  '  Fig. 8. 7  shows  a  plot  of  defect  density 
profiles    for    the   Ec-0.31   eV,    Ec-0.6   eV,    Ev+0.4   eV,    and    Ev+0.52   eV 
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Figure  8.7     The  defect  profile  of  several    traps   in  GaAs, 
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traps  observed  in  the  200  keV  and  1013  p/cm2  proton  irradiated  cell. 
From  the  result,  it  is  noted  that  defect  density  profile  for  electron 
traps  follows  the  Gaussian  distribution,  indicating  that  these 
electron  traps  are  vacancy  related  defects.  As  for  the  hole  traps, 
their  defect  profile  is  more  resembled  to  that  of  the  background 
dopant  density  profile  in  n-GaAs  indicating  that  most  of  these  hole 
traps  may  be  related  to  the  impurity  complex  defects. 

The  main  results  and  conclusions  deduced  from  the  DLTS  analysis 
of  the  proton  irradiated  GaAs  solar  cells  are  summarized  as  follows: 
(1).  The  50  keV  proton  damages  are  stopped  at  about  0.47  um  below  the 
surface,  and  thus  most  of  the  damages  are  in  the  (AlGa)As  window 
layer  and  the  p+GaAs  emitter  region.  No  deep-level  traps  were 
observed  in  the  n-GaAs  region  indicating  that  damages  were 
confined  in  the  proton  penetrated  region.  (2).  The  200  keV  protons 
penetrate  about  1.8  um,  and  hence  create  damages  and  deep- level  traps 
in  the  p-n  junction  region  of  the  cell.  (3).  The  290  keV  protons 
penetrate  deeper  (around  2.5  um) ,  causing  damage  throughout  the  active 
region  of  the  n-GaAs  layer.  (4).  For  proton  fluences  less  than  1010 
p/cm  ,  little  or  no  change  in  deep  level  defects  was  observed  in  these 
cells  indicating  that  proton  irradiation  at  this  fluence  level 
produces  no  damage  to  the  crystal  itself.  (5).  Increasing  proton 
fluence  will  produce  more  deep-level  defects  and  higher  defect 
densities   than   increasing   the  proton  energy.    In  general,    it   is   found 

that  the  density  and  number   of  defects   will    increase  with   increasing 

12  2 

proton  fluence.      (6).   For  proton  fluence  less  than  10       p/cm    ,   the 

dominant  hole  trap  is  due  to  the  F^+0.71  eV,   and  the  dominant  electron 

trap  is  due  to  the  Ec-0.71  eV  trap.     This  midgap  trap  level  serves  as 
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an  effective  recombination  center  which  controls  the  electron  and 
hole  lifetimes  in  the  junction  space  charge  region  and  in  the  n- 
GaAs  active  layer.  (7).  For  proton  fluence  equal  to  10  p/cm  , 
significant  carrier  removal  occurs.  The  defect  spectra  for  both 
electron  and  hole  traps  observed  at  this  fluence  differ  significantly 
from  the  rest  of  proton  irradiated  samples  studied.  The  complex  and 
large  DLTS  spectra  observed  in  the  10  p/cm  irradiated  cells  clearly 
illustrate  that  considerable  damages  have  been  produced  in  these 
cells.  From  the  above  observation,  it  is  shown  that  defect  levels 
created  by  low-energy  proton  irradiation  are  indeed  quite  complicated. 
This  indicates  that  higher  order  bombardment  due  to  heavy  particles 
such  as  protons  is  important  in  the  proton   irradiated  cells. 

Important  recombination  parameters  such  as  carrier  lifetimes, 
diffusion  lengths,  and  capture  cross  sections  can  be  determined  from 
the  DLTS  data  presented  above.  This  is  discussed  as  follows: 
Minority  carrier  lifetimes  are  directly  related  to  the  defect  density 
and  capture  cross  section  in  each  trap.  The  minority  diffusion  length 
is  derived  from  the  diffusivity  and  minority  lifetimes.  Since  protons 
can  not  penetrate  throughout  the  entire  solar  cell,  it  is  expected 
that  the  minority  carrier  diffusion  length  is  nonuniform.  In  the 
proton  penetration  range,  carrier  diffusion  length  in  the  irradiated 
cells  decreases  caused  by  the  radiation  damage.  Beyond  the 
penetration  range,  minority  carrier  diffusion  length  in  the  irradiated 
cells  has  the  same  value  with  that  of  the  unirradiated  cells.  Figure 
8.8  shows  the  physical  model  for  the  irradiated  (AlGa) As-GaAs  solar 
cells.  Wysocki  et  al.  I-115-!  presented  the  effective  diffusion  length 
L    as  follow 
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L     =   {L]_  /   [(L1+LQ)exp(x1/L1)    -    (Lj-Lq) exp^x-,/!^) } 
{(L1+LQ)exp(x1/L1)    +    (L1-LQ)exp(-x1/L1) 
+  2    [(L02  -  LL2)   /  L±]}  (8.7) 

where  LQ  is  the  minority  carrier  diffsion  length  for  the  unirradiated 
cell  and  L-^  is  the  carrier  diffusion  length  of  the  irradiated  cells, 
x-^  is  the  proton  penetrate  depth.  Figure  8.9  and  table  8.2  show  the 
calculated  diffusion  length  vs.  proton  fluence  for  different  proton 
energies.  The  value  of  Lj_  is  calculated  from  the  DLTS  data.  Note 
that  L  is  independent  of  fluence  at  low  fluence  level,  and  at  higher 
fluence  level,  values  of  L  decreases  rapidly  with  increasing  proton 
fluence.  As  the  proton  energy  increases,  the  decrease  in  L  occurs  at 
a  much  lower  value  of  proton  fluence. 

8.5     Defect  Parameters  for  the  One-MeV  Electron  Irradiated 
(AlGa)As-GaAs  Solar  Cells 


In  this  section,  we  present  the  results  of  our  study  of  the  deep- 
level  defects  induced  by  one-MeV  electron  irradiation  in  GaAs  solar 
cells  subject  to  different  electron  fluences  (1014,  1015,  5xl015,  and 
1016  e/cm2),  fluxes  (2xl09  and  4xl010  e/cm2.s)  and  annealing  times 
(10,  20,  30,  and  60  min.).  Table  8.3  shows  the  trap  levels  vs. 
electron  fluences  in  one-MeV  electron  irradiated  GaAs  solar  cells. 
From  the  results,  it  is  concluded  that  (1).  One-MeV  electron 
irradiation  on  GaAs  solar  cell  will  in  general  produce  three  to  four 
electron  traps  and  one  to  five  hole  traps  depending  on  the  electron 
fluence  and  flux  as  well  as  annealing  conditions.  (2).  Defect  density 
and  the  number  of  deep  level  defects  will  increase  with  increasing 
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electron  fluence  and  flux.  (3).  Increasing  annealing  temperature  and 
annealing  time  will  reduce  the  density  and  number  of  deep  level  traps 
induced  by  electron  irradiation.  (4).  Low  temperature  thermal 
annealing  is  effective  in  reducing  the  density  of  both  electron  and 
hole  traps.  (5).  Deep  level  defects  observed  in  the  undoped  GaAs 
cells  are  somewhat  different  from  those  observed  in  the  Sn-doped  GaAs 
cells  (e.g.,EF4  =  Ec-0.31  eV  for  the  Sn-doped  GaAs,  and  EF4  =  Ec- 
0.41  eV  for  the  undoped  GaAs).  (6).  A  comparison  of  the  hole  traps  in 
the  undoped  GaAs  cells  irradiated  at  room  temperature  and  at  200°C 
revealed  that  several  shallow  level  hole  traps  were  annealed  out  for 
cells  irradiated  at  200°C  (7).  A  combined  thermal  and  recombination 
enhanced  annealing  process  may  be  used  to  eliminate  some  of  the  defect 
levels  in  the  one-MeV  electron  irradiated  GaAs  solar  cells,  and  thus 
improved  the  performance  of  GaAs  solar  cells.  Table  8.4  summarizes 
the  deep- level  defects  and  recombination  parameters  deduced  from  the 
DLTS  data  for  the  undoped  and  Sn-doped  GaAs  irradiated  by  one-MeV 
electrons  for  different  electron  fluences.  Figure  8.10  shows  the  hole 
diffusion  length  vs.  electron  fluence  for  the  one-MeV  electron 
irradiated  GaAs  solar  cells,  as  calculated  from  the  DLTS  data.  The 
detailed  results  of  the  one-MeV  electron  irradiated  (AlGa)As-GaAs 
solar  cells   is  presented   in  next. 


8.5.1      Sn-doped  GaAs   Irradiated   at   Room  Temperature  and  Post  Annealed 
at  230UC 


Figures  8.11  and  8.12  show  the  DLTS  scans  of  electron  and  hole 
traps  for  the  Sn-doped  GaAs  cells  irradiated  at  room  temperature  for 
1016  e/cm2   fluence  and  subsequently  annealed  at  230°C  for  20,    30,   and 
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60  min.,  respectively.  Three  electron  traps  with  energies  of  Ec-0.31, 
0.71,  and  0.90  eV  and  one  hole  trap  with  energy  of  Ev+0.71  eV  were 
observed  in  these  cells.  With  the  exception  of  the  EFll  (Ec-0.90  eV) 
level,  the  DLTS  results  show  that  density  of  each  trap  decreases  with 
increasing  annealing  time  following  a  230°C  annealing.  The  results 
also  show  that  shallower  trap  such  as  the  EF4  (Ec-0.31eV)  level  can  be 
more  easily  annealed  out  by  the  low  temperature  annealing  than  the 
deeper  trap  such  as  the  EFll  level.  Significant  reduction  in  the 
density  of  E^+^.lleV  hole  trap  was  also  observed  by  the  230°C  thermal 
annealing,  as  is  shown  in  Fig.  8.12. 

A  si  mi  liar  DLTS  result  was  also  obtained  for  the  Sn-doped  GaAs 

1  c  9 

cells  irradiated  with  10x  e/cm  electron  fluence.  However,  the 
density  of  each  trap  level  was  found  to  be  lower  than  those  cells 
irradiated  with  10  e/cm  fluence.  Table  8.5  summarizes  the  defect 
and  recombination  parameters  deduced  from  the  DLTS  and  C-V  data  for 
the  Sn-doped  GaAs  solar  cells  discussed  above.  Note  that  hole 
diffusion  lengths  calculated  from  the  DLTS  data  were  found  to  vary 
between  0.5  and  2.5  um,  depending  on  the  fluence  and  annealing  time 
used.  Figures  8.13  and  8.14  show  the  defect  density  vs.  annealing 
time  for  both  electron  and  hole  traps  for  the  Sn-doped  GaAs  cells 
annealed  at  230°C  for  10  to  60  min.,  respectively.  The  annealing 
rates  for  each  electron  and  hole  trap  level  are  also  listed  in 
Figs.  8.13  and  8.14. 
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8.5.2        Sn-Doped  GaAs  Irradiated  at  150  and  200°C  Cell's  Temperature 
for  Different  Fluxes  and  Fluences 


To  study  the  effects  of  varying  the  flux  and  the  cell's 
temperature  during  electron  irradiation  on  deep  level  defects, 
irradiation  on  Sn-doped  GaAs  cells  was  made  for  fluxes  of  2xl09  and 
4xl010  e/cm.s  for  a  total  fluence  of  1015  e/cm2,  and  isothermally 
annealed  at  150  and  200°C  during  electron  irradiation.  The  DLTS 
results  show  that  the  dominant  electron  trap  is  due  to  the  Ec-0.71  eV 
level  and  the  dominant  hole  trap  is  due  to  the  Ev+0.71  eV  level. 
Density  of  both  traps  was  found  to  increase  with  increasing  incident 
flux  and  decreasing  cell's  temperature.  It  is  noted  that  cells 
irradiated  at  2x10^  e/cm  s  flux  have  only  one  electron  trap  observed 
as  compared  to  three  electron  traps  for  the  cells  irradiated  at  4xl010 
e/cm^s  flux.  Our  DLTS  results  clearly  indicate  that  isothermal 
annealing  process  used  during  electron  irradiation  is  as  effective  in 
reducing  the  electron  trap  density  as  the  post  annealing  treatment. 
Table  8.6  summarizes  the  defect  and  recombination  parameters  for 
electron  and  hole  traps  deduced  from  the  DLTS  and  C-V  data  for  the  Sn- 
doped  GaAs  cells  irradiated  with  different  fluxes  and  cell's 
temperature  during   irradiation. 

8.5.3      Undoped    GaAs    Irradiated    at    200-^C    for    1014    and    1015    e/cm2 
Fluences 

In  addition  to  the  DLTS  results  shown  above  for  the  Sn-doped  GaAs 
cells,  we  have  also  performed  electron  irradiation  on  the  undoped  GaAs 
cells  at  200°C  for  electron  fluences  of  1014  and  1015  e/cm2.  Fig.8.15 
shows  the  DLTS  scans  of  electron  traps  for  cells   irradiated  with  1014 
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and  10  e/cm  fluences  and  continuously  annealed  at  200°C;  four 
electron  traps  with  energies  of  Ec-0.31f  0.41,  0.71,  and  0.90  eV  were 
observed  in  these  cells.  Note  that  both  Ec-0.13  and  0.41  eV  electron 
traps  were  not  detected  in  the  Sn-doped  GaAs  shown  in  Fig.  8.11;  while 
the  Ec-0.31  eV  electron  trap  observed  in  the  Sn-doped  GaAs  was  not 
detected  in  the  undoped  GaAs.  Another  interesting  result  for  the 
undoped  GaAs  is  that  the  dominant  electron  traps  are  due  to  the  Ec- 
0.13  and  0.41  levels,  with  trap  density  one  to  two  orders  of  magnitude 
higher  than  that  of  Ec-0.71  and  0.90  eV  electron  traps.  DLTS  scans  of 
hole  traps  for  the  undoped  GaAs  cells  showed  that  two  hole  traps  with 
energies  of  Ev+0.29  and  0.71  eV  were  observed  in  cells  irradiated  with 
10  e/cm  fluence,  and  only  one  hole  trap  with  energy  of  Ev+0.71  eV 
was  observed  in  cells  irradiated  with  10  e/cm  fluence.  Note  that 
both  electron  and  hole  trap  density  were  found  to  increase  with 
increasing  electron  fluence.  However,  the  Ev+0.29  eV  hole  trap 
observed  in  the  undoped  cell  was  not  detected  in  the  Sn-doped  GaAs 
cells.  Table  8.7  summarizes  the  defect  and  recombination  parameters 
deduced  from  the  DLTS  and  C-V  data  for  the  undoped  GaAs.  In  addition 
to  the  low  temperature  thermal  annealing  study  discussed  above,  we 
have  also  conducted  a  combined  thermal  and  recombination  enhanced   (by 

applying  a  forward  bias  current  to  the  cell;  current  density  of  1.9 

9 
A/cm     was  used   in  this  case)    annealing  experiment  on  some   of   the 

undoped   GaAs   cells    irradiated   by  one-MeV  electrons   at    200°C.      A 

comparison   of   our   DLTS    results    for    the   electron    and    hole    traps 

observed    in   these  cells   with  and   without  recombination-enhanced 

annealing    (REA)    process   indicates    that   REA  process    is   very   effective 

in  further  reducing  the  density  of  the  electron  and  hole  traps  induced 
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Table  8.7  Electron  and  hole  trap  parameters  in  one-MeV  electron 
irradiated  (AlGa)As-GaAs  solar  cells  (undoped  n-GaAs) 


<De    Nd 

(e/cra2)  (era-3) 

Electron  traps        Hole  traps 
Efc(eV)   Nt(cm"3)   Efc(eV)   Nt(aiT3) 

Irradiation 
Temp. (°C) 

1014   1.45xl015 

Ec-0.13  3.2xl013  EL+0.29 
Ec-0.41  1.3xl0^-3   EL+0.71  6.4xl012 
Ec-0.71  1.2xl012 
Ec-0.90  1.6xl012 

200 

1015   1.05xl015 

Ec-0.13  2.2xl0^4  Ey+0.29  6.9xl012 
Ec-0.41  1.3xl014  F  +0.71  2.0X1013 
Ec-0.71  7.8xl012 
Ec-0.90  9.5xl012 

200 

1014   1.5xl015 

Ec-0.14  1.8xl013  EL +0.71  3.0X1012 

Ec-0.41  8.2xl012 

Ec-0.71 

Ec-0.90  8.8X101-1- 

25 

1015    1.1X1015 

Ec-0.41  2.1xl0|4   EL+0.13   2.2xl014 

Ec-0.71  1.7x10 |3  Ev+-.29  4.0X1014 

Ec-0.90  2.8xl013  EL+0.35  8.0X1013 

Ev+0.71  6.4xl013 

25 

5xl015  9.4xl014 

Ec-0.41  4.5xl0^4  Ey+0.13  8.9xl014 

Ec-0.71  3.2xl013  EL+0.29  1.5xl015 

Ec-0.90  5.0X1013  EL+0.35  1.0xl015 

Ey+0.71  1.9xl014 

25 

1016   8.4xl014 

Ec-0.41  8.8xl0^3  EL+0.13  8.4xl014 

Ec-0.71  5.0X1013  EL+0.29  1.6xl015 

Ec-0.90  6.5xl013   EL+0.35  1.0xl015 

Ev+0.71  2.7xl014 

25 
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by  one-MeV  electrons.  Furthermore,  this  technique  appears 
particularly  effective  for  reducing  the  density  of  deeper  hole  trap 
(e.g.,  Ey  +  0.71  eV).  Thus,  it  is  obvious  that  a  combined  low 
temperature  thermal  and  recombination  enhanced  annealing  process 
should  be  a  powerful  technique  for  eliminating  the  one-MeV  electron 
radiation  induced  defects  in  the  GaAs  solar  cells. 


8.5.4     Undoped  GaAs   Irradiated  at  Room  Temperature   for  Different 
Electron  Fluences 


Study  of  deep- level  defects  in  the  undoped  GaAs  cells  irradiated 
at  room  temperature  for  four  different  electron  fluences  has  also  been 
carried  out  and  the  results  are  discussed  next.  The  DLTS  scans  of 
electron  traps  showed  that  three  electron  traps  with  energies  of  E  - 
0.41,  0.71,  and  0.90  eV  were  observed  in  the  1016  e/cm2  irradiated 
cell.  For  1014  e/cm2  fluence,  only  two  shallow  electron  traps  with 
energies  of  Ec-0.13  and  -0.41  eV  were  observed.  As  electron  fluence 
increases  to  1015  e/cm2,  the  Ec-0.41eV  electron  trap  becomes  the 
dominant    level    while    deeper    traps    such    as    Ec-0.71    and    -0.90    eV 

starting  to  emerge   in  the  DLTS   scan,   while  becomes  more  prominant  at 

IS  7 

5xl0J-'  e/cm^  fluence.      The  DLTS  scans  for   the  hole  traps  in  these 

undoped  GaAs  cells  as  a  function  of  electron  fluence  have  also  been 

made.      For    10        e/cm      electron    fluence,    only  one   hole    trap    was 

detected  with  energy  of  Ev+0.71  eV.     As  fluence  increases   to  1015 

e/cm  ,    three  additional  hole  traps  with  energies  of  Ev+0.20,    +0.29  eV, 

and  +0.35  eV  were  observed.     The  density  of  these  hole  traps  continued 

to  increase  with  increasing  electron  fluence.     Note  that  at  1016  e/cm2 

fluence,    another   hole   trap  with  energy  Ev+0.40  eV  was   also   detected. 
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Thus,  at  10  e/cm  electron  fluence,  there  are  five  hole  traps  and 
three  electron  traps  observed  in  the  undoped  GaAs.  The  density  of 
shallow  hole  traps  such  as  the  Ey+0.29  and  0.35eV  level  is  very  large 
at  10  e/cm2  fluence,  and  is  comparable  to  the  background  carrier 
density  in  the  n-GaAs  active  layer.  Since  the  deep  level  traps  such 
as  the  E  -0.71,  -0.90eV  and  Ev+0.71eV  levels  are  the  dominant 
recombination  centers  in  these  electron  irradiated  cells,  the  defect 
parameters  deduced  from  the  DLTS  data  can  be  used  to  calculate  the 
recombination  parameters  such  as  lifetimes  and  diffusion  lengths  due 
to  these  levels.  Table  8.7  summarizes  the  defect  and  recombination 
parameters  deduced  from  the  DLTS  data  for  the  undoped  GaAs  irradiated 
at  200°C  and  at  room  temperature  by  one-MeV  electrons.  The  carrier 
removal  rate  and  the  defect  density  vs.  electron  fluence  for  the 
undoped  GaAs  cells  irradiated  at  room  temperature  are  shown  in  Figs. 
8.16  and  8.17,  respectively.  Figure  8.16  shows  the  carrier  density 
and  the  density  of  electron  traps  (Ec-0.31,  0.71  and  0.90eV)  vs 
electron  fluence  for  the  undoped  GaAs  cells.  Figure  8.17  shows  the 
density  of  hole  traps  (i.e.,  Ev+0.20,  0.29,  0.35,  and  0.71eV)  vs. 
electron  fluence  for  the  undoped  GaAs  cells  discussed  above.  It  is 
interesting  to  note  that  density  of  both  electron  and  hole  traps 
increases  rapidly  with  increasing  electron  fluence  initially,  and 
becomes  saturated  at  10x  e/cm  fluence  when  the  defect  densi  ty  equals 
to  the  background  dopant  density  in  the  undoped  n-GaAs  active  layer. 

A  comparison  of  the  DLTS  results  for  the  undoped  GaAs  cells 
irradiated  at  200°C  and  at  room  temperature  revealed  that  several 
shallow  hole  traps  such  as  Ev+0.20,  0.29,  0.35,  and  0.71eV  observed  in 
the  room  temperature  irradiated  GaAs  cells  were  not  detected  in  the 
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0p(cm' 


Figure     8.16     Density  of  electron  traps  and  background  electron 
density  vs.   electron  fluence  for  the  undoped  GaAs 
solar  cells  irradiated  by  one-MeV  electrons  at  room 
temperature. 
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Figure  8.17     Density  of  hole  traps  vs.   electron  fluence  for  the 
undoped  GaAs   irradiated  by  one-MeV  electrons   at 
room  temperature. 
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cells  which  were  continuously  annealed  at  200°C  during  electron 
irradiation.  The  results  suggested  that  these  hole  traps  were 
effectively  annealed  out  when  the  cell's  temperature  was  raised  to 
200°C  during  electron  irradiation.  However,  no  significant  different 
was  detected   in   the  DLTS   spectrum  for  the  electron  traps   in  both 


8.6     Conclusions  and  Summaries 

In  this  study,  DLTS  and  C-V  measurements  analysis  have  been  made 
on  both  one-MeV  electron  and  low-energy  proton  irradiated  GaAs  solar 
cells. 

For  low  energy  proton  irradiated  GaAs  cells,  six  electron  traps 
with  energy  of  Ec-0.14,  Ec-0.20,  Ec-0.31,  Ec-0.52,  Ec-0.60,  and  Ec- 
0.71  eV  and  seven  hole  traps  with  energy  of  F^+0.17,  Ey+0.27,  F^+0.29, 
E  +0.44,  Ey+0.52,  Ev+0.57,  and  Ev+0.71  eV  were  observed.  The  main 
results  are  summarized  as  follows:  (1)  No  new  trap  was  detected  for 
proton  energy  less  than  50  keV  and  fluence  less  than  1010  p/cm  ;  (2) 
increasing  proton  fluence  has  more  influence  on  trap  density  and 
the  number  of  trap  levels  than  increasing  the  proton  energy;  (3)  trap 
density  and  the  number  of  deep  level  traps  will  increase  with 
increasing  proton  fluence;  (4)  for  proton  fluences  higher  than  10 
p/cm2,  significant  carrier  removal  occurs  in  the  junction  space  charge 
region  of  the  solar  cell;  (5)  most  of  the  measured  electron  trap 
density  profile  has  a  Gaussian  distribution. 

For  one-MeV  electron  irradiated  GaAs  solar  cells,  five  electron 
traps   with  energies   of   Ec-0.14,    Ec-0.20,    Ec-0.31,    Ec-0.71,    and   Ec~0.90 
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eV,  and  four  hole  traps  with  energies  of  ev+0-13,  Ey+0.29,  Ev+0.35, 
and  E  +0.71  eV  were  observed.  The  main  results  and  conclusions  are 
summarized  as  follows:  (1)  The  Sn-doped  GaAs  cells  are  more  radiation 
hard  than  the  undoped  GaAs  cells;  (2)  high  temperature  irradiation 
(200°C)  produces  less  defect  density  and  few  defect  levels  than  room 
temperature  irradiation;  (3)  defect  density  and  the  number  of  trap 
level  will  increase  with  increasing  electron  fluence  and  flux;  (4) 
increasing  annealing  temperature  and  time  will  reduce  the  density 
and  the  number  of  deep  level  traps;  (5)  low  temperature  thermal 
annealing  is  effective  in  reducing  the  density  of  both  electron  and 
hole  traps. 

A  comparison  of  the  deep-level  defects  produced  by  the  one-MeV 
electron  irradiation  and  the  low-energy  proton  irradiation  reveals 
that  there  is  a  significant  difference  in  defect  profile  as  well  as 
types  of  defects  observed  in  these  GaAs  solar  cells.  Th  defect 
profiles  are  much  more  complicated  for  the  proton  irradiated  cells 
than  the  one-MeV  electron  irradiated  cells. 


IX.   CONCLUSIONS 

A  detailed  theoretical  and  experimental  study  of  the  grown-in 
deep-level  defects  and  radiation-induced  defects  in  GaAs  has  been 
carried  out  in  this  research.  The  main  contributions  from  this 
research  include:  (i)  development  of  a  new  defect  model  for  predicting 
defect  density  vs  growth  parameters  in  GaAs,  (ii)  providing  a  new 
interpretation  of  the  EL2  electron  trap  in  GaAs,  (iii)  determining  the 
energy  level  and  trap  density  for  a  large  number  of  deep-level  traps 
in  GaAs  grown  by  LBC,  VPE,  LPE,  and  MOCVD  techniques  under  growth  and 
annealing  conditions,  (iv)  development  of  the  nonexponential  DLTS 
response  for  determining  the  potential  well  of  each  trap  in  GaAs.  The 
research  described  in  this  dissertation  has  led  to  the  following 
conclusions: 

(1)  The  theoretical  model  for  the  nonexponential  capacitance 
transients  due  to  electric  field  dependent  emission  rate  of  trapped 
charge  has  been  developed.  Based  on  the  nonexponential  capacitance 
transient,  theoretical  calculations  of  DLTS  response  for  the  deep 
level  traps  in  GaAs  were  made  using  five  different  potential  wells; 
namely,  coulombic  well,  Dirac  well,  square  well,  polarization  well, 
and  dipole  well.  A  comparison  of  the  theoretical  calculation  of  the 
nonexponential  DLTS  response  with  the  DLTS  data  for  each  trap  level 
would  allow  us  to  determine  the  potential  well  of  the  trap  involved. 
This  method  has  been  applied  to  identify  the  physical  origins  of  the 
EL2  electron  trap   in  GaAs. 
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(2)  Theoretical  model  for  predicting  native  point  defects  in 
GaAs  has  been  developed.  During  crystal  growth,  the  chemical 
thermodynamic  principle  is  used  to  derived  the  density  of  vacancy  and 
interstitial  defects  under  thermal  equilibrium  condition.  After 
crystal  growth,  thermal  kinetic  is  employed  to  predict  the  ant i site 
defect.  Conclusions  are  listed  as  follows:  (i)  high  purity  material 
can  be  grown  for  the  low  arsenic  pressure  case  under  optimum  cooling 
condition,  (ii)  GaAs  grown  under  higher  arsenic  pressure  condition 
[e.g.,  LBC,  VPE  and  MOCVD  techniques]  will  produce  more  native  point 
defects  than  under  lower  arsenic  pressure  condition  [e.g.,  LPE 
technique],  (iii)  Arsenic  antisite  defect  only  observed  in  GaAs  grown 
under  As-rich  or  high  arsenic  pressure  condition;  this  defect  can  not 
be  produced  under  low  arsenic  pressure  and  Ga-rich  condition. 

(3)  A  new  model  for  the  EL2  electron  trap  in  GaAs  is  presented 
in  this  study.  EL2  electron  trap  is  attributed  to  two  different 
native  point  defects.  One  is  identified  as  the  EL2a  (Ec-0.83eV) 
electron  trap  and  the  other  is  denoted  as  the  EL2b  (Ec-0.76eV) 
electron  trap.  The  physical  origin  for  the  EL2a  level  is  attributed 
to  arsenic  antisite  (AsGa++)  defect,  whereas,  the  physical  origin  for 
the  EL2b  level  is  attributed  to  arsenic  anti site-plus-arsenic  vacancy 
(Asr  VA  ).  Based  on  this  model,  relationship  between  the  density  of 
EL2a  and  EL2b  trap  levels  vs  [As]/[Ga]  mole  fraction  ratio  in  the 
MOCVD  and  VPE  grown  GaAs  is  established.  The  result  shows  that 
density  of  EL2a  trap  is  proportional  to  the  mole  fraction  ratio  of 
(r-1)1/2;  while  the  density  of  EL2b  trap  is  proportional  to  the  mole 
fraction  ratio  of  (r-1)1/4,  where  r  =  [As]/[Ga].  This  conclusion  is 
supported  by  experimental  data  for  the  MOCVD  and  VPE  grown  GaAs. 
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(4)  From  the  results  of  our  study  of  the  grown-in  defects  in  the 
LEC,  VPE,  LPE,  and  MOCVD  GaAs  samples,  it  is  clear  that  grown-in 
defects  and  the  background  concentration  are  sensitive  to  the  growth 
condictions,  stoichi ometry,  and  substrate  orientations.  The  main 
conclusions  are  listed  as  follows:  (i)  EL2a  electron  trap  is 
observed  in  VPE  GaAs  layers  and  in  LEC  grown  GaAs  annealed  at  500°C  in 
H2  ambient,  (ii)  EL2b  electron  trap  was  observed  in  unannealed  LEC 
GaAs  layers,  and  it's  density  can  be  reduced  after  thermal  annealing, 
(iii)  Ey+Q.QQeV  hole  trap  has  been  observed  in  the  LEC  p-GaAs  samples. 
The  physical  origin  of  this  hole  trap  is  believed  to  be  due  to  gallium 
antisite  (Ga,  )  defect,  (iv)  reducing  the  temperature  cooling  rate 
during  crystal  growth  is  beneficial  for  reducing  the  density  of  grown- 
in  defects  in  the  LPE  GaAs  layers.  These  experimental  results  are 
supported  by  the  theoretical  model  of  native  point  defects  and  EL2 
electron  trap  in  GaAs  described  in  this  work. 

(5)  For  low  energy  proton  (i.e.,  E  =  50  to  290  keV)  irradiated 
GaAs  solar  cells,  six  electron  traps  with  energy  of  Ec-0.14,  Ec-0.20, 
Ec-0.31,  Ec-0.52,  Ec-0.60,  and  Ec-0.71eV  and  seven  hole  traps  with 
energy  of  Ey+0.17,  Ev+0.27,  Ey+0.29,  Ey+0.44,  Ev+0.52,  Ev+0.57,  and 
EL>0.71eV  were  observed.  For  one-MeV  electron  irradiated  GaAs  solar 
cells,  five  electron  traps  with  energies  of  Ec-0.14,  Ec-0.20,  Ec-0.31, 
Ec-0.71,  and  Ec-0.90eV,  and  four  hole  traps  with  energies  of  Ev+0.13, 
Ey+0.29,  Ey+0.35,  and  E^+QJIeV  were  observed.  The  radiation  induced 
deep-level  defects  are  believed  to  be  Frenkel  type  defects  such  as 
vacancies,  interstitials,  antisite  defects  and  their  complexes.  Our 
study  reveals  that  low  temperature  thermal  annealing  or  periodic 
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annealing  process  is  effective  in  reducing  or  eliminating  these 
radiation  induced  defects  in  GaAs  solar  cells.  The  shallow  traps  are 
much  easier  to  anneal  out  than  the  deep- level  traps. 
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